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HAUPTAUFSÄTZE 


Über eine neue T'heorie der freien Turbulenz. 
Von H. Reichardt in Göttingen. 
(Aus dem Kaiser-Wilhelm-Institut für Strömungsforschung.) 


Impulsmessungen in turbulenten Vermischungsgebieten zeigen, daß zwischen 
den turbulenten und den molekularen Ausgleichsvorgängen eine Analogie be- 
steht, die eine sinngemäfße Anwendung der Differentialgleichung der Wärme- 
leitung auf die turbulente Impulsausbreitung ermöglicht. Diese Tatsache 
bildet die (rundlage einer neuen phänomenologischen Theorie der freien 
turbulenten Strömungen. 


1. Leistung und Grenzen der bisherigen Theorie des Mischungsweges. 


Während die Molekulartheorie von bestimmten Voraussetzungen ausgehen kann, weil 
(die Elemente ihrer Statistik, die Moleküle, vorgegeben sind, ist die Turbulenztheorie in der 
schwierigen Lage, daß sie zunächst keine sichere Grundlage hat, auf der sie eine Statistik 
aufbauen kann. Denn die Elementarteile der Turbulenz, die Wirbelballen, sind sehr wechsel- 
volle Gebilde, die sich in schwer erkennbarer Weise gegenseitig beeinflussen. Man kann da- 
her die Ergebnisse oder die Methoden der Molekularphysik nieht ohne weiteres auf die Tur- 
bulenzmechanik übertragen. 

Trotz beträchtlicher Unterschiede haben aber die molekularen und die turbulenten Vor- 
gänge doch wesentliche Züge gemeinsam. Beispielsweise besteht eine Ähnlichkeit in der 
Schubübertragung. In beiden Fällen kommt die Schubspannung dadurch zustande, daß Ele- 
mentarteile mit einem bestimmten Impulsinhalt der Hauptgeschwindigkeit durch eine Quer- 
bewegung um eine „freie Weglänge“ in eine neue Umgebung mit anderer Hauptgeschwindig- 
keit gelangen und dort durch Stoß bzw. Vermischung ihren Impulsüberschuß an die benach- 
barten Teile abgeben oder sie durch ihren kleineren Impuls bremsen. 

Es ist daher naheliegend, den molekularen Begriff der freien Weglänge sinngemäß auf 
die turbulente Austauschbewegung zu übertragen und diese Größe zur Darstellung der turbu- 
lenten Scheinreibung zu verwenden. In der Tat ist aus solchen Überlegungen heraus ein der 
[freien Weglänge entsprechender Begriff, der sogenannte Prandtlsche Mischungsweg, für die 
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Turbulenz definiert worden'). Dieser Prandtlsche Mischungsweg hat erstmalig die rechne- 
rische Behandlung turbulenter Strömungen ermöglicht. 

Auf Grund der Vorstellung, daß ein Wirbelballen einen gewissen Weg !’ ohne Ver- 
mischung mit seiner Umgebung zurücklegt, erhält Prandtl für die Schwankungsgeschwin- 
digkeiten «’ (in der Hauptströmungsrichtung) und »’ (in der Querrichtung) die Ausdrücke 


,‚dü 
dy dy 
und somit für die turbulente Scheinreibung = ow’»’ den Ansatz 
ydy u). 


Hier ist @ die mittlere Geschwindigkeit in der Hauptströmungsrichtung =, y die Querrichtung, 
in der der Austausch erfolgt, o die Dichte des Mediums und ! der aus der Mittelwertsbildung 
von und hervorgegangene Begriff des Mischungsweges. 

Um die Berechnung der mittleren Geschwindigkeit einer turbulenten Strömung durch- 
führen zu können, bedarf es noch geeigneter Aussagen für die Größe I. Hier war es nötig, 
sich auf die Erfahrungstatsachen zu stützen, die den Weg zu brauchbaren Ansätzen lieferten. 

Die Theorie des Mischungsweges ist seit ihrer ersten Formulierung im Jahre 1925 in 
zahlreichen Fällen angewendet worden, was schon rein äußerlich in einer umfangreichen 
Literatur zum Ausdruck kommt. Diese Untersuchungen beziehen sich teils auf die wand- 
nahen Strömungen, hauptsächlich aber auf die turbulenten Ausbreitungsvorgänge ohne Be- 
hinderung durch Wände. Diese Vorgänge der sogenannten freien Turbulenz sind nämlich 
wegen ihrer relativen Einfachheit noch verhältnismäßig leicht zu behandeln. 

Die bisherigen Leistungen der Mischungswegtheorie lassen sich etwa wie folgt kurz um- 
reißen: In zahlreichen Fällen wurde die mittlere Strömungsbewegung näherungsweise be- 
rechnet. Das in Reibungsschichten auftretende logarithmische Geschwindigkeitsprofil konnte 
erklärt werden. Auf dem Gebiete der freien Turbulenz ist es gelungen, die beobachteten 
Gesetzmäßigkeiten der Ausbreitung der Vermischungsgebiete in der Hauptströmungsrichtung 
in den Hauptzügen richtig wiederzugeben. 

Die bisherigen Ergebnisse sind aber nicht derart verbindlich, daß sie nicht auch anders 
erklärt werden könnten. Man kann daher nicht eigentlich sagen, der \ Mischungswegbegriff 
habe eine Bestätigung durch die Erfahrung gefunden und es sei durch ihn ein wesentlicher 
Einblick in die Turbulenzstruktur erfolgt. Im einzelnen ergeben sich hier folgende Einwände. 

Die Untersuchungen beziehen sich ausschließlich auf die Berechnung der mittleren 
Geschwindigkeiten und deren Vergleich mit den gemessenen Mittelwerten. Diese Mittelwerts- 
untersuchungen sind aber allein nicht ausreichend für die Bestimmung einer Größe, die für 
die Struktur der Turbulenz kennzeichnend sein soll. Eine Klärung des Mischungsweg- 
begriffes wird daher wesentlich von Schwankungsuntersuchungen abhängen, die den Zu- 
sammenhang zwischen den Schwankungsgrößen und den Mittelwerten zum Gegenstand haben. 

Aber auch die bisher durchgeführten Mittelwertsuntersuchungen standen von vornherein 
im Zeichen einer grundsätzlichen Schwierigkeit. Der ursprüngliche Ansatz (1) führt nämlich 
zu Geschwindigkeitsprofilen, die verglichen mit den gemessenen Profilen in den Maxima 
wesentlich zu spitz sind (vor allem zeigen Geschwindigkeitsgradient und zweite Ableitung, 
daß der Verlauf der gerechneten Profile nicht richtig sein kann). Für eine physikalisch be- 
friedigende Darstellung turbulenter Profile wurde deshalb von Prandtl schon in seiner ersten 
Veröffentlichung?) ein komplizierterer Ausdruck vorgeschlagen, der die Form hat 


Durch diesen Ausdruck wird außer dem eigentlichen Mischungsweg I! noch eine zweite 
Mischungslänge !, definiert. Da der erweiterte Ansatz zwei zur freien Verfügung stehende 
Größen enthält, kann sich die Theorie den Experimenten besser anpassen. Tatsächlich ist 
es in einzelnen Fällen gelungen (wenn auch unter großem mathematischem Aufwand), mit dem 
Ansatz (2) wirklichkeitsgetreue Profile zu berechnen. 

Physikalisch betrachtet ist allerdings das Verfahren mit zwei freien, den jeweiligen Be- 
dürfnissen angepaßten Größen zu arbeiten, wenig befriedigend. Ist es schon nicht leicht, die 
Bedeutung des Weges I klar zu erkennen, so erscheint die aus weiteren qualitativen Über- 
legungen eingeführte Länge I, noch schwerer deutbar. Es drängt sich daher die Frage auf, 
ob man auf dem bisherigen Wege physikalisch überhaupt weiterkommen kann. 


1) L. Prandtl: Bericht über Untersuchungen zur ausgebildeten Turbulenz. Z. angew. Math. Mech. Bd. 5 
(1925), S. 136. 
2) Siehe Fußnote 1, S. 2. 
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2. Die Zweckmäßigkeit der induktiven Methode. 

Als eine wesentliche Ursache für die derzeitige Lage der Theorie erscheint dem Ver- 
fasser die bisherige Überbewertung der deduktiven Betrachtungsweise in der Turbulenz- 
forschung. Der theoretische Hydrodynamiker ist gewohnt, von den Bewegungsgleichungen 
als einem physikalisch gesicherten Fundament auszugehen. Für große Teile der Hydrodynamik 
ist dies sicher berechtigt. Auf vielen Gebieten ist die Theorie bestimmt richtig, wenn sie 
sich mathematisch durchführen läßt. Bei der Lösung schwieriger Aufgaben, die nur nähe- 
rungsweise zu bewältigen sind, braucht man zwar häufig experimentelle Daten. Diese sind 
dann aber doch organisch in eine an sich zuverlässige Theorie eingebaut. 

Ganz anders dagegen liegen die Verhältnisse in der Turbulenzforschung. Hier hat man 
so gut wie keine theoretischen Grundlagen. Diese Grundlagen sollen ja erst geschaffen 
werden. Man ist daher im wesentlichen auf die Aussagen des Experimentes angewiesen. 
Trotzdem besteht auch hier die Tendenz, die Theorie auf Grund allzemeiner Annahmen deduktiv 
zu entwickeln, was auf das Erproben hypothetischer Ansätze hinausläuft. Diese Arbeit des 
Anpassens eines Ansatzes an die experimentellen Tatsachen ist aber dadurch erschwert, daß 
zwischen dem Ansatz und dem Experiment die ganze, zumeist komplizierte Rechnung liegt, 
die bei jeder Abänderung der Voraussetzung wieder von neuem durchgeführt werden muß. 

Zweckmäßiger dagegen erscheint für eine aus der Erfahrung schöpfende Turbulenz- 
forschung das induktive Verfahren einer unmittelbaren Fundierung der Theorie auf den Ver- 
suchsergebnissen. Wenn es gelingt, die tatsächlichen Vorgänge mathematisch darzustellen, 
so erhält man durch Umkehrung des bisherigen Rechnungsganges bestimmt den richtigen 
„Ansatz“. 

Natürlich hat auch das induktive Verfahren seine Schwierigkeiten. Denn man muß 
durch die streuenden Meßpunkte die „richtigen“ Kurven legen. Dabei sind die Meßfehler 
sowie die Einflüsse überlagerter Effekte höherer Ordnung (die in einem noch wenig erforschten 
Gebiet zunächst gar nicht interessieren) richtig abzuschätzen. Diese Aufgabe, innerhalb des 
Rahmens einer unvermeidbaren Unsicherheit die wesentliche Gesetzmäßigkeit aus den ver- 
schiedensten Meßergebnissen herauszufinden, ist aber verhältnismäßig einfach zu lösen, wenn 
es gelingt, durch eine Gesamtschau der Erscheinungen zu einer brauchbaren Hypothese zu 
gelangen. 

Man kann also auch bei rein phänomenologischem Vorgehen nicht ganz auf spekulative Betrach- 
tungen verzichten. Gegenüber der deduktiven Methode besteht aber der grundsätzliche Unterschied, 
daß sich die idealisierenden Annahmen direkt auf die Meßergebnisse beziehen. Dadurch wird 
eine unmittelbare Kontrolle von Hypothesen möglich. 


3. Über die Möglichkeit eines Vergleiches zwischen den molekularen und den 
turbulenten Ausgleichsvorgängen. 


Bei den Strömungen der freien Turbulenz kann man nun unschwer eine allgemeine 
Hypothese aufstellen. Da es sich hier um statistische Ausgleichsvorgänge handelt, ist ein 
Vergleich der turbulenten Vorgänge mit den molekularen Erscheinungen der Diffusion und 
der Wärmeleitung sehr naheliegend. Wenn auch die molekularen und die molaren Prozesse 
in den Einzelheiten ganz verschieden verlaufen, so ist es doch andererseits recht wahrschein- 
lich, daß sich auch bei den turbulenten Vorgängen im Endeffekt rein statistische Verteilungen 
ergeben. In der Tat haben die in zahlreichen Fällen gemessenen Geschwindigkeitsverteilungen 
in Freistrahlen und in Windschattengebieten eine große Ähnlichkeit mit der Gaußschen 
Fehlerkurve. Die Fehlerkurve gibt vor allem den Verlauf der Profile in der Nähe des 
Maximums (der mit dem Ansatz (1) nicht richtig herauskommt) sehr befriedigend wieder’). 

Es ist daher die Vermutung naheliegend, daß Impulsunterschiede*) in einer freien turbu- 
lenten Strömung in ähnlicher Weise ausgeglichen werden wie Konzentrationsunterschiede oder 
wie Wärme- bzw. Temperaturdifferenzen. Beispielsweise sind die bereits vor zehn Jahren von 


3) An den äußersten Rändern kann die Fehlerkurve natürlich nicht mehr gelten, da die Potentialbewegung des 
Anßenraumes der turbulenten Ausbreitung eine endliche Grenze setzt. 

Die von Förthmann im ebenen Freistrahl gemessenen Geschwindigkeitsverteilungen wurden kürzlich von 
Okaya und Hasegawa durch die Gleiehung dargestellt 


42 


wobei /=y/x und k=konst. bedeutet. Die Verfasser diskutieren die Verträglichkeit dieser Geschwindigkeitsformel 
ınit dem Mischungswegansatz (2) (Proc. Physieo-Mathematical Soe. Japan, 3. Ser., Bd. 22, 1940). 


4) Physikalisch wirksam ist nicht die Geschwindigkeit als solche, sondern der Impuls gu?. Wenn 0 u? durch eine 
Feblerkurve dargestellt ist, so ist die Verteilung von u ebenfalls eine Exponentialfunktion. — In Windschatten- 
gebieten ist der Impulsverlust näherungsweise —2U u, (U =konstante Hauptbewegung, u = Geschwindigkeitsverlust). 
Hier sind also (bei nicht zu großen Untergeschwindigkeiten «,) Impulsverlust und Geschwindigkeitsverlust einander 
proportional. 
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H. Muttray’) gemessenen Impulskurven, die das allmähliche Abklingen einer Impulsdelle 
hinter einem Widerstandskörper zeigen, mit den bekannten abklingenden Kurven der Wärme- 
leitungstheorie zum Verwechseln ähnlich. (Man braucht bloß den Abstand x von der Störungs- 
stelle mit dem Zeitparameter der Temperaturkurven in geeigneter Weise zu vertauschen.) 

Die Hypothese, daß die turbulente Impulsverteilung der Hauptbewegung durch eine 
partielle Differentialgleichung nach Art der Wärmeleitungsgleichung darstellbar sei (wobei der 
Örtsparameter & an die Stelle des Zeitparameters tritt), wird also einem unbefangenen Be- 
trachter der Muttrayschen Meßergebnisse sehr nahegelegt. Wir setzen deshalb 


_ 


Wenn nun die Wärmeleitungsgleichung auf die Vorgänge der freien Turbulenz bislang 
noch nicht angewendet worden ist (dem Verfasser ist jedenfalls ein solcher Versuch nicht 
bekannt), so müssen andererseits gewisse Hemmungen gegen eine derartige Lösung bestanden 
haben. Tatsächlich erweckt ja auch die Analogiebetrachtung das grundsätzliche Bedenken, 
daß eine Größe, der der Wärmeleitfähigkeit entsprieht, im Impulsfall nicht vorgegeben ist. 
Die physikalische Bedeutung der Größe A der Gl. (3) ist daher (zunächst wenigstens) nicht 
klar ersichtlich. 

Dieser Einwand ist aber nicht entscheidend, denn jede neue Theorie bringt zunächst 
neue Unbekannte. Es muß daher ein anderer Grund gewesen sein, der einer molekularen 
Analogiebetrachtung entgegenstand. Anscheinend waren es die bisher vorliegenden Unter- 
suchungen der sogenannten freien Strahlgrenze (wo ein Freistrahl mit konstanter Geschwindig- 
keit mit ruhender Flüssigkeit in Berührung tritt), die einer solchen Betrachtungsweise 
widersprachen. 

Die ebene Berührung zweier Strahlen mit verschiedener Geschwindigkeit entspricht im 
Wärmefall der ebenen Berührung zweier Körper mit verschiedener „Temperatur. Die Theorie 
der Wärmeleitung führt bekanntlich zu Temperaturverteilungen im Übergangsgebiet, die durch 
Fehlerintegralkurven dargestellt werden. Wenn die Körper gleiche Leitfähigkeit haben, hat 
die Berührungsebene stets die mittlere Temperatur, und die beiden Äste der Fehlerintegral- 
kurve liegen (mit vertauschtem Vorzeichen) „symmetrisch“ zur Grenzfläche. Wenn also die 
Analogie zwischen den molekularen und den turbulenten Ausgleichsvorgängen wirklich brauch- 
bar ist, so müssen die Impulskurven im Vermischungsgebiet zweier Strahlen Fehlerintegral- 
kurven sein, die ebenfalls (in dem- angegebenen Sinn) symmetrisch zur Grenzfläche liegen. 


Eine fehlende Symmetrie könnte natürlich durch die Annahme verschiedener „Impulsleitfähigkeiten“ 
zu beiden Seiten der Grenze erklärt werden. Ein solcher Ausweg ist aber nicht erlaubt. Denn unsere 
Analogiebetrachtung, die ja nur das Wesentliche des turbulenten Ausgleichsvorganges erfassen soll, 
ist überhaupt nur sinnvoll bzw. kann nur durchgeführt werden, wenn man annimmt, daß das Impuls- 
ausbreitungsvermögen vom Impulswert selbst bzw. von % (wenigstens in erster Näherung) unabhängig 
ist. Daß die turbulenten Vorgänge in zweiter Näherung sehr wahrscheinlich ganz anders und erheblich 
komplizierter sind als die entsprechenden molekularen Prozesse, daß also ein Vergleich in zweiter 
Näherung recht problematisch wäre, braucht nicht näher erörtert zu werden. 


Die bisherigen Untersuchungen an der freien Strahlgrenze beziehen sich auf den am 
leichtesten zu verwirklichenden und praktisch wichtigen Grenzfall, wo die «-Geschwindigkeit 
auf der einen Seite verschwindet, wo also ein Strahl konstanter Geschwindigkeit mit ruhen- 
der Flüssigkeit in Berührung kommt. Hier müßte man also auf der Grenzfläche den halben 
Impuls der Maximalgeschwindigkeit erwarten, und die Impulsverteilungen diesseits und jen- 
seits der Grenze müßten einander genau entsprechen. In Wirklichkeit sind aber bisher nur 
deutlich unsymmetrische Verteilungen des Gesamtdrucks (der praktisch dem Impuls entspricht) 
gemessen worden ®). Der Gesamtdruck der Grenzfläche wurde stets kleiner als der halbe 
Maximaldruck festgestellt, und die gemessenen Verteilungen verlaufen auf der Seite der 
ruhenden Flüssigkeit flacher und in den Strahl hinein auffallend steiler. 

Die ausgesprochene Unsymmetrie der gemessenen Impulskurven mag auch plausibel 
erschienen sein, weil aus dem Gebiet der Ruhe ständig neue Flüssigkeit angesaugt wird, von 
der anderen Seite dagegen nicht. Es liegen also auf beiden Seiten der Grenze zweifellos 
verschiedene hydrodynamische Verhältnisse vor. Dazu kommt noch, daß man die festgestellte 
Unsymmetrie des Gesamtdruckverlaufes durch den Mischungswegansatz (1) recht gut erklären 
konnte’). Es ist daher sehr wohl verständlich, daß ein Vergleich zwischen den turbulenten 


5) H. Muttr ay: Die experimentellen Tatsachen des Widerstandes ohne Auftrieb. Handb. d. Experimentalphysik 
von Wien-Harms, IV, 2, S. 325 (1932). 

6) 5. erh "der Aerodynamischen Versuchsanstalt Göttingen, 2. Lieferung (1923) S.69. Ferner: G. Cordes, 
Ing.-Arch. Bd. 8 (1937) S. 245; v. Bohl, Ing.-Arch. Bd. 11 (1940) S. 295. (Die Meßpunkte von Cordes weichen nur wenig 
von der Fehlerintegralkurve ab.) 

?), W. Tollmien: Berechnung turbulenter Ausbreitungsvorgänge. Z. angew. Math. Mech. Bd. 6 (19:6), S. 468. 
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und den molekularen Ausgleichsvorgängen gar nicht versucht worden ist, zumal das Studium 
der Strahlgrenze (deren Verhalten für die Konstruktion von Windkanälen von erheblichem 
Interesse war) zu den ersten Untersuchungen auf dem Gebiete der freien Turbulenz gehörte. 

Diese Argumente gegen die Anwendung einer molekularen Analogiebetrachtung bestehen 
aber doch nur scheinbar. Zwar spielt die Menge der einseitig aus der Ruhe angesaugten 
Flüssigkeit eine wesentliche Rolle. Aber der Querpuls des einseitigen Zustromes ist dennoch 
ganz geringfügig. Vom Standpunkt des Impulsausgleiches (und nur darum handelt es sich 
bei der Analogiebetrachtung) ist daher die erwähnte Unsymmetrie vernachlässigbar gering. 

Aber sprechen nicht doch die vorliegenden Messungen gegen eine Übertragung der 
molekularen Betrachtungsweise? Auch dies ist ein Trugschluß. Die Messungen sind nämlich 
gar nicht an wirklich ebenen Strahlgrenzen ausgeführt worden. Die alten Göttinger Messungen 
beziehen sich auf den runden Freistrahl und die ursprünglich rechteckigen Strahlen von 
Cordes und von v. Bohl dürften in einer gewissen Entfernung von der Düse auch nicht 
mehr ganz eben gewesen sein. Beim runden Freistrahl ist es aber ziemlich selbstverständlich, 
daß die Impulskurve nach außen hin flacher und in den Strahl hinein steiler verläuft. (Das 
entsprechende Wärmeleitungsproblem (erwärmter Zylinder in einer kälteren Umgebung) ist 
ja auch unsymmetrisch.) 

Die Entscheidung der Frage, ob eine Anwendung der molekularen Analogie auf die 
. Vorgänge der freien Turbulenz grundsätzlich möglich ist oder nicht, hängt also von Gesamt- 
druckmessungen an einer wirklich ebenen Strahlgrenze ab. Der Verfasser hat daher 
solche Untersuchungen angestellt. Durch entsprechende Führung des Strahles war dafür 
gesorgt, daß auch noch in größerer Entfernung von der Düse praktisch ebene Verhältnisse 
vorlagen. 

Die unter diesen Bedingungen durchgeführten Untersuchungen zeigen nun, daß die 
Impulsverteilungen tatsächlich in erster Näherung symmetrisch sind und sich überraschend 
gut durch Fehlerintegralkurven darstellen lassen. Diese Tatsache wurde durch zahlreiche sorg- 
fältige Messungen in den verschiedensten Abständen von der Düse nachgewiesen. 


In Bild 1 ist eine solche Meßreihe 
dargestellt. Die schwarzen Me&ßpunkte 


bedeuten die in den einzelnen Quer- 
abständen y bestimmten Gesamtdruck- 

werte für einen vorgegebenen x-Ab- 

stand von der Düse (diese Werte Taen 


sind auf den maximalen Gesamtdruck 
(max der Düsenströmung bezogen). 
Die ausgezogene Linie ist die Fehler — - 
integralkurve, die auf der Grenze 


y=0 durch den Punkt geht, 


max 
die also mit ihren beiden Ästen (in 
dem angegebenen Sinn) genau sym- 
metrisch zur Grenzfläche liegt. Die 
Meßpunkte liegen auf oder in un- ’ 
mittelbarer Nähe der Fehlerintegral- Perg; 
kurve. Bei genauem Betrachten er- Bild 1. Gesamtdruckverteilung in der Vermischungszone der ebenen 
kennt man, daß eine schwache Un- a ui 
symmetrie zweifellos vorhanden ist. Diese Abweichung ist aber so geringfügig, daß sie nur 
bei sehr guter Meßgenauigkeit (die durch Mittelwertbildung über längere Zeiträume erzielt 
wurde) deutlich in Erscheinung tritt. 

Wegen der Bedeutung dieser Frage wurde für einen x-Abstand auch der Differential- 
quotient des Impulsverlaufes aus der Druckdifferenz zweier nahe benachbarter Stauröhrchen 
bestimmt. Diese Differentiationsmessungen sind durch die hellen Kreise dargestellt (wobei 
die Einzelwerte wieder auf den Größtwert bezogen sind). Bei genau symmetrischen Ver- 
hältnissen müßten die Meßpunkte auf der in Bild 1 eingezeichneten Fehlerkurve liegen, 
die aus der Differentiation der ausgezogenen Fehlerintegralkurve hervorgeht. In Wirklichkeit 
sind die Meßpunkte etwas nach links verschoben. Die Neigung der Gesamtdruckkurve ist 
also auf der Außenseite tatsächlich etwas flacher als auf der Strahlseite. 

Die Abweichungen von der Symmetrie sind aber derart gering, daß man sie für eine 
Theorie der ersten Nährung mit gutem Gewissen vernachlässigen kann ®). Damit ist das 


8) Allzu große Genauigkeitsanforderungen zu stellen, hat schon deswegen keinen Sinn, weil bei starker Turbulenz 
(besonders in den Randgebieten) mit dem Auftreten systematischer Meßfehler zu rechnen ist. Außerdem ist der Unter- 
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wesentliche Argument gegen die Anwendung der molekularen Analogiebetrachtung auf die 
freie Turbulenz beseitigt. 

Die sinngemäße Übertragung der Theorie der Wärmeleitung auf die turbulenten Aus- 
breitungsvorgänge erhält noch dadurch ihre besondere Rechtfertigung, daß die statistischen 
Funktionen der Fehlerkurve und der Fehlerintegralkurve die gemessenen Impulsverteilungen 
wesentlich besser wiedergeben als die mit dem Mischungswegansatz (1) erhaltenen Funk- 
tionen. Außerdem ist die rechnerische Behandlung unvergleichlich einfacher als die mathe- 
matische Durchführung der Mischungswegtheorie 


4. Entwicklung einer phänomenologischen Theorie der freien Turbulenz 
auf Grund der molekularen Analogiebetrachtung. 


Im folgenden wollen wir nun unsere Theorie kurz wiedergeben. 


Die Impulsverteilungen in ebenen Freistrahlen (Schlitzstrahlen) und in ebenen Nachlauf- 
gebieten sind nach den vorliegenden Messungen in guter Näherung durch die Gleichung 


C, 


darstellbar (hier ist die Länge b ein mit «= veränderliches Maß für die Breite der Ver- 
mischungszone und e, und ce, sind Konstanten). 

Diese Gleichung ist eine Lösung der Differentialgleichung (3), wenn für die „Über- 
tragungsgröße*“ gesetzt wird, 

Da die Impulsverteilung an der Grenze eines ebenen Strahles, der an ruhende Flüssig- 
keit grenzt. näherungsweise durch eine Fehlerintegralkurve darstellbar ist, so dürfen wir auch 
ganz allgemein für den Fall der ebenen Berührung zweier Parallelstrahlen ansetzen: 


Auch diese Gleichung ist eine Lösung der Gl. (3), wenn wir für die Übertragungsgröße 
wiederum die Beziehung (5) einführen. 

Die Impulsverteilung in größerer Entfernung hinter einem ebenen Gitter können wir 
in guter Übereinstimmung mit vorliegenden Messungen '") durch die Gleichung wiedergeben: 


(Hier ist 5 die Gitterteilung. Der Exponent f(x) scheint nach den vorliegenden Messungen 
proportional x zu sein.) 
Die Gl. (7) ist ebenfalls eine Lösung der Gl. (3), wenn für .1 gesetzt wird: 


Damit haben wir die Übertragungsgröße für die beiden Fälle der veränderlichen und 
der konstanten Vermischungsbreite bestimmt. Allerdings wird durch diese rein formale Be- 
stimmung von 4 keine tiefere physikalische Aussage erzielt. (Die physikalische Bedeutung 
von A wird man wohl erst dann klar erkennen, wenn es gelungen sein wird, die Vorgänge 
der freien Turbulenz durch eine Statistik der Wirbelballen darzustellen.) 

Nachdem wir nun festgestellt haben, daß die Gl. (3) die Gesetzmäßigkeit der Impuls- 
übertragung für das uns bekannte Gebiet der freien Turbulenz in sehr guter Näherung wieder- 
gibt, haben wir nach dem „Ansatz“ zu fragen, der erforderlich ist, wenn aus der Bewegungs- 
gleichung die Differentialgleichung (3) folgen soll. 


schied zwischen dem gemessenen Gesamtdruck p + : u? und dem Impuls 2 #2 von der Größenordnung der Abweichun- 


gen der Meßpunkte von der Fehlerkurve bzw. der Fehlerintegralkurve. 

Es wäre übrigens auch denkbar, daß nieht 0%2, sondern p+0«2? statistisch verteilt wird (die Gl. (10) wäre dann 
dementsprechend abzuändern). Solche Feinheiten, die an sieh von grundsätzlicher Bedeutung sind, können aber zur 
Zeit noch nieht beurteilt werden, da ein hinreichend genaues Verfahren zur Bestimmung des statischen Druckes noch 
fehlt. 

%) Eine ausführliche Darstellung der Gesetzmäßigkeiten der freien Turbuleuz auf Grund der an dieser Stelle 
entwickelten Theorie erscheint demnächst in einem VDI-Forsehungsheft. Dort wird auch die Brauchbarkeit der stati- 
stischen Funktionen für die Darstellung des Impulsverlaufes an Hand zahlreicher Messungen im einzelnen nachgewiesen. 

10) R.Gran Olsson: Gesehwindigkeits- und Temperaturverteilung hinter einem Gitter bei turbulenter Strömung. 
Z. angew. Math. Mech. Bd. 16 (1936), S. 257. 
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Die Bewegungsgleichung des ebenen Problems lautet in. der Impulsform 


/(P dur _ 
wenn Zähigkeitskräfte zu vernachlässigen sind. Wenn wir auch noch die Druckunterschiede 
in erster Näherung vernachlässigen (was bei den Vorgängen, die zur Gl. (3) führten, sicher 
erlaubt ist), so ersehen wir, daß die Gl. (3) nur dann aus der Impulsgleichung erhalten wird, 
wenn für «® die Beziehung gilt 

du? 


Diese Beziehung über die Abwanderung des w-Impulses in die »-Richtung ist also das 
Ergebnis unserer induktiv entwickelten Theorie der freien turbulenten Strömungen. 

Das durch Gl. (10) dargestellte „Impulsübertragungsgesetz“ besagt: Der Strom von 
u-Impuls,derindie Querrichtung yabwandert, ist proportional dem «-Impuls- 
gefälle in dieser Richtung. 


Die Impulsübertragung erfolgt in Ebenen, die durch die «-Richtung und den Gradienten des 
u-Impulses beschrieben werden. Das Übertragungsgesetz lautet daher in Vektorschreibweise 


wenn i bzw. j die Einheitsvektoren der Längs- bzw. der Querrichtung bedeuten (die EEE wird 
durch (10a) definiert). 

Die Darstellung (10a) ist vom Bezugssystem unabhängig. Denn i ist kein willkürlicher Vektor, 
sondern © beschreibt die tatsächliche Bewegungsrichtung der Flüssigkeit vor ihrem Eintritt in den Be- 
reich der turbulenten Vermischung. 

Das Ergebnis unserer Untersuchungen ist nun insofern recht merkwürdig, als wir wider 
Erwarten keinen Ausdruck für die turbulente Scheinreibung gefunden haben. Denn die Gl. (10) 
gibt lediglich an, wie der Gesamtimpuls out=ow"-+our ausgebreitet wird. Über das 
Verhalten von ou’ v’ allein wird nichts ausgesagt, bzw. der Mechanismus der freien Turbulenz 
verfügt lediglich über die Summe des stationären und des instationären Anteil des Quer- 
transportes, ohne daß ein Unterscheidungsmerkmal zwischen diesen Gliedern erkennbar wäre’). 

Natürlich kann man Gl. (10) nach w?’ auflösen. Es ist dann aber noch die Frage, ob 
eine solche Gleichung für «r’ ein für die freie Turbulenz geltender Sonderfall einer univer- 
sellen Schubspannungsformel darstellt. 


Eine mögliche Form für einen allgemeinen Schubspannungsausdruck wäre beispielsweise 


Der Faktor x ist bei freier Turbulenz 1 (wodurch (11) mit (10) identisch wird) und bei Nicht- 
turbulenz oder bei Turbulenz ohne Se 

Bekanntlich ist das Auftreten einer turbulenten Schubspannung an das Bestehen einer Korrelation 
zwischen den Längs- und den Querbewegungen geknüpft. Der Faktor k dieser Korrelation verschwindet 
nun nicht nur für gewisse turbulente Strömungen (beispielsweise für die abklingende, isotrope Turbu- 
lenz), sondern auch für periodische Wirbelbewegungen (Karman-Wirbel)'?). Da sich die freie Turbu- 
lenz aus groben periodischen Störungen entwickelt, ist der Versuch naheliegend, den Faktor x als 
Funktion des Korrelationskoeffizienten k anzusetzen beispielsweise durch x<—=Zahl A). 


5. Die Ausbreitung der Vermischungszonen in der x-Richtung. 


Bei den vorangegangenen Betrachtungen ist der Verlauf der Breite der Vermischungs- 
zone in ihrer Abhängigkeit von der Koordinate « noch offengelassen worden, dieser ist 
daher noch zu klären. In den bisherigen mit dem Begriff des Mischungsweges operierenden 
Betrachtungen wurden von Anfang an nur solche Fälle in Betracht gezogen, wo die Ge- 
schwindigkeitsprofile in den verschiedenen Entfernungen x affin zueinander waren. Es er- 
geben sich hier durch Einsetzen des Ausdrucks für die scheinbare Schubspannung in die 
hydrodynamischen Differentialgleichungen jeweils automatisch die Potenzen von x, nach 
denen die Breite des Vermischungsgebietes zunimmt und die Geschwindigkeit abnimmt (soweit 
diese nicht wie beim Strahlrand konstant bleibt). 


1) em statistischen Standpunkt aus sind die Anteile #7’ und x offenbar gleichberechtigt bzw. die Kenn- 
zeichen dieser Glieder sind lediglich hydrodynamischer Natur. Die Gleichartigkeit der Anteile % # und 77 erscheint 
insofern plausibel, als ja.auch die stationäre v-Bewegung eine unmittelbare Folge der Einwirkung der Turbulenz auf 
die Flüssigkeit ist. (Ohne die Turbulenz wäre bei der vorausgesetzten Reibungslosigkeit keine ®-Komponente vor- 
handen und es würde sich nur die u-Gesehwindigkeit an den Strahlrändern diskontinuierlich verändern.) 

12) W.Tollmien: Über die Korrelation der Geschwindigkeitskomponenten von periodisch schwankenden Wirbel- 
verteilungen. Z. angew. Math. Mech. Bd. 15 (1935), S. 96. 
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Da in der neuen Theorie über das Impulsprofil von vornherein schon verfügt ist und 
zwar so, daß auch hier die verschiedenen Impulsprofile ein und derselben Strömung zuein- 
ander affin sind, genügt es, diese Affinität formelmäßig zu erfassen, wie es in den Gl. (4) 
und (6) geschehen ist, und zuzusehen, durch welche Potenzen man die Gl. (3), (5) und (10) 
befriedigen kann. Hierzu ist nur noch die Annahme nötig, daß die Strömungen einschließ- 
lich der Schwankungen in den einzelnen Querschnitten zueinander ähnlich verlaufen, daß 
also beispielsweise die Schwankungsgeschwindigkeiten in entsprechenden Querabständen be- 
stimmte Bruchteile des maximalen Geschwindigkeitsunterschiedes «, „.„, sind und daß auch 
“ev proportional zu u» gesetzt werden darf. Dann ergibt sich aus Gl. (10) unter Ver- 
wendung der Gl. (5) für .1, daß 


nur von y/b abhängt. 
Für die Geschwindigkeit ist im allgemeinen 
u=U+u, 


zu setzen, wobei [’ eine konstante Geschwindigkeit bedeutet und die Geschwindigkeits- 
unterschiede gegen durch ausgedrückt sind. 


Beim ebenen Nachlauf kann man für kleine ir näherungsweise «= — 2 U u, schreiben. 
Es folgt daher aus (12) |da aus Ähnlichkeitsgründen - "nur von rn abhängt ‚daß 
1max 
U db 
max dr 


ist. Daraus folgt weiter (da nach dem Impulssatz bier « b = konst. ist) 


1 11aX 


Ein solcher Gang der Vermischungsbreite ist für den ebenen Windschatten seit den Messungen 
von Schliehting bekannt 


Beim ebenen Freistrahl in ruhender Umgebung sind « und «, identisch. Bei der ebenen 
Berührung zweier Strahlen ist «, „., eine Konstante. In beiden Fällen ist aus Ähnlichkeits- 


db 
gründen u?/w} „a, nur von y/b abhängig. Es folgt somit aus (12), daß 2 — konst. bzw. 


ist. Damit haben wir die für Strahlen bekannte lineare Ausbreitung der Vermischungsgebiete 
richtig wiedergegeben. 


6. Zusammenfassung. 


1. Es wurde festgestellt, daß die Impulsverteilungen in turbulenten Ausbreitungs- 
bewegungen durch die statistischen Funktionen der Fehlerkurve und der Fehlerintegral- 
kurve in guter Näherung beschrieben werden. Die turbulenten und die molekularen 
Ausgleichsvorgänge können daher in gewisser Weise miteinander verglichen werden, so dal; 
die smngemäße Anwendung der Differentialgleichung der Wärmeleitung auf die Impulsaus- 
breitung in freien turbulenten Strömungen möglich ist. 


2. Aus dieser Differentialgleichung des Impulses und der hydrodynamischen Bewegungs- 
gleichung folgt, das der gesamte, in die Querriehtung abwandernde Strom des w-Impulses 
(der teils durch die Schwankungsbewegungen, teils durch die stationäre (uerbewegung her- 
vorgerufen wird) proportional ist dem « Impulsgefälle in dieser Richtung. Als Proportionalitäts- 
faktor erscheint die „Übertragungsgröße* 1 des Impulses, die als Funktion der Vermischungs- 
breite durch die Theorie selbst bestimint ist. 

3. Mit Hilfe der Formel für die Übertragungsgröße können die bekannten Ausbreitungs- 
gesetzmäßigkeiten der Vermischungsgebiete auf Grund von Ähnlichkeitsbetrachtungen richtig 
wiedergegeben werden. 313 


13) H, Schlichting: Über das ebene Windschattenproblem. Ing.-Arch. Bd. 1 (1950), S. >33. 
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Turbulente Ausbreitung eines Stromes erhitzter Luft. 
Von Wilhelm Schmidt in Bochum. 


(Aus dem Kaiser Wilhelm-Institut für Strömungsforschung.) 
I. Teil®) 


Der Prandtlsche Mischungswegansatz, dessen Verwendung bei der Berechnung 
einer Reihe turbulenter Ausbreitungsvorgänge yute Übereinstimmung mit 
den Experimenten ergab, wird auf die Berechnung der Ausbreitung eines 
ebenen und eines räumlichen Strahls erhitzter Luft angewandt. Die Aus- 
hreitung des ebenen Strahls wurde auch nach dem Ansatz von Taylor be- 
rechnet. Über die Ausbreitung des runden Strahls und seine e.rperimentelle 
Untersuchung wird in dem im nächsten Heft erscheinenden II. Teil der 
Arbeit berichtet. 


Einleitung. 

In turbulenten Strömungen treten neben den durch die unregelmäßigen Molekular- 
bewegungen hervorgerufenen Schubspannungen die durch die ungeordneten Mischbewegungen 
molarer Massen verursachten „scheinbaren Schubspannungen“ auf. Beide genügen dem New- 
tonschen Reibungsansatz : 

Schubspannung proportional dem Geschwindigkeitsgefälle. Im besonderen schreibt man 
für die von der Zähigkeit herrührenden Schubspannungen 


Try > u USW, 
xy Yy 


und für die scheinbaren Schubspannungen 


xy 


USW., 

wobei A eine zunächst unbekannte Austauschgröße ist. In den meisten Fällen kann r,, als 
klein gegen r/,, vernachlässigt werden. Um für die Austauschgröße A und damit auch für 
die scheinbare Schubspannung r/,„ einen Ansatz zu erhalten, denkt man sich die augenblick- 
liche Geschwindigkeit zerlegt in einen zeitlichen Mittelwert und die darüber gelagerten 
Schwankungen: 

u=ü+uw, 

(u, Mittelwerte; «', w’ = Schwankungen). 


Die Schwankungen bewirken einen Impulstransport, als dessen Rückwirkung die scheinbaren 
Schubspannungen anzusehen sind. Das Gesamtsystem der Spannungen bildet einen Spannungs- 
tensor, dessen Komponenten sich in folgender Form darstellen lassen: 


Ty=-owWv usw. 


Da die Schwankungen «’, v’, w’ aber unbekannt sind, lassen sich die Bewegungsgleichungen 
mit diesen Ansätzen nicht integrieren. Man muß vielmehr die Spannungskomponenten durch‘ 
die zeitlichen Mittelwerte der Strömung ausdrücken. Prof. Prandtl [4]') gelang es, solche 
Ansätze zu machen, indem er den Begriff des Mischungsweges einführte. 

Um zu diesen Ansätzen zu gelangen, denken wir uns eine ebene Strömung in Richtung 
der x-Achse mit einem Geschwindigkeitsgefälle senkrecht zur Strömungsrichtung. Den in der 
Strömung stattfindenden Vermischungsvorgang erklärt man dann durch die Annahme, daß 
innerhalb der Hauptbewegung einzelne Flüssigkeitsballen eine Eigenexistenz führen. Diese 
Eigenexistenz dauert nur für eine bestimmte Wegstrecke an, dann vermischen sich die 
Flüssigkeitsballen mit der neuen Umgebung. Wir betrachten nun einen einzelnen Flüssigkeits- 
ballen, der sich mit der Eigengeschwindigkeit »’ senkrecht zur Strömungsrichtung bewegt, 
und der von seiner Ausgangsstelle die «-Komponente der dort herrschenden mittleren Ge- 
schwindigkeit mitbringt. Diese Geschwindigkeit kann er nur bis zu seiner Vermischung mit 
einer neuen Umgebung, also nur längs einer gewissen Wegstrecke I, beibehalten. Die Weg- 
strecke 1, längs der der Ballen die Geschwindigkeit des Ausgangsortes ungeändert beibehäit, 
nennt man nach Prandtl den „Mischungsweg“. Der Geschwindigkeitsunterschied, den der 
Ballen durch seine Verschiebung gegenüber der neuen Umgebung gewinnt, stellt sich mit 


Hilfe des Mischungsweges dar durch 4 und in erster Näherung kann man daher die 


*, Der 2. Teil wird im nächsten Heft erscheinen. 
1) Die Zahlen in eckigen Klammern verweisen auf das Sehrifttumverzeichnis am Schluß der Arbeit, 
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mittlere Geschwindigkeitsschwankung «’ diesem Ausdruck proportional setzen 


Die Quergeschwindigkeit »’ kann man sich dadurch entstanden denken, daß zwei Flüssigkeits- 
ballen mit verschiedenem «’, die sich voreinander befinden, zusammenprallen oder sich von- 
einander entfernen. Die entstehenden Geschwindigkeiten r’ können also proportional «’ ge- 
setzt werden. Für die Schubspannung erhält man demnach, wenn die Proportionalitätsfaktoren 
in « und v’ und der Korrelationsfaktor in w’v’ in das noch unbekannte I hineingenommen 


werden, 


und wenn man noch berücksichtigt, daß r/;„ mit Er sein Vorzeichen ändert, 


Der Vorteil dieses Ansatzes besteht darin, daß man in manchen Fällen Aussagen über das 
Verhalten von ! machen kann, während dies für A nicht möglich ist. Insbesondere nimmt 
Prandtl in Fällen der freien Turbulenz an, daß 
1. für genügend große Reynoldssche Zahlen in Querschnitten senkrecht zur Strö- 
mungsrichtung die Vorgänge geometrisch und mechanisch ähnlich verlaufen, und daß 
2. ! über den ganzen Querschnitt jeweils konstant ist. 
Aus der zweiten Annahme folgt dann, daß ! der Breite b des Strahles proportional ist. 
Einen ähnlichen Ansatz für die scheinbare Schubspannung gewinnt G. I. Taylor [9] 
durch eine vom Prandtlschen Gedankengang abweichende Überlegung. 
Taylor setzt voraus, daß in der betrachteten turbulenten Strömung nicht nur die 
Hauptbewegung, sondern auch die turbulente Mischbewegung als zweidimensional angesehen 
werden kann. Führt man nun die Rotation 


1 =) 


0x 
- sie ist nach Taylor die Größe, die dem Austausch unterliegt - in die Eulersche Be- 
wegungsgleichung ein, so läßt sich diese 
Zur Vereinfachung der Theorie nimmt Taylor eine Strömung mit den mittleren Geschwin- 
digkeiten 
U=fy und V=0 an. 
Bildet man nun die zeitlichen Mittelwerte und drückt die Schwankung mit Hilfe des Mischungs- 
weges L aus, so lautet die Gleichung Taylors 
Der Mischungsweg 7, ist nicht mit dem der Prandtlschen Theorie identisch, sondern ist 


eine für jedes einzelne individuelle Teilchen verschiedene positive oder negative Länge. Mit 
einer anderen ebenfalls veränderlichen Länge 7’ wird 


v=L' 
dy 
gesetzt. Damit erhält man 


Diese Gleichung läßt sich, wenn I von y unabhängig ist, integrieren, und man erhält für die 
scheinbare Schubspannung in 


1 „faUV 


einen Ansatz, der sich von dem Prandtlschen durch den Faktor 1/2 unterscheidet. 


2) — bedeutet proportional. 
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In ähnlicher Weise wie der Impuls werden auch andere Eigenschaften der Flüssig- 
keiten, wie Temperatur, Farbe, Staubgehalt usw. ausgetauscht. Beim Temperaturaustausch 
bekommt man nach Aufspaltung der augenblicklichen Temperatur in einen zeitlichen Mittel- 
wert 9 und die Schwankungen #’ für die turbulente Wärmeleitung den Ansatz 


| 


wobei ec, die spezifische Wärme bei konstantem Druck ist. In Analogie zu der Ableitung 
des Ausdrucks für r,, erhalten wir hier 

Damit wird der Wärmetransport durch die Flächeneinheit in der Zeiteinheit z durch turbu- 
lente Wärmeleitung mit AV =1? 
du 
Erikr (3). 
Auch hier werden die Proportionalitätskonstanten und der Korrelationsfaktor in « ©’ wieder 
in ! hineingenommen. Neben der turbulenten Wärmeleitung tritt immer die durch Molekular- 


bewegung hervorgerufene Wärmeleitung 


auf, die aber als klein gegen x vernachlässigt werden kann. 

Während man allgemein annimmt, daß die Austauschgrößen in der scheinbaren Schub- 
spannung und der turbulenten Wärmeübertragung vollkommen gleich sind, daß also auch die 
Mischungswege identisch sind, ergibt sich aus der Taylorschen Theorie ein Unterschied. 
Der Mischungsweg des Temperaturaustauschs ist nach Taylor y2mal so groß wie der des 
Impulsaustauschs. 

Mit Hilfe des Prandtlschen Schubspannungsansatzes wurden Beispiele von Toll- 
mien [8], Schlichting [6] und L. M. Swain [7] behandelt, wobei sich verhältnismäßig 
gute Übereinstimmung mit den Experimenten ergab. 

Aufgabe dieser Arbeit ist es, die Geschwindigkeits- und Temperatur- 
verteilung in einer durch Übertemperatur gegenüber der ruhenden Luft 
hervorgerufenen Strömung im ebenen und rotationssymmetrischen Fall 
nach den Ansätzen von Prandtl und der Annahme von der Identität der 
Mischungswege und im ebenen Fall auch nach G. I. Taylor zu ermitteln. 


Das ebene Problem. 


Längs einer waagerechten geraden Linie, die so lang gewählt wird, daß ein ebenes 
Problem vorliegt, wird der Luft Wärme zugeführt. Der entstehende Luftstrom vermischt sich 
an seinem Rande mit der umgebenden ruhenden Luft. Die auftreibende Kraft wird hervor- 
gerufen durch die durch die Temperaturdifferenz gegen die ruhende Luft bedingten Dichte- 
unterschiede. Die x-Achse legen wir in Richtung der Strömung und die y-Achse senkrecht 
zur Strömungsrichtung und senkrecht zur erwärmten Linie. 

Bezeichnen wir mit 7 die absolute Temperatur des Strahles, mit 7, die der ruhenden 
Luft und mit o und o, die entsprechenden Dichten, dann erhalten wir für die archimedische 
Auftriebskraft 


Die in dieser Beziehung auftretende Dichtedifferenz läßt sich mit Hilfe des Gay-Lussac- 
schen Gesetzes umformen in 

Dabei sind $ der Ausdehnungskoeffizient und # die als klein angenommene Temperatur- 
differenz T— T,. Die Auftriebskraft wird damit 


Mit dieser Größe für die Auftriebskraft lautet die Bewegungsgleichung 
gu du 
00x 


gegen 


wenn die Zähigkeitsschubspannung r,, und die Änderung der Normalspannung 


dx 


dy 
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T,„ vernachlässigt werden. Diese Gleichung enthält die drei unbekannten Funktionen u, v 
und #9, zu deren Ermittlung noch zwei weitere Gleichungen, nämlich die Kontinuitätsgleichung 


ou 


und die Wärmetransportgleichung 


zur Verfügung stehen. In der letzten Gleichung müßte die Wärmeleitung 


4 
ebenfalls vorkommen, kann aber als klein gegen x vernachlässigt werden. 

Bei der Untersuchung der Strömung in ihren Einzelheiten machen wir von den S. 266 
angegebenen Annahmen über den Mischungsweg Gebrauch. Wenn wir mechanische und 
geometrische Ähnlichkeit der Vorgänge in den Querschnitten senkrecht zur Strömungsrichtung 
annehmen, folgt für den Mischungsweg, daß er der Strahlbreite 5 proportional sein muß 


I=eb (e=konst.). 


Den Zusammenhang zwischen der Strahlbreite und der Entfernung vom Spalt ermitteln wir 


durch die folgende Abschätzung : Es ist v’ x Für einen mittleren Wert von v’ ergibt 
sich Proportionalität zwischen »’ und 

‚la 

Mit I=eb folgt daraus 


Außerdem folgt aus der Aussage 


(4, substantieller Differentialquotient) 


also 
b=r-.x, wobei r eine reine Zahl ist. 


Für den Mischungsweg ! können wir also schreiben: 


Wenn wir zu dieser Gleichung den Impulssatz und den Wärmetransport der Hauptbewegung 
hinzunehmen, bekommen wir den Zusammenhang der Geschwindigkeit und der Temperatur 
mit dem Abstand von der erwärmten Linie. 

Wir setzen voraus, daß der Wärmetransport durch die Querschnitte senkrecht zur 
Strömungsrichtung konstant ist, daß also 


und damit 
onst. 


ist. Aus dem Impulssatz, der in unserem Falle aussagt, daß die Änderung des Impulses 
längs des Weges dx gleich der Auftriebskraft längs dieser Strecke ist, gewinnen wir die 


Gleichung 


Mit 9=Dx " (m zunächst noch unbekannt: D ist eine Konstante, die eingeführt wird, da- 
mit die Gleichung dimensionsrichtig wird) 


Konst. und 5b=r-x wird der Impuls 
r x - Konst.? 
J= 


db 

u= odx . . . . . . . . . . . . . . . . 
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Nehmen wir im folgenden 59 als so klein gegen 1 an, daß es gegen I vernachlässigt werden 
kann, dann kann an Stelle von o überall o, geschrieben werden. Es wird somit 


Konst.’r x 


0 D? +? 
o 
und 


Für den Impuls gilt dann die Beziehung 


und für die Impulsänderung erhält man 


” 
Gl. (9a) geht über in 
9 
. (9b). 
Ersetzt man auch in der Auftriebskraft o durch o, gemäß 
P9) 
und begeht dieselbe Vernachlässigung wie beim Impuls, so wird diese 


Setzen wir nun (10) und (11) einander proportional, dann erhalten wir 


- Mm 
x 


konet., 


also 
w—i. 


Die Temperatur # ist also der ersten Potenz des Abstandes von der Wärmequelle umgekehrt 


proportional. Wird nun 9= un in (9b) eingesetzt, so zeigt sich, daß die Geschwindigkeit 


von dem Abstand von der erwärmten Linie unabhängig ist. Mit der Dimensionslosen 7 = F 
ergeben sich demnach für « und d die Ansätze: 


1 
und 


Die neben « auftretende Geschwindigkeitskomponente » ist mit dieser durch die Kontinuitäts- 


gleichung verknüpft. Man genügt der Kontinuitätsgleichung durch Einführung einer Strom- 
funktion y. Ausgehend von 

erhält man durch Integration 
v=\fo)dy 
Bezeichnen wir \f(y)dn mit F (n), dann wird 


und mit 


y=zF(n). 


Durch Differentiation nach y und « gewinnt man daraus die Geschwindigkeitskomponenten 
u und v. 


Konst. 
o—. 
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wobei ’ Differentiation nach », bedeutet. 
Um die einzelnen in den Gleichungen vorkommenden Glieder in dimensionsloser Form 
zu erhalten, setzen wir 
u=U®d; v=U(-®+n®) 
und 


U ist die von « unabhängige Geschwindigkeit in der Mitte des Strahls. 


IE... 


Lösung des Problems mit dem Prandtlschen Schubspannungsansatz. 


und 
PP" U? 


erhält die Bewegungsgleichung die Form 

Man erkennt, daß r/,, für positive 7 das negative Vorzeichen erhalten muß, wenn man be- 
rücksichtigt, daß die zu erwartende Geschwindigkeitsverteilung in der Mitte ein Maximum 


haben und nach dem Rande zu auf Null absinken muß. Da die Rechnung nur für die rechte 
Hälfte durchgeführt wird — für die linke Hälfte erhält man das Ergebnis durch Spiegelung —, 


muß T und damit auch r/, negativ sein. Das Entsprechende gilt auch für den Wärme- 


transport #. Multiplizieren wir die Wärmetransportgleichung mit %g, dann erhalten wir mit 
den Ansätzen (12a) und (3) und den Differentialquotienten 


x? 
die Wärmetransportgleichung 
Diese Gleichung läßt sich sofort einmal integrieren und lautet nach Umordnung der Glieder 


Um den in (15a) und (16a) auftretenden, unbekannten Faktor ec? zu beseitigen, machen wir 
die Substitution 


> 


mit noch unbekanntem a. Dann wird, wenn jetzt die Striche Differentiationen nach 1, bedeuten, 
d® 


Eingesetzt in (15a) und (16a), erhält man damit 


und 
=0. 
Setzen wir nun die Zahlenfaktoren der linken Seiten von (15b) oder (16b) gleich, dann folgt 


1 


1 


so daß die beiden Differentialgleichungen die Form annehmen: 
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Multiplizieren wir nun (16) mit (2 eayt3 und bezeichnen das Produkt (2 ec?) sg, mit ©,, dann 
lauten unsere Differentialgleichungen 

und 


Zu diesen Differentialgleichungen treten die Randbedingungen hinzu. Diese sind 


1. für die Mitte des Strahls 
9,=1 durch Normierung erreichbar, 
“„=endlich d.h. ®'= endlich, 


d.h. P=0; 
2. für den Strahlrand 


0, —0 wegen der Stetigkeit des Anschlusses nach außen, 
P"’—0 ebenfalls wegen der Stetigkeit des Anschlusses. 


Die Lösung unseres Differentialgleichungssystems in geschlossener Form ist nicht mög- 
lich. Es müssen für ® und ©, Reihenansätze gemacht werden. Dabei muß aber eine Eigen- 
tümlichkeit des Prandtlschen Schubspannungsansatzes beachtet werden. Die Lösungen mit 
Hilfe dieses Ansatzes weisen nämlich gewisse Singularitäten auf. Man bekommt keinen 
asymptotischen Übergang zu konstanter Geschwindigkeit, sondern einen Anschluß nach außen 
mit endlicher Krümmung des Geschwindigkeitsprofils, und beim Geschwindigkeitsmaximum 
erhält man sogar unendliche Krümmung. Die Funktion ®’ verhält sich dort, d. h. bei e=0 


wie y=4+ Ba:. Entsprechendes gilt auch für die Temperaturverteilung. Die Singularitäten 
treten dort auf, wo wegen aym der Austausch verschwindet. In Wirklichkeit verschwindet 
dort der Austausch aber nicht vollständig. Der Vergleich der mit diesem Ansatz gewonnenen 
Ergebnisse mit den experimentellen Ergebnissen zeigt aber hinreichend gute Übereinstimmung, 


so daß die Verwendung dieses Ansatzes gerechtfertigt ist. Um die Singularität ganz zu ver- 
meiden, hat Prof. Prandtl! einen verbesserten Ansatz angegeben, der die Form 


’ „ou u\? 


hat. Dieser Ansatz soll jedoch mit Rücksicht auf die Durehführbarkeit in unserer Rechnung 
nicht verwandt werden. 


Dem Auftreten der Singularität trägt man in den Lösungsansätzen dadurch Rechnung, 
3 
daß man Reihen ansetzt, die nach Potenzen von n,? fortschreiten. 


5 ı1 


Wegen 9, 0)=1 ist ,—1. 
Durch Einsetzen der Ansätze (20) und (21) in die Differentialgleichungen und durelı 


Vergleich gleich hoher Potenzen von n, errechnen sich die Koeffizienten, die zum Teil noch 
von a, abhängen, zu 


6 
| 
a, = 3402 do; 4,=— 33856525 V a, d,= 0,000000435 do; 
2 8 /2 


2 608 431 
V 0,0001023 a} + 0,0001085 usw. 


Eine Angabe über die Größen a, und den Randpunkt 7,x erhalten wir aus den Rand- 


bedingungen 
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Da die Konvergenz der Reihen, besonders die der Reihe für ©, in der Nähe des Randes 
schlecht wird, setzen wir eine weitere Reihenentwicklung um den Randpunkt n,z an, die 
gerade in dessen Nähe gut konvergiert. Wir setzen 


(22) 

und erhalten die Lösungsansätze: 


Wegen der Randbedingungen ist 
a, =4,=b,=b, 


Geht man mit diesen Ansätzen in die Differentialgleichungen, so liefert der Vergleich gleich 
hoher Potenzen von 7, die Koeffizienten 


Die Größe von @, läßt sieh angenähert angeben. Es ist nämlich 
P 


Da aus P’(n,R)= P’(n,R)=O ıR angenähert bekannt ist, man also auch einen Näherungs- 
wert für 7, besitzt, ist b, die einzige Unbekannte. Über b, erhält man eine Angabe, wenn 
man die Werte beider Entwicklungen für ®, wie sie sich mit den Näherungswerten 7, x, ds, 
ä,, 7, ergeben, in einem Anschlußpunkt übereinstimmen läßt. 


Verbesserung der Näherungswerte. Zur Verbesserung der zunächst nur an- 
genäherten Werte wenden wir das Newtonsche Verfahren in einer Form für 4 Variable an. 
Wenn wir verlangen, daß sich die Entwicklung um den Nullpunkt und den Randpunkt mit 
stetiger Tangente und stetiger Krümmung aneinander anschließen sollen, dann bekommen wir 
8 Entwicklungen, deren Differenz zu je zweien im Anschlußpunkt Null sein muß. Wegen der 
nur angenäherten Bestimmung von «a,, @,, b, und ») ist diese Differenz aber nicht 0, sondern 
nimmt die verschiedenen Werte &; an. Es ist also 


wobei die &; Funktionen der vier Unbekannten sind 
Die ersten Näherungswerte der Unbekannten sind 
5=04; V=06. 
Die verbesserten Werte «a,. @,, b,, 7 stellen sich mit Hilfe der ersten Näherungen und der 
Verbesserungen h, k, I, m dar in der Form 


Entwickeln wir nun die &; nach Potenzen von h, k, I, m und vernachlässigen die Glieder von 
höherer als der ersten Ordnung, so können wir die Verbesserungen aus dem System der 
folgenden vier Gleichungen bestimmen: 


Die an die &; angehängten Indizes «,, @,. b,, 7 bedeuten Differentiation nach diesen Größen. 
Für die Verbesserungen bestimmen sich aus dem Gleichungssystem (25) die Werte 
h—= 0,0002; k=0,00178;  1=0,00051;  m= — 0,0035. 
Die verbesserten Werte sind dann «a, 
a,=1139%8: 4,—=0,%822; 02451; 7=0,5%5. 


3) Die angehängten Indizes geben an, ob es sich um die Entwicklung um den Nullpunkt oder den RKandpunkt 
handelt. 
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Die Reihen für ® und ©, um den Nullpunkt sind damit 
5 1 
D = 1,13982 7, — 0,377129,? +0,063322 — 0.005079 74° 
8 
11308 0,9287, +0,253288 — 0,0279 . (26a). 
8 
9, 1 — 1,0746035 n,? + 0,6108258 — 0,1154203 +... | 
Diejenigen um den Randpunkt (7,=n,r — n,) lauten 
D=0,%822  — 0,16124 7} + 0,010356 7} + 0,001596 73 - - | 
u 
—= 0,48372 73 0,041424 7} — 0,00798 7} . (6b). 
9, = 0,2404 73 — 0,061767 7; + 0,0193 74 + - | 
Die Werte der Funktionen ©,, ®, Ü und 77 | ” | | 
finden sich in Zahlentafel 1. Durch geeignete Z [6, | | 
Maßstabswahl für die Geschwindigkeiten kann 
man erreichen, daß die Geschwindigkeits- und 2 | | 
Temperaturprofile im Scheitel übereinstimmen. ı | | 
Die auf diese Art dargestellten Geschwindig- __ _ Geschwindigkeit -— —— Temperatur — 
keits- und Temperaturprofile finden sich in \ | 
| 
Bild 1. | | | 
Die einzige Unbekannte der Theorie ist | | | vll | | | 
die Konstante c, der Proportionalitätsfaktor r 1 Be 
zwischen dem Mischungsweg und der Ent- | | | 
fernung von der erwärmten Linie. Sie kann 
nur durch Vergleich mit dem Experiment er- | | 
mittelt werden. r a 04 N 
Lösung des Problems mit dem | er Bad 
Taylorschen Ansatz. a | | | 
Wenn wir für die scheinbare Schub- a 7 
spannung den Taylorschen Ansatz (2) Bild ı. 
Zahlentafel 1. 
1 0 1.1398 0 0 (0) 
0,1 0,9665 0,1128 1,1202 (0.9826) 0,11202 0.0007 0,00061) 
0,2 0,90666 0,2215 1.0582 (0,9284) 0,21164 — 0,00986 - 0,00865) 
0,3 0,8311 0,3154 0,9888 (0,8675) 0,29664 — 0,01879  (— 0,01649) 
0,4 0,75703 0,41995 0.91698 (0,8045) 0,3668 — 0,05256  (— 0,04611) 
0,5 0.067354 0,507 0,8369 (0,7343) 0,4185 — 0,089 (— 0,0781) 
0,6 0,59068 0,5866 0,7607 (0,6674) 0.4564 - 0,1302 (— 0,1142) 
0,7 0,51964 0,65775 0,66914 (0,587) 0,4684 0,1893 (— 0,1661) 
0,8 0.4324 0,7204 058575 (0,519) 0,4686 — 0,2518 0,2209) 
0,9 0,387 0,7749 0,5032 (0,4415) - 0,4529 - 0,322 (— 0,2825) 
1 0,294 0,8205 0,4246 (0,3725) 0,4246 0,3959  (— 0,3473) 
3 0,2365 0,8604 0.3525 (0,3093) 0,3878 — 0,47265  (— 0,4147) 
12 0,1818 0,8915 0,2812 (0,2467) 0.3375 - 0,554 (— 0,48605) 
1.3 0,1364 0,9166 0,2185 (0,1917) 0,2841 — 0,6324 (— 0,5548) 
1,4 0,0983 0,9356 0,1622 (0,1423) 0,22708 0,7085  (— 0,6216) 
1.5 0.0678 0,9495 0,1139 (0,0999) 0,16915 0,77865  (— 0,6831) 
1,6 0,0496 0.9584 0.0733 (0,0643) 0.11728 — 0,8411  (— 0,738) 
1,7 0.0233 0,9639 0,0415 (0,0364) 0,0705 — 0,8934  (— 0,7838) 
1,8 0,009 0,9667 0,0178 (0.01562) 0,032 — 0,9347 (0,8201) 
1,9 0,0025 0,9678 0.0048 (0.00421) 0,00912 0,9587 (— 0,8411) 
2 0 0,9682 0 (0) 0 0,9682  (— 0,8494) 
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verwenden, lautet unser Differentialgleichungssystem : 


Mit der Substitution 


im 


„=ay, a=c 
geht dieses über in $ 


Multiplizieren wir Gl. (28b) mit (e)t3 und setzen 


dann lauten die Taylorschen Gleichungen 
Dazu kommen wieder die Randbedingungen. Diese sind: 
1. für die Strahlmitte 

endlich d. h. ®= endlich, 

v 


d.h. $=0: 
2. für den Strahlrand 


.P"=0 (wegen der Stetigkeit des Anschlusses). 


Die Randbedingung 9; —=0 ist bei der Taylorschen Lösung des Problems nicht erfüllt, so 
daß die Temperatur mit einer Unstetigkeit in der Tangente nach außen anschließt. Im Null- 
punkt tritt auch hier die schon früher erwähnte Singularität, also unendliche Krümmung, auf. 
Die Lösung der Differentialgleichungen erfolgt wieder mit den Reihenansätzen (20) und (21) 


+bmM+...... 


Durch Einsetzen dieser Ansätze in die Differentialgleichungen und durch Vergleich gleich 
hoher Potenzen gewinnen wir für die Koeffizienten die Bestimmungsgleichungen: 


225 15 
34-744 = b, 
28a — 12 a,a, — = 
3465 3267 


63 
+ > %4,-7%4- = b; 


11475 13005 35343 255 183 147 


ma, + + —— 4,0, — 4 - - = 
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ab, 
180,1, 
— a,b, — a,b, — a,b, =0 
a,b, +54a,0,+ 05, a,b, a,b, — a,b, — a,b, 


225 1 76 
a,b, +126a,b,+ a,b, —a,b, — a,b, —a,b,—a,b, —a,b,=0 


135 297 765 
+ +7 0b,+135a,b, — a,b, —a,b, 


— 4,b, — a,b, — a,b, — a,b, =0 


und daraus die ersten Koeffizienten 


2 1 4 4 
13 137 68 


Die übrigen Koeffizienten wurden durch Einsetzen der bereits berechneten in die Bestimmungs- 
gleichungen gewonnen. Aus den Bedingungen 


Zahlentafel 2. 


| u u \9| v \» 
1 | 198 (0) 
012756 | 1,2509 (0,9716) | 0,1258 —0,00177 (— 0,001375) 
0.2 0.947 0,2508 | 1,2029 (0.9355) 024438 0,0064 0,00497) 
0,3 0902, 08678 (0,8823) 0,34087° 0,0269 (0,0209) 
| 04788 | 1,0682 (0.8256) 04258 --0,053 (0,0412) 
(07 | — 0,0897 (0,0697) 
06 | | 0,8964 (0,6962) | 0,5378 | —0,13965 (-- 0,1085) 
0,7589 0,809 (0,6283) | 0,5663 — 0,1985 (0,1508) 
| 072 (0,5608) | 0,5776 — 0,2567 (0,199) 
09 | 0,5658 005 0,6358 (04938) | 05722 |-—0,3328 (0,2585) 
1 05062 0.9686 | 0,5519 (0,1287) | 05519 — 0,4117 (0,3198) 
11 | | 1082 | 0497 (03671) | 051997  —-0,4952 (0,3846) 
12 | 088793 | 1,0582 | 0,3955 (0,3072) | 047456 | —-0,5836 (0,4533) 
13 | 0,3362 1.0941 0,3357 (02607) . 04864 0,6577  (-- 0,5108) 
14 | 08817 11234 | 0,25% (02016) 0,3633 0,7601 (0,5904) 
15) 080 | | 0,2012 (0,1563) | 0,3018 | —-0,8449 (0,6562) 
1.6 | 01858 | 11584 | 0,1467 (0,1139) | 02347 — 0,923 (0,7169) 
17 | 0145 1.1715 0,10298 (0,08) 0,9864 (0,7739) 
1,8 | 0,095 1.1799 (0052) 0,1207 1,0592 (0,8285) 
19 | 0,0689 1.1858 0,0405 (0.0315) | 0.07695 | —1,10903 (— 0,8614) 
1.188552 | 0,0212 (001647) | 004 | 1,146 (0,8901) 
2,09 0 ' 1,1894 (0) |— 1.1894 (— 0,9165) 


1) Die eingeklammerten Zahlen sind die auf die Scheitelhöhe 1 umgerechneten Geschwindigkeiten. 
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Diese Entwicklungen konvergieren in der Nähe des Randpunktes noch gut. Als Lösungen 
für ®, w und ©, haben wir 


7 
1,285 7, -- — 0,0042869n,* + 0,0002759 7} 


0,125583 73 — 0,014507 + 0,0011363 
Der Randpunkt 7,r liegt bei 2,09. 
u 


Die Werte der Funktionen 2, a und ©, finden sich in Zahlentafel 2, die dazu 


gehörigen Kurven in Bild 2. Die Profile sind wieder so gezeichnet, daß sie in den Scheiteln 
übereinstimmen. 


Der Vergleich des Taylorschen Ergebnisses mit dem nach Prandtl gewonnenen 
zeigt, daß bis auf Unterschiede in der Form und der Lage des Randpunktes das Temperatur- 
profil und das Geschwindigkeitsprofil ihre Rollen vertauscht haben. Während bei der 
Taylorschen Lösung das Geschwindigkeitsprofil schneller absinkt als das Temperaturprofil, 
ist bei der Prandtlschen Lösung gerade das Umgekehrte der Fall. Eine Entscheidung 
darüber, welche der beiden Theorien in unserem Fall richtig ist, könnte nur durch das 


Experiment erbracht werden. 

YA 720% EN \ 


08 Bild 9. 
Po 
Bild 2. 


Stromlinien. Die Gleichung der Stromlinien lautet: 


Läßt man Y nacheinander die Werte 7 r 3 .... 8... annehmen, so zeigen Zahlentafel 3 und 
das zugehörige Stromlinienbild, Bild 3, die Stromlinien nach Prandtl, Zahlentafel 4 und 
Bild 4 die Stromlinien nach Taylor. 


Berechnung der Druckunterschiede. Bisher war die zwar nicht genau gültige, 
aber für die erste Näherung ausreichende Voraussetzung gemacht worden, daß.der Druck in 
Querschnitten senkrecht zur Strömungsrichtung konstant ist. 


In zweiter Näherung sollen nun aber auch die durch die berechneten Geschwindigkeiten 
und Spannungen hervorgerufenen Druckunterschiede mit berücksichtigt werden. Wir gehen 
von der zweiten Bewegungsgleichung aus: 


dv dv (dr, 10 
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Zahlentafel3. 
vg vg v=3 
02 05565 0.118 1,18 0.226 1,6695 0,339 226 0,152 | 2,788 |0,5625 3,39 0,678 
0402982 0,1194 0.5964 | 0.2385 | 0.8946 | 0.358 | 1.1924 | 0,477 1,491 0,5962 1,789 0,715 
0,6 02131 0,128 |0,4262 | 02557 0.6393 | 03835 0,8524 | 05114 1,0655 0,6392 1.279 0.707 
08 0174 0,139 0,3475 10278 0.521 0417 0,695 | 0,556 | 0.8687 0,695 11.0425 0,834 
10,158 0,1577 031545 | 0,315 0,473 0,4132 0,631 | 0,631 0,7886 0,789 0.9464 0,9464 
12 0,1408 0,1688 0,2805 | 0.337 0421 0,5049 0,561 | 0,673 | 0,7012 0,842 0.8415 1.0098 
14 0,1323 0.1851 0.2645 | 0.3703 0,397 0,555 0,529 | 0,741 0,6612 0,926 0.7935 1,119 
1.6 0,1304 | 0,2082 0,26085 0.4174 0.3912 0.6245 0,522 0.837 0.6521 1.04085 0,7826 1.252 
18 0,1298 0.2335 0,2595 0,467 0,3893“ 0,7005 0,519 0.934 0.6486 1.1675 10.7785 1.381 
20.1291 0.258 0.2588 0,516 0.3874 0,714 0517 1.088 0,6456 1.29 0,7748 1,548 
02 3,8055 | 0,791 4.452 0904 904 1,818 
06 1,4917 0,895 11,705 1.0228 | 2131 12785 | 2,567 1.584 2,088 1.79 3,41 2,046 
08 1216 0978 11,19 1112 1.2375 139 2.085 1,668 2,433 1.946 2.776 2,2208 
11,108 | 1,104 51262 1.262 | 1,577 1,577 1,898 | 1,892 2,208 12,208 2,43 2,43 
12 09818 1.178 1,122 1,346 1,403 1,683 | 1.688 | 2.02 1.964 12,356 2.244 2,6928 
1.4 0,9258 | 1,296 1,058 | 1,481 | 1,33 | 1,852 1,587 | 2222 1.852 12,59 2.116 29624 
1.6. 0,91297 1,475 1,0834 | 1,669 | 1.304 2,087 1,565 2,504 1,826 |2,9215 2,087 |3,3392 
18 0,9083 1,635 1.038 | 1,868 | 1.208 2,325 1,557 | 2,802 1.817 3.269 2.076 | 3,7368 
209088 1,806 1.033 12,066 1.201 2,58 1,5495 | 3.099 1.808 3.612 2.066 14,132 
Zahlentafel 4. 
02 0,4984 0.09968 0.9968 0,1994 | 1,495 0,299 1,9936 0,399 2,192 | 0,498 12,99 0,598 
04 02614 0.1045 | 0,5227 0,209 0,784 0314 | 1,0454 | 0,118 1,307 0,523 11,568 10,627 
0,6 0,1844 0,11063| 0,3688 0,221 0,553 0,332 | 0,7375 | 0,4425 0,922 0,553 1,106 ‚0,664 
0,8. 10.1498 | 0,1109 0,2996 | 0.2397 0,4495 0,3506 | 0,599 | 0,479 | 0.749 0,599 10.899 0,719 
1 0,1287 0,1287 0,2595 | 0,2595 0,389 0,389 0,519 0.519 0.649 0,649 0,773 0,778 
12 0,1181 0,1417 0,2362 0,2834 0,354 0,25 0,472 0,567 0,591 0,7086 \0,709 0,85 
14 0,1113 0,1558 0.2226 0,3116 | 0.334 0,467 0,451 | 0,623 0,556 0,779 0,676 0,9347 
1,6. 0,10735| 0,1718 0.2147 | 0,3435 0,322 | 0,515 0,499 0,687 0,537 0,850 0,644 1,0306 
1,8 0,10544| 0,1898 | 0,2109 0,3796 | 0,316 0.500 0,422 10.759 0,527 0,949 0.633 1,1387 
0,10517| 0,2034 0,2103 0419 | 0.3155 | 0631 | 0,4207 | 0,8414 0,526 | 1,052 0,031 11,262 
2.09 0.2197 021026 0.430 03154 | 0.659 | 0.4205 | 08780 05257 1.099 10.629 |1.318 
7 1 1 
| y=1 | y=1- | y=?2 
r 
02 13,480 | 0,698 | 3,987 | 0,797 4984 0,0968 | 5,981 1190 | 6978 | 1,396 7,974 1,595 
0,4 1,8205 | 0,7318 2,0908 0,836 2,614 | 1,0454 13,136 | 1,254 3,659 1,464 4,182 1,673 
0,6 | 1,2907 | 0,774 1,475 | 0,885 1,844 | 1,106 12,218 | 1,328 | 2,581 1,549 12,95 1,77 
0,8 |1.049 | 0,889 1,1986 | 0,9589 1.498 | 1.1986 11,798 | 1,438 2008 1,678 12,397 11,918 
1 0,908 | 0,908 1,0378 1.0378. 1.297 1.297 11,557 1,557 | 1,816 | 1.816 2.076 |2,076 
1.2 |0,8267 | 0,992 0,9449 | 1,1339 1,181 1,417 1,417 1,701 11.658 | 1,984 11,89 2,27 
14 | 0,7789 | 10908 0.8902 | 1246 1.118 | 1558 11,335 | 1860 1,578 | 2181 1,7804 2,408 
1.6 107515 | 1,2028 08588 | 1,974 10785 1,718 | 1.288 | 2,061 11,503 | 2,405 1,718 2,748 
18 10,738 | 1,3285 0,8435 | 1,518 | 1.054 1.898 1,265 2,2775 1,476 | 2,657 1,687 13,037 
2 0,7362 1.4724 0,8414 | 1,683 | 1,0517 | 2,103 1,962 |2,5%4 | 1,4724 2,945 1,688 3,365 
2,09 0,7359 | 1,538 10,841 | 1,758 | 1,0513 | 2,197 | 1,2616 | 2,697 1.472 3,076 11,682 13,516 
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Diese Gleichung formen wir um, indem wir schreiben: 


dv 


Dann ist 


JaT3y ist aber im Bereich unserer Näherung wegen der Kontinuitätsbedingung gleich Null, 


so daß wir nach Integration nach y erhalten 


Setzen wir für r/,,, TT und o, ihre Werte ein — es ist 


du dv 
—, 


so lautet die Gleichung für die Druckverteilung 


-[2n — \ U? n.Prandtl . . (32a). 


2 ) J 
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-5-% 27 +\n U® n. Taylor. 


Für die Druckdifferenz zwischen Strahlmitte und Strahlrand ergibt sich daraus 


nach Prandtl 
. (33a). 


TER -U? nach Taylor 


Nach Prandtl besteht also zwischen Strahlmitte und Strahlrand keine Druckdifferenz, nach 
Taylor dagegen ein kleiner Überdruck. 


Um die Druckdifferenz gegenüber der ruhenden Luft zu erhalten, muß man berück- 
sichtigen, daß am Strahlrand ein Unterdruck gleich dem Staudruck der radial zuströmenden 
Luft herrscht. Mit p, (Druck der ruhenden Luft) ergibt sich für diese Druckdifferenz 


Pm— Pr = — 0,306 nach Prandtl 


. (33b). 
Pmn—m=— nach Taylor 


Es ergibt sich also in der Strahlmitte gegenüber der ruhenden Außenluft nach Prandtl ein 
kleiner Unterdruck, nach Taylor ungefähr Gleichdruck. 283 
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Einfluß der Bodentopographie auf Strömungen 


über der rotierenden Erde. 
Von H! Görtler in Göttingen. 


(Aus dem Kaiser-Wilhelm-Institut für Strömungsforschung.) 


Für den Einfluß von Bodenerhebungen auf den Verlauf der Strömungen im 
Meer und in der Atmosphäre wird eine quantitative Theorie der ersten 
Näherung für schwache Erhebungen entwickelt. Dabei wird der Einfluß 
der Sekundärströmung in der bodennahen Reibungsschicht mit erfaßt. 


Einleitung. 

Der Verlauf der Strömungen im Meer und in der Atmosphäre ist durch vielerlei Ein- 
flüsse bedingt (Erddrehung, Erdkrümmung, Bodengestalt, Veränderlichkeit der Dichte, Reibung, 
Wärmekonvektion, Strahlung, Feuchtigkeit der Luft, bei Meeresströmungen auch Windver- 
teilung an der Meeresoberfläche sowie Verdunstung und Niederschlag usw.). Einer allge- 
meinen Theorie begegnen daher nicht zu bewältigende Schwierigkeiten. Man ist gezwungen, 
sich mit der Untersuchung wesentlich vereinfachter Modelle zu begnügen, die nur einzelnen 
Einflüssen Rechnung tragen und die übrigen ignorieren. Diese Zwangslage führt aber anderer- 
seits gerade zu einem tieferen Verständnis der Rolle, welche die einzelnen Einflüsse beim 
Zustandekommen der Strömung übernehmen, wenn auch selbstverständlich das wirkliche 
Zusammenwirken dieser vielerlei Einflüsse dann im allgemeinen nicht einfach durch Super- 
position erschlossen werden kann. 

In der vorliegenden Arbeit wird die Frage nach dem Zusammenwirken von Erddrehung 
und Bodentopographie auf den Verlauf einer stationären (laminaren oder turbulenten) Strömung 
behandelt. Nur der Einfluß von Erdrotation, Bodenbeschaffenheit und Reibung steht zur 
Debatte, von allen anderen Einflüssen, auch von dem Einfluß der Erdkrümmung, wird abge- 
sehen. Es wird also die Frage aufgeworfen: Wie strömt eine Flüssigkeit konstan- 
ter Diehte und endlicher Höhe stationär über weit ausgedehnte Bodener- 
hebungen auf der rotierenden Erde (ersetzt durch ihre Tangentialebene) 
und zwar unter Berücksichtigung der Reibung und damit der durch den 
Ausfall der Corioliskräfte in der bodennahen Schicht in Gang gesetzten 
Sekundärströmung? 

Die Lösung dieses vereinfachten Problems kann zur Klärung der Anschauung über die 
Strömungserscheinungen sowohl im Meer als auch in einer stabil geschiehteten Atmosphäre — 
hier ersetzt durch die gleichförmige Atmosphäre mit freier Oberfläche — beitragen. Im 
Rahmen von Betrachtungen über atmosphärische Strömungen findet sich eine Erörterung 
unseres Problems unter Vernachlässigung der Reibung in dem Beitrag von L. Prandtl zum 
Meteorologenkongreß in Edinburgh 1936 [5]'). Für Meeresströmungen liegt in den verdienst- 
vollen Arbeiten von V. W. Ekman eine weitentwickelte Theorie vor, die neben der Boden- 
topographie auch die Erdkrümmung und den Einfluß des Windes (Erzeugung des Oberflächen- 
stromes) berücksichtigt. In seiner ersten Theorie [1] 1922 hat Ekman eine weitgehende 
Vereinfachung der hydrodynamischen Gleichungen vorgenommen, indem er die longitudinalen 
Trägheitskräfte gegenüber den Corioliskräften vernachlässigte. Dies ist aber selbst bei den 
kleinen im Meer vorkommenden Strömungsgeschwindigkeiten nicht zulässig. Ekman hat 
daher seine Theorie verbessert ([2] 1932) durch Mitnahme der zuvor vernachlässigten Träg- 
heitsglieder, mußte aber dafür in Kauf nehmen, daß dies wesentliche Erschwerungen der 
Rechnung und eine erheblich verwickeltere Darstellung zur Folge hatte, die das Vordringen 
zu expliziten Ergebnissen nur unter weiteren Einschränkungen und in speziellen Fällen 
gestatteten. 

Grundsätzlich stimmen die der folgenden Theorie zugrunde liegenden Anschauungen mit 
dem Standpunkt der neuen Ekmanschen Theorie der Meeresströmungen 1932 überein. Erst 
in der mathematischen Formulierung des Problems glauben wir einen Fortschritt erzielt zu 
haben. Wir gelangen zu Gleichungen, die einer expliziten Lösung zugänglich sind und zu 
sehr einfachen und in allen Phasen durchsichtigen Ergebnissen führen. Nach einer von der 
Ekmanschen Theorie unabhängigen Entwicklung und mathematischen Formulierung der all- 
zemeinen Grundgleichungen unserer Theorie (Kapitel I) wird der zweidimensionale Fall 
(Strömung über langgestreckte Gebirgskämme) in allen Einzelheiten gelöst und an einem 
Beispiel diskutiert (Kapitel II). Auf Grund unseres heutigen Wissens auf dem Gebiet der 
turbulenten Strömungen kann die Größenordnung der für den Strömungsvorgang charakteris- 
tischen Parameter für die Verhältnisse im Meer und in der Atmosphäre abgeschätzt werden. 


1) Die Zahlen in eckigen Klammern verweisen auf das Schrifttumverzeichnis am Schluß der Arbeit. 
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(Kapitel II). In Verbindung mit dieser Abschätzung lehrt die vorliegende Theorie, daß der 
Reibungseinfluß bei atmosphärischen und ozeanischen Strömungen für den Strömungsverlauf 
im allgemeinen von völlig untergeordneter Bedeutung ist. 

Auch die Behandlung des dreidimensionalen Falles räumlich beliebig verteilter Boden- 
erhebungen sowie die Einbeziehung des Erdkrümmungseinflusses dürfte bei der vorliegenden 
mathematischen Formulierung der Theorie ohne erhebliche neue Schwierigkeiten zugänglich 
sein. Die Untersuchung dieser Fragen wird zur Zeit in Angriff genommen. 


I. Voraussetzungen der Theorie und Ableitung der Grundgleichungen. 


$ 1. Die hydrodynamischen Grundgleichungen. Die auf ein mit der Drehgeschwindigkeit 
© rotierendes Koordinatensystem bezogenen Bewegungsgleichungen einer idealen Flüssigkeit 
lauten, wenn vd den Geschwindigkeitsvektor des Massenteilchens bezüglich des drehenden 
Systems, p den Druck, o die Dichte und ® ein Kräftepotential bedeuten, in Vektorschreib- 
weise zusammengefaßt: 


‚dv 
d 


Im Falle der Bewegung eines Flüssigkeitsteilchens auf der rotierenden Erdkugel setzt sich 
— grad ® additiv aus der Beschleunigung des Teilchens durch die Erdanziehung und der durch 
die Erdrotation hinzutretenden Zentrifugalbeschleunigung zusammen, stellt also die Schwere- 
beschleunigung 9 am Ort des Teilchens dar. Gegenüber dem ruhenden System ist ferner die 


Coriolisbeschleunigung 2 vXxo hinzugekommen. Die Winkelgeschwindigkeit » der Erddrehung 
hat den Betrag 


— grad grad p + 20x@. 


2n 
Während die Zentrifugalbeschleunigung als nur ortsabhängige Größe eine rein planetarische 
Wirkung der Erddrehung darstellt, äußert sich die von v abhängende Coriolisbeschleunigung 
in einer durch den jeweiligen Strömungszustand bedingten Weise. Ihrem Einfluß ist die 
vorliegende Untersuchung in erster Linie gewidmet. 
Zu den drei Eulerschen Impulserhaltungsgleichungen (1.1) tritt als vierte Bestimmungs- 
gleichung für die unbekannten drei Geschwindigkeitskomponenten und den Druck die Massen- 
erhaltungsgleichung 


7292.10’ 


$ 2. Koordinatenwahl. Den folgenden Rechnungen 

| werde ein rechtwinkliges kartesisches Koordinatensystem 
Ss x, 9,2 zugrunde gelegt. Sein Ursprung O0 möge mit irgend: 
einem festen Punkt der Erdoberfläche zusammenfallen. 
Während die z-Achse radial nach außen weisen soll, 
mögen die in der durch © gehenden Tangentialebene an 
die Erdkugel liegenden x- und y-Achsen so festgelegt 
sein, daß die positive «-Achse mit der von O aus in der 
Tangentialebene nach Norden weisenden Richtung den 
Winkel a bildet (Bild 1). An späterer Stelle soll die 
«-Richtung in besonderer Weise ausgezeichnet werden, 
wir halten daher zu diesem Zwecke den Winkel a frei 
verfügbar. Es sei ferner x die geographische Breite 
des Ursprungs 0. 


Die Komponenten des Drehvektors %) bezüglich der 
so definierten Achsenrichtungen lauten dann 


Bild 1. Koordinatenwahl. @;,=msing. 


Die entsprechenden Komponenten der Strömungsgeschwindigkeit d seien u,v, w. 


$ 3. Die Voraussetzungen der Theorie. 1 Wie bereits einleitend bemerkt, liegt die 
planetarische Wirkung der Erddrehung außerhalb des Interessenbereichs der vorliegenden 
Untersuchung. Wir beschränken uns in üblicher Weise auf eine gewisse Umgebung des 
Ortes 0, in welcher wir die geographische Breite durch ihren konstanten Wert in O und 
die Erdoberfläche in diesem Bereich durch ihre Tangentialebene z=0 in Q ersetzen. 
Dort ist dann ferner grad ® durch 
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zu nähern. 

® Eine in der Theorie der ozeanischen und atmosphärischen Strömungen angebrachte 
Vernachlässigung beruht auf der verhältnismäßig geringen Höhe der betreffenden Flüssig- 
keitsschichten im Vergleich zu ihren horizontalen Ausdehnungen. Da die vertikalen Be- 
schleunigungen nämlich auf nur kleinen Wegen wirksam sind relativ zu den langen Strecken, 
auf denen die horizontalen Beschleunigungen erfolgen, vernachlässigt man die ersteren und 
zwar einschließlich der vertikalen Coriolisbeschleunigung 2(uo,—v®,). In den strengen 
Eulergleichungen (1.1) für stationäre Strömungen in Komponentenschreibweise: 


du du du 107 
u. oda +20 {vsing — wcosgpsina}, 
w 10 
COSY Cosa} 


sind dann die mit darunter stehender Ziffer 2 versehenen Glieder zu streichen, womit sich 
die dritte Eulergleichung auf die statische Grundgleichung reduziert. 

® Wir wollen ferner annehmen, daß die Höhendifferenzen der Bodenerhebungen im 
betrachteten Strömungsbereich klein seien gegenüber der Höhe H der Flüssigkeitsschieht. Ist 
etwa 2={(x,,) der Boden und an einer Stelle «,, 4, des Bereichs (etwa an einer Randstelle 
an der Einströmungsseite) 9%) =0 definiert, so soll im ganzen Bereich |{|<H sein. 
Ist die Länge B ein Maß für die horizontale Ausdehnung der betrachteten Bodenunebenheiten, 


so ist, insofern es sich um einen „glatten“ Berg handelt, daun |grad £| <5- Wir wollen dies 


aber nicht fordern, sondern nur die schwächere Voraussetzung machen, daß grad {|<1 sei. 
Damit ist dann die Vertikalgeschwindigkeit »w im Strömungsbereich klein gegenüber der 
horizontalen Geschwindigkeit u,v. Das hat zur Folge, daß in den Eulergleichungen (3.1) die 
mit darunter stehender Ziffer ® gekennzeichneten Glieder ebenfalls zu vernachlässigen sind. 

Die Ergebnisse der Rechnung werden die beiden letztgenannten Vernachlässigungen 
für alle interessierenden Strömungen über flache Bodenerhebungen sowohl bei ozeanischen 
wie bei atmosphärischen Verhältnissen bestätigen. 

Die verbleibenden Gleichungen für «,», w,p lauten nun: 


du 10» 


dv 
u, 
(3.2). 
du 
dx dy d2 


Hierin ist zur Abkürzung 
sin @ (= const) 


gesetzt. Die Abhängigkeit von der Wahl von «a ist in Fortfall gekommen. 

Wir werden im folgenden immer annehmen, daß die Strömung bei Eintritt in den 
betrachteten Strömungsbereich eine von z unabhängige Verteilung der Horizontalgeschwindig- 
keit u,» besitzt, in diesem Sinne an der Eintrittsstelle also zweidimensional ist. Das ist 
etwa der Fall, wenn die Strömung als reine Parallelströmung über ebenem Boden ankommend 
den Strömungsbereich der Bodenerhebung erreicht. Wir erkennen aber sofort, daß die 
Strömung bei einer solchen Einströmungsbedingung dann im ganzen Strömungsbereich diesen 
Charakter bewahrt. Die GlIn. (3.2) werden dann nämlich gelöst durch Ausdrücke u,v, w 
der Gestalt 


| 
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u=u(x,Y), vr=v(e,Yy), 


Stromlinien, die an einem Orte durch eine Parallele zur z-Achse gehen, tun dies also an 
jedem Orte des Strömungsbereiches, und die Aufwärtsgeschwindigkeit hängt linear von z ab. 
In der Tat, nehmen wir der ersten Zeile von (3.3) gemäß «u und v als von z unabhängig an, 
so folgt die zweite Zeile von (3.3) mit unbestimmten a,b aus der Kontinuitätsgleichung (1.2), 
2m 

weil dann 0 sein muß. Die gemachte Annahme rechtfertigt sich aber, indem nun « 
und » wie folgt als Funktionen von « und ,y allein bestimmt werden können: Durch Elimi- 
nation von p aus den beiden ersten Eulergleichungen (3.2) erhält man eine erste, und aus 
der Kontinuitätsgleichung mit za (x,9) eine zweite Differentialgleichung für « und v 
allein bezüglich der Veränderlichen x und y. Mit Hilfe der Randbedingungen bestimmen 
sich dann «, » und w in der Gestalt (3.3). Den Druck p gewinnt man aus den Euler- 
gleichungen durch Integration in der Gestalt p=—09z2-+p,(x,y). Die hier geschilderten 
Rechenschritte werden in der nachfolgenden Theorie durchgeführt werden. 

Es ist hier noch folgende Bemerkung am Platze. Hat man es nicht mit dem vorstehend 
geschilderten Problem zu tun, sondern mit einer Strömung über eine sich in ihrer Ebene 
drehende Scheibe, so fallen in den beiden ersten Eulergleichungen (3.1) gerade dieselben 
Glieder weg, die wir oben vernachlässigt haben (hier, weil &.=»,=0, oben dagegen, weil 
w klein sein sollte). Macht man auch hier die Annahme, daß die vertikalen Beschleunigungs- 
glieder zu vernachlässigen sind, und wählt man die Drehgeschwindigkeit =o'’=wsing, so 
hat man es also mit denselben GIn. (3.2) zu tun, die auf (3.3) führen. Man kann dann im 
Laboratoriumsversuch die Ergebnisse der Theorie mit befriedigender Genauigkeit nachprüfen. 
So waren uns die in der ersten Zeile von (3.3) zum Ausdruck kommenden Verhältnisse 
(Unabhängigkeit der Komponenten « und v von 2) aus Strömungsversuchen im rotierenden 
Laboratorium des Kaiser-Wilhelm-Instituts für Strömungsforschung, Göttingen, geläufig. Man 
kann sie durch Farbbahnen in einer mitrotierenden flüssigkeitserfüllten Schale bestätigen. 

Im Hinblick darauf, daß die Horizontalkomponenten « und ® eine zweidimensionale 
Geschwindigkeitsverteilung darstellen, zerlegen wir den Geschwindigkeitsvektor d zweck- 
mäßigerweise in 


mit B= fu, », 0}, 
w={0,0,"} 


und behandeln diese Vektoren gesondert. Um Verwechslungen zu vermeiden, schreiben wir 
im Gegensatz zu den dreidimensionalen Differentialoperatoren die entsprechenden zweidimen- 
sionalen mit großen Anfangsbuchstaben: 


Grad 
‚day 
Rota= (skalar!) 


Wir bringen nun noch zwei bekannte Tatsachen in Erinnerung. Soll eine Flüssigkeit 
gezwungen werden, über ebenem horizontalem Boden mit konstanter Geschwindigkeit 
v=WB,=fu,,v,) zu strömen, so müssen die durch ®, erzeugten Coriolisbeschleunigungen 
mit einem entsprechenden Druckgefälle im Gleichgewicht stehen. Aus (3.2) ergibt sich durch 


Integration 


demzufolge muß die freie Oberfläche (p=const) die bekannte Spiegelneigung besitzen. 
Jeder freien oder erzwungenen Bewegung entspricht vermöge der entstehenden Coriolis- 
kräfte ein seitwärts gerichteter Drang. Bei freien waagrechten Bewegungen ergibt sich eine 
Ablenkung (nach rechts auf der Nordhalbkugel, nach links auf der Südhalbkugel), deren 
Größe man unter Heranziehung des Begriffs des „Pendeltages“ zu veranschaulichen pflegt 


(a. i. die Zeit, die ein Foucaultsches Pendel braucht, um eine volle Drehung von 360° 


zurückzulegen, also a), Ein Massenpunkt beschreibt dann in einem halben Pendel- 
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tag einen vollen Kreis, den „Trägheitskreis“ (insoweit die Veränderlichkeit von 9 bedeutungs- 
los bleibt; in der Äquatorgegend ergeben sich mit Berücksichtigung dieser Veränderlichkeit 
in erster Näherung Kurven vom Typ der Elastica, s. L. Prandtl [4]). Der Radius dieses 
Trägheitskreises ist also 


wenn V, der konstante Betrag der Geschwindigkeit des Massenpunkts ist. 


(3.7), 


$ 4. Berechnung der Aufwärtsgeschwindigkeiten einer Strömung über unebenem rotieren- 
dem Boden unter Berücksichtigung der Sekundärströmung in der bodennahen Reibungsschicht. 
Gegenüber den bisher geschilderten Verhältnissen treten ganz neue Erscheinungen auf, wenn 
man die Reibung berücksichtigt. Die Reibungswirkung entfaltet sich in der relativ dünnen 
bodennahen Grenzschicht, in der die Strömungsgeschwindigkeit wegen des Haftens am Boden 
auf Null absiukt. Das Druckfeld in der darüberliegenden Hauptströmung, in der die Reibung 
vernachlässigt werden kann, setzt sich, wie man aus den Grundgesetzen der Grenzschicht- 
theorie weiß, praktisch unverändert durch die Grenzschicht fort. Hier in der Grenzschicht 
gehen aber die Geschwindigkeiten gegen Null, die ihnen entsprechenden Corioliskräfte können 
nicht mehr denselben Anteil zum Kräftegleichgewicht wie in der Hauptströmung beitragen. 
Demzufolge wird in der Grenzschicht eine Sekundärströmung in Gang gesetzt. Diese verläuft 
umgekehrt zur Richtung der ausfallenden Corioliskräfte und damit in der Grundrißprojektion 
senkrecht zur Hauptströmung. Sie erfolgt in solcher Stärke, daß die durch sie hervorgerufenen 
Reibungskräfte am Boden das Gleichgewicht gerade wiederherstellen. 

Liegt etwa eine geradlinig-gleichförmige Hauptströmung vor, so bleiben die Vorgänge 
in der Grenzschicht ohne Einfluß auf die Hauptströmung. Es erfolgt in der Grenzschicht 
durch die zusätzliche dazu senkrechte Sekundärströmung lediglich insgesamt eine 
Strömung, die schräg zur Hauptströmung verläuft. Liegt jedoch etwa eine nicht 
drehungsfreie Hauptströmung vor, so läßt schon die Anschauung erwarten, daß die horizontale 
Sekundärströmung, da sie nicht mehr divergenzfrei sein wird, aus Kontinuitätsgründen nur 
dann bestehen kann, wenn ein entsprechender vertikaler Massenzufluß von bzw. zu der 
Hauptströmung erfolgt. So kann die Sekundärströmung in der relativ dünnen Bodenschicht 
unter Umständen eine stark eingreifende Wirkung auf die Hauptströmung ausüben. Es ist 
daher von praktischem Interesse, den Einfluß der Reibung zu untersuchen. Wir wollen zu- 
nächst zeigen, wie der Reibungseinfluß für die Rechnung zweckmäßig erfaßt werden kann. 
Wir werden dabei zugleich die eben gemachten mehr anschaulich-orientierenden Bemerkungen 
in präzisierter Form bestätigt finden. 

Ist v, der waagrechte Geschwindigkeitsvektor der Sekundärströmung in der Grenzschicht, 
so ist der durch diesen bewirkte sekundliche Massenstrom M durch 


gegeben, wo Ö die Grenzschichtdicke und 2={ (x, y) den Boden bezeichnen. Zur Erhaltung der 
Masse muß eine Vertikalgeschwindigkeit mw, vermöge 


Om; 
Div Is =( 


vorhanden sein, also 


wo »,=[w;]|=:;+9s die durch die Sekundärströmung in Gang gesetzte Aufwärtsbewegung am 
Rande der Grenzschicht darstellt ([w,];=:—=0 wegen des Haftens am Boden). 

Die Grenzschichtströmung am Boden können wir uns nun in folgender Weise ersetzt 
denken: Die Hauptströmung fließt reibungslos über einen um die Verdrängungsdicke der 
Grenzschicht nach oben versetzten Boden. Die von der Sekundärströmung in der Grenz- 
schicht zu liefernde, jeweils von der Hauptströmung abhängige Vertikalgeschwindigkeit an 
diesem Ersatzboden kann man sich dann als von einer geeignet dort anzubringenden Quellen- 
und Senkenverteilung herrührend denken. Unter 2={(x,y) soll im weiteren dieser so gegen 
den wahren Boden vesetzte neue Boden verstanden werden. 

Unter Berücksichtigung der Richtung des Massenflusses schreiben wir 


| 
XuM 
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wo B={-—v,u} aus ® durch Drehung um einen rechten Winkel nach links hervorgeht und 
h einen Proportionalitätsfaktor von der Dimension einer Länge bezeichnet, dessen Vorzeichen 
durch sign h = sign sin @ festgelegt ist. Wir werden, um Vorzeichenunterscheidungen zu ver- 
meiden, unsere Aussagen im folgenden immer für die nördliche Halbkugel formulieren 
(h>0). hist die Dieke jener Flüssigkeitsschicht, durch die mit der festen Geschwindigkeit 
® dieselbe Menge pro Zeiteinheit befördert würde, wie jene, die in Wirklichkeit durch die 
Sekundärströmung am betreffenden Ort in der Grenzschicht horizontal versetzt wird. 

Bei laminarer Strömung über ebenem Boden mit konstanter Geschwindigkeit ®,—= 
fu,, v,, außerhalb der bodennahen Reibungsschicht läßt sich die Geschwindigkeitsverteilung 
berechnen (siehe [6]), und es kann in elementarer Weise gefolgert werden, daß in diesem Falle 


Zahl 


also von der jeweiligen Geschwindigkeit ®, unabhängig ist. Dasselbe gilt natürlich auch 
für entsprechende turbulente Strömungen insoweit, als man in üblicher Weise, den erhöhten 
Austausch durch eine vergrößerte konstante Zähigkeit in Rechnung setzen kann. Für unsere 
Theorie des Einflusses schwacher Bodenerhebungen genügt es für eine erste Näherung, auch 
hier h als geschwindigkeitsunabhängig anzunehmen. Streng aber muß man h als von ® 
abhängig, und also bei gegebenem ®(x,y) als h=h(x, y) annehmen. Um uns die Möglichkeit 
dieser Verbesserung offenzuhalten (siehe $ 11), sehen wir h im folgenden als veränderlich an. 

Es ist 
DvM=oDiv(hV)= —-oRot(h B). 


Somit ist die von der Quellen- und Senkenbelegung geliefert zu denkende Aufwärtsgeschwindig- 


keit am Boden 


Fließt nun die Hauptströmung über den durch z2={ (x, y) gegebenen Boden, so erfährt sie 
hierdurch die Aufwärtsgeschwindigkeit 


Somit ergibt sich für die Aufwärtsgeschwindigkeit w in der ER mit Rücksicht 
auf (3.3) insgesamt: 


wo H(x,»y) die freie Oberfläche der Strömung und my (x, y) die Vertikalgeschwindigkeit an 
dieser Oberfläche bezeichnen. Die Rechnung wird bestätigen, daß man für |Ü|<H in einer 
befriedigend großen Umgebung von O die freie Oberfläche H in allen uns interessierenden Fällen 
durch ein mittleres konstantes // ersetzen kann. Wir werden dies noch ausführlich diskutieren 
(siehe $ 14). Wir haben dann 


$ 5. Ableitung der Ausgangsgleichungen der weiteren Theorie. Wir haben nun alle Vor- 
bereitungen getroffen, um die grundlegenden Gleichungen der nachfolgenden Theorie zu for- 
mulieren. Die beiden ersten Eulergleichungen (3.2) liefern für die Hauptströmung nach 
Elimination des Druckes die Beziehung 


BoGrad Rot B+DivB(Rt B+2W0)=0.. N), 


und die Kontinuitätsgleichung ergibt mit Rücksicht auf (4.7) und wegen Div” = -. 


. . 2.222.082). 


Diese beiden Gleichungen für « und v bilden die Grundlage der weiteren Rechnung, in 
der zunächst einmal der Fall der Strömung über langgestreckte Gebirgskämme behandelt 
werden soll, für den eine Erniedrigung der Dimensionszahl möglich ist. ®, das im allgemeinen 
Fall von x und y abhängt, wird dann nur von einer einzigen Veränderlichen abhängen, wo- 
mit (5.1) und (5.2) sich auf gewöhnliche Differentialgleichungen für « und v reduzieren. Sind 
« und v bestimmt, so kann man nachher w und p in elementarer Rechnung ermitteln, und 
zwar w am bequemsten aus 


 ° 
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Es sei nun noch bemerkt, daß man zunächst versucht sein wird, die Gln. (5.1) und (5.2) 
für sehr kleine Geschwindigkeiten ® zu behandeln. Dazu würde man (5.1) unter Streichung 
der in ® und dessen Ableitungen quadratischen Glieder durch 


Dv®=0 
Roth®=— BoGrad? 


erhalten. Physikalisch bedeutet dies die Vernachlässigung der longitudinalen Trägheitsglieder 
gegenüber den in den Geschwindigkeiten linearen Coriolisgliedern in den Eulerschen 
Gleichungen. V. W. Ekman ist ursprünglich in seiner Theorie der Meeresströmungen 1922 
sinngemäß von diesem Gedanken ausgegangen (s. [1]), da bei den kleinen Geschwindigkeiten 
im Meer die vorkommenden Trägheitskreisradien klein sind im Vergleich zu den in Frage 
kommenden horizontalen Abmessungen der Bodenerhebungen. Allein es hat sich gezeigt — 
und dies wird auch in der vorliegenden Untersuchung klar hervortreten —, daß selbst bei 
Meeresströmungen diese Trägheitsbeschleunigungen eine wesentliche Rolle spielen und die 
Strömungen gänzlich von denen verschieden sind, die bei einem bedeutenden Überwiegen 
der Coriolisbeschleunigungen resultieren würden. Selbst eine iterative Verbesserung des oben 
skizzierten Näherungsansatzes reicht nicht aus, um zu praktisch interessierenden Ergebnissen 
zu gelangen. Unsere Rechnungen werden also keine diesbezüglichen Voraussetzungen machen, 
sondern auf den vollständigen Gln. (5.1) und (5.2) gegründet sein. 

Wesentliche Voraussetzungen unserer Theorie werden neben dem 
bewußten Absehen von planetarischen Einflüssen (Ersatz der Erdober- 
fläche durch ihre Tangentialebene und der variablen Breite durch einen 
mittleren konstanten Wert) also allein die relative Flachheit der Strö- 
mungsräume und in diesen Räumen die relative Kleinheit der Boden- 
erhebungen £ gegenüber der Schichthöhe AH sein. 


ersetzen und somit nach (5.2) 


Il. Vollständige Lösung des Problems für den zweidimensionalen Fall. 


In diesem Kapitel II soll zunächst die allgemeine mathematische Berechnung und Dis- 
kussion des zweidimensionalen Falles durchgeführt werden, worauf in Kapitel III die Ergeb- 
nisse speziell auf Strömungen im Meer und in der Atmosphäre angewendet werden sollen. 


$ 6. Vereinfachung des Problems im zweidimensionalen Fall. Wir nehmen an, daß die 
Bodenfläche { und damit die Geschwindigkeit ® von y unabhängig sind. Es bezeichne 
B,={u,,v,, die Strömungsgeschwindigkeit an der Eintrittsstelle = x, der Strömung in den 
betrachteten Strömungsbereich 2, =x2=<x,, in welchem sich die uns interessierende Boden- 
erhebung befindet. Wir setzen fest, daß £ (x,)=0 sein soll, was lediglich die Festlegung des 
Niveaus der x,y-Ebene bedeutet. Ferner soll, um hydrodynamisch saubere Bedingungen zu 
schaffen, angenommen werden, daß das Medium für << x, über ebenem Boden ströme. 
Hier ist also auch noch ungestört B=%2,. Es ist £ (x) als überall stetig und stetig-differenzier- 


bar vorausgesetzt, und es soll also auch z (z,)—=0 sein. In «x, (oder weiter stromabwärts) 
beginnt die Bodenerhebung oder -senkung. 
Die Gl. (5.1) vereinfacht sich nun auf 


also nach einmaliger Integration unter Beachtung der Randbedingungen 


Hinzu kommt die Gl. (5.2), welche nun in der Gestalt 


geschrieben werden kann und integriert 


(c, Integrationskonstante) ergibt. 
Wie bereits in $ 4 bemerkt, ist die für den Reibungseinfluß maßgebende Länge h im 
allgemeinen abhängig von der Geschwindigkeit ®. Unsere Kenntnisse über die Gesetze turbu- 
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lenter Strömung reichen aus, um für diese Abhängigkeit wenigstens eine Abschätzung zu 
geben. Dies soll an späterer Stelle (Kapitel III) geschehen. Zunächst aber denken wir uns 
die Länge h durch einen konstanten mittleren Wert im Strömungsbereich ersetzt. Mit h = const 
erzielen wir einfache Formeln. Es wird Aufgabe einer dann folgenden Untersuchung ($ 11) 
sein, den Einfluß der Veränderlichkeit von h auf die Ergebnisse der Rechnung zu untersuchen. 
Es wird sich zeigen, daß unser vorläufiger Ersatz von h durch einen Mittelwert in weiten 
Grenzen eine gute Näherung gibt. Wenn wir hier vorwegnehmen, daß in Verbindung mit 
den folgenden Rechnungen die oben angekündigte Abschätzung des Kapitels III den Einfluß 
der Reibung als eine bei ozeanischen und atmosphärischen Strömungen im allgemeinen be- 
langlose Größe nachweisen wird, wird der Leser den Ersatz von k durch einen Mittelwert 
mit dem Ziele einer ersten Näherung billigen. 

Mit h=const wird in (6.4a) wegen ®=®%, und (=0 in <=, die Integrationskonstante 

— Hu —hv,, also 


Aus den Gln. (6.2) und (6.4) mit kA = const stellen wir eine Differentialgleichung für die 


Geschwindigkeitskomponente « allein her. Indem wir aus beiden Gleichungen 1 berechnen, 


erhalten wir nach Gleichsetzen der gewonnenen Ausdrücke 


Hat man aus dieser nicht-linearen Differentialgleichung « ermittelt, so ergibt sich durch Ein- 
setzen in (6.4) sofort auch v. Das Problem ist alles in allem auf eine mathematisch ver- 
hältnismäßig einfache Fragestellung zurückgeführt worden. 

Die Behandlung der verbleibenden Differentialgleichung (6.5) wird PR nun noch wesent- 
lich erleichtert durch folgende Bemerkung. Daß wir gerade für die Komponente « allein eine 
Differentialgleichung herstellten, hat seinen Grund darin, daß dies dem Problem besonders 
angepaßt ist. Es wird (siehe $ 8 und allgemein für variables h auch $ 11) leicht einzusehen 
sein, daß mit |Zly. <H die Störung von « durch die Bodenerhebung, d. h. die Differenz 
ua — u,, dem Betrage nach klein gegenüber «, von mindestens derselben Ordnung ist. Eine 
bezüglich £/H linearisierte Theorie, wie sie in unseren bisherigen allgemeinen Entwicklungen 
des ersten Kapitels vorliegt, hat also konsequenterweise stets auch eine Linearisierung bezüg- 
lich der relativen «-Störung zur Folge. Es wäre dagegen unzweckmäßig gewesen, eine (nicht- 
lineare) Differentialgleichung für v allein herzustellen, denn die Störungen » — v, bleiben in 
den physikalisch interessierenden Fällen im allgemeinen durchaus nicht klein gegenüber dem 
Betrag der ungestörten Anströmgeschwindigkeit ®,. 


$ 7. Einführung dimensionsloser Größen. Bevor wir fortfahren, wird es zweckmäßig 
sein, die auftretenden Veränderlichen auf charakteristische Größen zu beziehen, um dann in 
dimensionsloser Schreibweise bequemer die Größenordnung der sich ergebenden Ausdrücke 
überblicken zu können. 

Es sei B eine charakteristische Länge, durch welche die waagrechte Breitenerstreckung 
der betrachteten, in «=, beginnenden Bodenerhebung und damit des zugrunde liegenden 
Strömungsbereichs in «-Richtung definiert ist. Auf diese beziehen wir alle horizontalen Längen- 
messungen als Einheit. Die vertikalen Abmessungen beziehen wir dagegen auf A. Ferner 
erweist es sich als zweckmäßig, alle Geschwindigkeiten auf die zum Bergkamm senkrecht 
gerichtete Komponente «u, der ungestörten Anströmungsgeschwindigkeit in =x, zu beziehen. 

Um die Schreibweise einfach zu halten, wollen wir vereinbaren, daß für alle folgenden 
Ausführungen unter «*, y*, z*, H*, w*, die bisher benutzten dimensions- 
behafteten Größen zu verstehen sind, während für das Weitere die entsprechenden nicht 
mit einem Stern versehenen Größen ausschließlich die folgenden Dimensionslosen sein sollen, 
mit denen fortan gerechnet wird: 


y* 
H* 


u* v* w* 
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Die Größen B, H und h bleiben nach wie vor dimensionsbehaftet. Die Dimensionslosen x 
und z sind im interessierenden Bereich des Berges und seiner Umgebung von der Größen- 
ordnung O(1). Nach Voraussetzung ist ferner |{!<1. Behalten wir dies im Auge, so ist es 
leicht, sich in jedem Stadium der folgenden Rechnung einen Überblick über die auftretenden 
Größenordnungen zu verschaffen. 

Als charakteristisch für die hier behandelten Vorgänge der Stromablenkung durch Boden- 


erhebungen auf der rotierenden Erde erweisen sich die Dimensionslosen MT und 5 wo 


den zu u, gehörigen Trägheitskreisradius bezeichnet. (Damit hängt es zusammen, daß man 
gerade mit der «-Komponente der Anströmgeschwindigkeit als Bezugsgeschwindigkeit eine 
einfache Schreibweise der dimensionslosen Gleichungen erzielt.) 

Neben den eben genannten Dimensionslosen verwenden wir noch ihr Verhältnis, die Zahl 


Dieser Zahl @ wird im weiteren eine wesentliche Rolle zufallen. 

Die Bestimmungsgleichungen (6.5), (6.4) und (5.3) für die Geschwindigkeitskomponenten 
schreiben sich in den eingeführten Dimensionslosen wie folgt: 


R 


d 
H du 


Nebenbei, der Standpunkt der älteren Ekmanschen Theorie der Meeresströmungen aus 
dem Jahre 1922 [1], nämlich: Vernachlässigung der longitudinalen Trägheitsbeschleunigungen 
gegenüber den Coriolisbeschleunigungen, würde in unseren Gleichungen dadurch zum Aus- 
druck kommen, daß R=0 und also auch @=0 zu setzen wäre. In diesem Grenzfall ergibt 
sich nach unserer Theorie bereits das fertige Resultat: 


und als Projektion der Stromlinien der Hauptströmung auf die «, y-Ebene demnach 
x 
H 
xo 


Diese Stromlinien haben über den Extrema von £ (Bergkamm, Talboden) Wendepunkte, wie 
dies V. W. Ekman [1], S. 109, auch nach seiner damaligen Theorie am speziellen Beispiel 
eines sinusförmig gewellten Bodens in einer Figur demonstrierte. Die für normale Meeres- 
strömungen und atmosphärische Strömungen anzunehmenden Werte von @ sind aber nicht 
klein, sondern im allgemeinen sogar so groß, daß die resultierenden Strömungen in ihrem 
Verlauf einen vom Falle @=0 gänzlich abweichenden Charakter haben. Dieser wird aber 
durch einen anderen Grenzfall, nämlich den der reibungslosen Strömung (h=0, also G=»), 
genähert, wie wir noch sehen werden. 


$ 8. Geschwindigkeitsverteilung und Stromablenkung. Für das in x stetige und stetig- 
differenzierbare Bodenprofil £ (x) bestehen auf Grund unserer allgemeinen Theorie des ersten 


Kapitels die Beschränkungen, daß im Strömungsbereich |{|<1 und ferner i<z sein 


soll (letzteres als Garantie für das Kleinbleiben der Vertikalkomponente w der Strömungs- 
geschwindigkeit). Bei den ins Auge gefaßten Anwendungen ist B/H eine gegen 1 sehr große 
Zahl (z. B. 10°), wie dies ja auch der gemachten Voraussetzung relativ flacher Strömungs- 
räume entspricht. Für die folgenden Betrachtungen bedeutet daher die genannte Beschränkung 


| 
von Fa keine weitere Vereinfachung. Glücklicherweise kommt man aber auch ohne schärfere 


Beschränkung von a2) aus, da sich im folgenden die Kleinheit von |{| bereits ohne Annahmen 
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über - als für die Durchführbarkeit der Rechnung ausreichend erweisen wird. Es werden 


also im weiteren keinerlei Voraussetzungen über Größenordnungen gemacht, die nicht schon 
der allgemeinen Theorie des ersten Kapitels zugrunde liegen. Im übrigen würde eine der 
Forderung an £ entsprechende For- 


7 
derung an also eine 
| ' den Gegebenheiten der geplanten 
Anwendungen nicht entsprechen- 
de Beschränkung bedeuten (siehe 
Bild 2a). Die von uns gemachte, 
wesentlich schwächere Voraus- 


setzung: wird diesen je- 


doch schon eher gerecht (siehe 
Bild 2b, beachte aber, daß der 
-—-r wirkliche Strömungsraum im Ver- 
hältnis B:H gestreckt ist, d. h. 


Vorbehaltlich der Bestätigung durch die weitere Rechnung nehmen wir nun, wie bereits 
in $ 6 angekündigt, an, daß die Störung der «-Komponente der Strömungsgeschwindigkeit 
durch die Bodenerhebung, d.h. also die Differenz « — 1, dem Betrage nach klein gegen 1 ist 


von (mindestens) derselben ersten Ordnung, wie nach Voraussetzung |{!y., Klein im Vergleich 
zu 1 ist. Wir schreiben 


Die Bezeichnung „Z;* für die dimensionslose Störungsgröße ist einer später folgenden Deutung 
derselben bereits im voraus angepaßt. 
Gehen wir mit diesem Ansatz (8.1) in die Differentialgleichung (7.4) für « ein, so er- 
gibt sich: 
dd; 


Bild 2a.° <1. Bild 2b. didx 


Für eine Berechnung der Störungsgröße {; in erster Ordnung — nur diesen Anspruch kann 
unsere Theorie überhaupt stellen — erhalten wir somit wegen |{|<1 und unter der oben 
ausgesprochenen Annahme, daß daher auch |{;!<1 ist, die folgende lineare Differential- 
gleichung erster Ordnung: 


Berücksichtigt man die Randbedingung, daß in «=, die Störungsgröße {;=0 sein muß (dort 
ist v=1), so erhält man für Z; als Lösung: 


Durch partielle Integration folgt: 


X 
Die Bodenerhebung bedeckt nach Vereinbarung eine Strecke <—x,=0(l). Man erkennt 
ohne weitere Erörterung, daß die bezüglich der Größenordnung von £; gemachte Annahme 
bestätigt wird, und zwar für alle möglichen Werte von 6 (0). 
Indem wir die G]. (8.2) als Differentialgleichung für £ (x) auffassen, erhalten wir als Um- 
kehrung von (8.4) die weitere wertvolle Beziehung 


x 


xo 
Verwenden wir (8.6) gleich bei der Berechnung von » aus der — bezüglich £ und {; 
entsprechend zu linearisierenden — Gl. (7.5), so ergibt sich insgesamt folgende Geschwindig- 
keitsverteilung: 
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(8.6). 


at, 
x; 


Es erhellt unmittelbar, daß die Störung ® — v, der v»-Komponente im allgemeinen durch- 
aus nicht klein gegen 1 ist. Bei Meeresströmungen etwa interessiert gerade der Fall großer 
Werte von B/R. 


Die Stromlinien der Hauptströmung ergeben in ihrer Projektion auf die x, y-Ebene ein 


Bild von der horizontalen Ablenkung der Hauptströmung. Diese Projektion berechnet sich 
x 


v IR 
aus y—y,=\ bei Linearisierung bezüglich zu 
y-y 


xo 


x 
Xo 


B x Jq . . . . . . . . . (8.7). 
x) 


Der von diesen Kurven mit der positiven «-Achse gebildete Winkel y berechnet sich aus 
x 


xo 
Ist tgy,=r,, so ist also die Ablenkung y— y, gegenüber der Anströmungsrichtung gegeben 
durch 


x 


x 
B Br 
... 
x 
Der Grenzfall verschwindender Reibung wird durch h>0, G >» gegeben. 
Es ist 
h>V 


Die Gln. (8.6) bis (8.9) sind bewußt in einer solchen Form geschrieben, daß man diesen 
Grenzfall unmittelbar erhält, indem man einfach dureh ersetzt. 

Das führt, umgekehrt aufgefaßt, zu einer wertvollen Deutung der Größe L;: Es ist 
&; jenes gedachte Bodenprofil, durch welches man das wirkliche Profil £ er- 
setzen muß, wenn man bei entsprechender Rechnung ohne Berücksichti- 
gung der Reibung dieselbe Hauptströmung erhalten will. Diese Deutung von £; 
liefert uns in der Abweichung der Größe £, von £ direkt ein anschauliches Maß für den 
Reibungseinfluß. 

Wenn man die Antwort auf die obige Fragestellung nach jener Bodenerhebung, welche 
bei Rechnung ohne Berücksichtigung der Reibung dieselbe Hauptströmung liefert wie der 
wirkliche Berg £ (x) mit Berücksichtigung der Reibung, nicht direkt aus (8.6) nebst (8.10) ab- 
liest, sondern auf Grund der Gl. (5.2) zu beantworten sucht, so ist zu fordern, daß der Ersatz- 
berg — nennen wir ihn jetzt vorübergehend & — der Gleichung 


VoGradö= Vo Grad Rot ® 


genügt. Im vorliegenden zweidimensionalen Fall vereinfacht sich diese auf 
dz ulld«' 
Es ist also die Steigung des Ersatzberges dort kleiner (größer) als die wahre Bergsteigung, 
wo Rot (h »-n4"<0> 0) ist, im Einklang mit unserer Quellen-Senken-Vorstellung, wo- 


nach der abgeänderte Berg £ die verringerte (vermehrte) Aufwärtsgeschwindigkeit in Boden- 
nähe veranlassen muß, welche durch das zusätzliche w, hereinkommt. Ersetzt man in obiger 
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Gleichung « und ® durch die berechneten Werte (8.6), wobei für unsere Theorie der 1. Ord- 
nung hier « im Nenner durch 1 zu ersetzen ist, so ergibt sich als Bestätigung natürlich 
wieder, daß :{; sein muß. 

Der sich im reibungslosen Grenzfall ergebende Ausdruck für v stimmt, wie es sein muß, 
mit dem von L. Prandtl [5] angegebenen Ergebnis überein, welches dort einfach auf Grund 
des Helmholtzschen Wirbelsatzes erschlossen wurde. (Zum reibungslosen Fall, dessen 
Theorie, selbständig durchgeführt, natürlich wesentlich einfacher wird, als sie in der obigen 
Behandlung als Grenzfall in Erscheinung tritt, sei der Leser auch auf den Schluß unseres 
$ 15 verwiesen.) 

Für die numerische Rechnung ist die Gestalt der Gln. (8.6) bis (8.9) mit den darin auf- 
tretenden Integralen über [; weniger geeignet als eine solche, in der durch Verwendung von 
(8.5) diese Integrale beseitigt sind. So erhält man etwa durch Umformung 


Zwar ist in dieser Form der Grenzfall verschwindender Reibung nicht mehr unmittelbar 
ersichtlich, dafür aber tritt ein anderer Grenzfall, nämlich der Grenzfall sehr kleiner Ge- 
schwindigkeiten, einfach zutage: Sind die Geschwindigkeiten so klein, daß man die longitudi- 
nalen Trägheitsglieder gegenüber den Coriolisschen vernachlässigen kann, so ist in obigen 
Gleichungen R=0 zu setzen (siehe $ 5). Damit ist auch @=0 und es ist 
Wie es sein muß, erhält man für diesen Grenzfall die bereits in (7.7) und (7.8) angegebenen 
Gleichungen für Geschw indigkeitsverteilung und Stromlinienprojektion. 

Für große Werte von @ ist es angebracht, den Ausdruck (8.4) in eine Reihe nach 
fallenden Potenzen von @ zu entwickeln. Indem man die Entwicklung der Exponential- 
funktion im Integranden vornimmt, ergibt sich die für alle x und @=F0 absolut und gleich- 
mäßig konvergente Reihe: 


| 


x 
Von dieser Reihenentwicklung sind jeweils nur so viele Glieder zu berücksichtigen, als dies 
für eine Berechnung von {; in erster Ordnung notwendig ist; Beiträgen, die dem Betrage 
nach klein im Vergleich zu |Z/ya, sind, kommt im Rahmen unserer Theorie der ersten Ord- 
nung natürlich keine Bedeutung zu. 

Die gleiche Entw icklung erhält man übrigens durch einen Iterationsprozeß in (8.5), aus- 
gehend von der Näherung {;={ für G=x. Man erhält sie aber auch einfach durch wieder- 
holte partielle Integration von (8.3). Beidemal ist nur noch zur Identifizierung der erhaltenen 
Reihen mit (8.12) zu berücksichtigen, daß 


Einen Vergleich unserer Ergebnisse mit jenen der neueren Ekmanschen Theorie der 
Meeresströmungen ([2] 1932) können wir dort durchführen, wo Ekman sich ebenfalls auf die 
Berücksichtigung des topographischen Einflusses beschränkt und es ihm gelingt, für den zwei- 
dimensionalen Fall Aussagen über v» herzuleiten. Er nimmt wie wir an, daß die Strömung 


? 
x 
x x 
xo xo 
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über ebenem Boden ankommend an einer Stelle — sagen wir wieder #—=.x, — einen Boden- 
anstieg oder -abfall erreicht und zwar — durch seine Darstellung bedingt — speziell senk- 
recht anströmt. Wenn nun vorausgesetzt wird, daß nur eine schwache Stromablenkung er- 
folgt, so daß eine Linearisierung sowohl in der Störung u*— u, als aueh in »* (siehe unsere 
diesbezüglichen Bemerkungen oben, etwa Schluß von $ 6) wegen |u*— u5|< u), |v*| <u‘, vor- 
genommen werden kann, so gelangt Ekman zu einer bezüglich v linearen Differentialgleichung. 
Diese führt zu dem folgenden Ergebnis ([2], Abschnitt „Vollständigere Theorie“), das wir hier 
nunmehr in unseren Dimensionslosen schreiben: 

x 


xo 
Dies würde mit unserem unter allgemeineren Voraussetzungen gültigen Resultat (8.6a) mit 
v,=0 dann exakt übereinstimmen, wenn man A mit @”' identifizierte. Es gilt in Wirklich- 
keit nun folgende Uinrechnungsformel: 
1 h 2n 
Dsin2a’ 
wo D die für die Bodenströmung gebildete „Ekmansche Reibungstiefe“ 
D=al/ -“, 
00 
bezeichnet, welche etwa der Dicke ö der Bodengrenzschicht entspricht, und wo ferner a den 
Ablenkungswinkel der Strömung am Boden gegenüber der Hauptströmung bezeichnet. Wie 
eine in Kapitel III durchgeführte Abschätzung lehren wird, ist 2-3 jedenfalls im all- 
gemeinen von der Größenordnung O(1) bei Meeresströmungen, so daß wir uns insoweit in 
Übereinstimmung befinden mit diesem Ekmanschen Grenzfall. (Man kann hier einwenden, 
daß auch unser Ergebnis nur für kleine Störungen von ® gilt, da wir bisher die Ver- 
änderlichkeit von h nicht berücksichtigt haben. Dazu verweisen wir jedoch auf $ 11 unten.) 
Die obigen Ergebnisse für Geschwindigkeitsverteilung und Stromlinienverlauf stellen wir 
nun noch in den dimensionsbehafteten Ausgangsgrößen zusammen (die wir vereinbarungsgemäß 


A 


— siehe (7.1) — nun mit einem Stern versehen). Zuvor definieren wir noch, der oben an- 
gegebenen Deutung von £; entsprechend, die Länge {7 (Höhe des „Ersatzberges“) vermöge 
= 
Dann ist wegen G = So Bh 
x” 
& 
x* 
n = 
u 
„ * v, * 
(8.8°): 


x* 


wur 
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x® x* 
20 
. (8.12%). 
% 
Si Bild 3. Ersatzberg |; $ 9. Ein Beispiel. Die Bodener- 
hebung sei beispielsweise durch das 
folgende Profil gegeben: 
für|e]s1j2, (MD, 


0 1/2 


N 
“ -a25 wo £, eine Konstante (voraussetzungs- 
gemäß |£,]<1) bedeutet (siehe Bild 3, 
Kurve G =»). 


Der Ersatzberg £; berechnet sich zu 


0 
25,0 sin?ax 
n? 2 


Er ist in Bild 2 für verschiedene Werte von @ dargestellt und liefert uns eine sehr anschau- 
liche Vorstellung vom Reibungseinfluß. 

Setzt man (9.2) in unsere Formeln ein, so erhält man nach elementarer Rechnung Ge- 
schwindigkeitsverteilung und Stromlinienverlauf. Für die Grundrißprojektion der Stromlinien 
außerhalb der Reibungsschicht ergibt sich insbesondere: 


für 
x +05 
H 
für |«]<1/2; 
für 1/2 
Im reibungslosen Grenzfall ergibt sich 
B 1” 1 2 + (938). 
für > 1/2 


Wir haben diese Formeln für den Fall der senkrechten Berganströmung (v, = v}/u} = 0) 
ausgewertet. Ist »,=F0, so erhält man, solange v,<B/R ist, die Stromlinien hinreichend 
genau, indem man den Stromlinienordinaten des Falles einfach v, (x —x,) additiv 
überlagert. 

Um zunächst den Einfluß des Verhältnisses B/R zu untersuchen, haben wir A/h 


als endlichen festen Wert angenommen und eine Stromlinie, genauer (y — eine 
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reine Funktion von x und @, in Bild 4 für verschiedene Werte von @ aufgetragen. Für den 
Grenzfall G>0 (R]B>0, völliges Überwiegen der Corioliskräfte gegenüber den longi- 
tudinalen Trägheitskräften) erhält man, wie schon am Schluß des $ 7 allgemein diskutiert 
wurde, eine Stromlinie mit Wendepunkt über dem höchsten Punkt des Berges. Dieser Wende- 
punkt wandert nun mit wachsenden Werten von @ (d. h. von R/B) rasch stromabwärts, um 
bald ganz außerhalb des Interessenbereiches der Bergumgebung zu geraten. (Es muß beachtet 
werden, daß unser Bild nur die Abhängigkeit von der Bergbreite bei irgendeiner festen An- 
strömgeschwindigkeit (festem R) wiedergibt, da ja auch h streng genommen und jedenfalls 
in einem größeren Variabilitätsbereich der Geschwindigkeit als geschwindigkeitsabhängig an- 
gesehen werden muß.) € 


Bild 5 (rechts). Stromlinien in Abhängigkeit von h/A. | 


Bild 4 (unten). Stromlinien in Abhängigkeit von B/R. -05 


3 


- 
! | 


—ı 0 05 5 


Der Einfluß des Verhältnisses A/H, also, was uns hier in erster Linie interessiert, 
der Einfluß der Sekundärströmung in der bodennahen Reibungsschicht, wird in Bild 5 wieder- 


—1 
gegeben. Hier ist R/B als fester endlicher Wert gedacht und demgemäß (y — Y,) (% SR) 


B 
R ’ 
aufgetragen. Es ist also, wenn wir H, B und R einzeln fest denken (also auch die Anström- 
geschwindigkeit «;) der Einfluß verschieden starker (laminarer oder turbulenter) Zähigkeit der 
bodennahen Schicht wiedergegeben. Man sieht, daß die für =» (h/HT—=0, reibungsloser 
Grenzfall) sich ergebende Kurve bereits durch die Kurven für endliche, nicht sehr große @ 
gut approximiert wird. Für @= 10 liefert die Reibung schon keinen wesentlichen {Beitrag 
mehr und bereits bei G=3 ist die Abweichung von G=x praktisch bedeutungslos. (Das 
entspricht auch den Verhältnissen des Bildes 3 gemäß unserer Deutung von £;.) Für solche 
Werte von @ wird man also bedenkenlos unter Vernachlässigung der Reibung rechnen dürfen. 

Die Versetzung Ay eines Flüssigkeitsteilchens in y-Richtung an der Stelle x gegenüber 
der verlängerten Anströmlinie ist im Falle »,=0 allgemein durch 


wieder eine reine Funktion von x und @ wegen Ya G: für verschiedene Werte von @ 


x 4 
% &o 
gegeben (und nach einer Bemerkung im Anschluß an (9.3) also auch für niehtverschwindende »,, 


solange |v,)<B/R, worauf wir uns der Einfachheit halber beschränken). Für hinreichend 
große @ ist dafür genähert 


Xo 


zu schreiben. Der gemachte Fehler läßt sich leicht unter Heranziehung von (8.5) abschätzen. 
Im behandelten Beispiel lehrt das Bild 3, daß man für e—x,=1, also nach Zurücklegung 
des Weges B, schon für den noch verhältnismäßig kleinen Wert @=3 einen prozentualen 
Fehler von weniger als 8°/o macht, wenn man (9.4) durch (9.5) ersetzt, für G=10 einen 
solchen von 3°Jo. 


| | 02 
o 04 
| 
| | 
| 3 | 
| 
-05 o 
. 
XUM 
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In unserem Beispiel ist nach der Näherungsformel (9.5) für hinreichend große @ 


für >1/2 


Die mit »,—0 ankommende Strömung ist nach Überschreiten des Berges um einen gewissen 
Winkel y abgelenkt worden, der durch 


gegeben ist. = 

Es sei etwa B=200km, 5 =10"'H, =17,292.10”° sec'. Ist dann 
= 10 cm sec”! (20 cm sec”'), so ist R= 1,37 km (2,74 km), y= — 8212’ (— 74°40’), 
= —7,29(— 3,65) für so daß unmittelbar hinter dem Berge («= 1/2) also 
Ay*=—3,65 B(— 1,82 B) ist. Die zum Überschreiten des Berges benötigte Zeit t beträgt 22 Tage 
(11 Tage). Die gewählten Geschwindigkeiten entsprechen dabei jener Größenordnung, wie sie 
im Meere vorkommt. Wählen wir im Hinblick auf die Anwendungen bei atmosphärischen 
Strömungen etwa «,—20 m sec”', was einem frischen Wind in Bodennähe entspricht, so ist 


R=274km, y=- 95, _ _.0,0365 (für «>1/2), Ay*—= —00182B für = 1/2, 


t= 2,64 Stunden. (Dies gilt alles natürlich nur unter der Voraussetzung, daß @ hinreichend 
groß ist.) 


$ 10. Isobaren und freie Oberfläche. Für den Druck p* in der Hauptströmung schreiben 
wir mit Rücksicht auf die dritte Eulergleichung (3.2) (statische Grundgleichung) und auf (3.6) 
p; ist dann der durch die Bodenerhebung hinzukommende Störungsanteil des Druckes. 
Setzt man in den beiden ersten Eulerschen Gln. (3.2) die berechneten Geschwindigkeits- 
komponenten u, v gemäß (8.6), (8.6a) ein, so ergibt sich mit p, =p}/lo ws? in dimensionsloser 
Schreibweise 


dd B 
somit unter Verwendung von (8.5) - 
x 4 
2, ()=— \ 


Es ist also 
x 4 
.... 


In einer bestimmten Höhe 2=const sind die Isobaren demnach durch 


H B 


gegeben. Ein Vergleich mit den Stromlinien (siehe $ 8) 


x x A 
x x 


ergibt eine Abweichung Ay = Yıson. — Yatr, VON 
x 


xo 


Z 


B o B 
| 
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Nehmen wir etwa den Fall v»,=0 der senk- 
rechten Berganströmung, so beträgt die auf die 
Stromlinienabszisse bezogene prozentuale Abwei- 
chung 


Für R/B<<1 (Meeresströmungen über ausge- 
dehnte Berge, etwa R/B= 10°?) ist diese vernach- 
lässigbar klein (es sei denn, daß durch besondere 
Verhältnisse das Doppelintegral im Nenner an der 
betreffenden Stelle x klein vom Betrage gegenüber 
I&;| wird; bei einem normalen Berg und großen 
Werten von @ (£;>{) kommt dies nicht vor). Für 
Werte von R, die mit B vergleichbar sind (atmo- 
sphärische Strömungen mit ihren wesentlich größe- 
ren Geschwindigkeiten), kann dagegen die Isobare 
stark von der Stromlinie abweichen. Für - 
und R/B=1/2 sind für das Beispiel des $ 9 Strom- en 
linien- und Isobarenfeld in Bild 6 eingezeichnet. Bild 6. Stromlinien und Isobaren einer Strömung 
Es ist in diesem Fall für 1/2 mit u=10 msee-!, BIA=10? undG=«. 


1 1 
(y— Yılzır. =7 +1) — + 5) 


1 1 
1) — |, 


und in beiden Fällen y — für 1/2. 

Starke Abweichungen können aber auch bei allen Geschwindigkeiten entstehen, wenn 
der Berg unter flachem Winkel zum Kamm angeströmt wird (Riehtungszahl v,—= »}/u) der 
Anströmung von großem Betrag). 

An der freien Oberfläche ist p—p,=0. Verstehen wir unter H, die ungestörte Ober- 
fläche der reinen Parallelströmung ®, (©=0), so ist nach (10.2) mit £;=0 


Dies ist der zu ®, gehörige ebene, senkrecht zur Strömungsriehtung geneigte Spiegel. 
Bei Vorhandensein einer Bodenerhebung stellt sich die freie Oberfläche H gemäß 


x x 
ein. Für die praktische Rechnung wird wiederum folgende Umformung von Nutzen sein 
x 
Die von uns entwickelte Theorie machte wesentlich Gebrauch von der Voraussetzung, 
daß H—1 (d. h. H*— H) hinreichend klein sei, um für alle beabsichtigten Anwendungen bei 
ozeanischen und atmosphärischen Strömungen in dem Ausdruck (4.6) für die Vertikalgeschwin- 
digkeit w* die freie Oberfläche H* durch H zu ersetzen. Nur unter dieser Voraussetzung 
bestehen also auch die eben abgeleiteten Ausdrücke für die Spiegeleinsenkung zu Recht. Wir 
werden auf diese Frage erst dann eingehen, wenn wir eine Untersuchung der im Meer und 
in der Atmosphäre auftretenden Größenordnung der charakteristischen Parameter durchgeführt 
haben (siehe $ 14). 


$ 11. Berücksichtigung der Veränderlichkeit von h. Wie bereits bei Einführung der 
Größe h in $ 4 betont wurde, zeigt diese für den Reibungseinfluß maßgebende Länge einen 
merklichen Gang mit der Geschwindigkeit ® der Hauptströmung. Sie hängt im übrigen natür- 
lich auch von der Beschaffenheit der turbulenzerzeugenden Bodenrauhigkeiten ab und von der 
Dicke der Reibungsschicht. Hierüber erfolgen in Kapitel III nähere Angaben. 


| 
= 
| 
” 
| 
.| x 
| 
_&B 
Yolısob. 
XUM 


. Z.a . Math. Mech. 
296 Görtler, Strömungen über der rotierenden Erde Baur Nr.5 Okt. 1941 


Wenn wir bisher in diesem Kapitel zur Vereinfachung der Rechnung die Länge h durch 
einen konstanten mittleren Wert in den Ausgangsgleichungen ersetzt dachten, so kann man 
vom Rechenergebnis nur eine erste Näherung für nicht zu starke Veränderungen von h im 
Strömungsbereich erwarten. Die folgende Rechnung wird zeigen, daß dies zutrifft, und zwar 
in dem Sinne sogar, daß für die Gültigkeit der bisherigen Näherungsrechnung wohl der 
Schwankung von h eine obere Schranke gesetzt ist, daß aber h innerhalb dieser Schranken 
durchaus schnellen Veränderungen unterworfen sein darf. Solche plötzlichen Änderungen 
von h über kurze Strecken kommen in der Natur etwa vor, wenn sich die Bodenbeschaffenheit 
plötzlich ändert. 

Bei im einzelnen bekanntem Strömungsverlauf ist h eine bekannte Funktion der Orts- 
koordinaten, im vorliegenden zweidimensionalen Fall insbesondere h=h(x). Wir wollen 
nun im folgenden h zunächst von vornherein als gegebene Funktion von x ansehen. Wir 
umgehen damit die Schwierigkeit, eine komplizierte Geschwindigkeitsabhängigkeit der Länge 
h auf Grund einer — bis heute allein möglichen — groben Abschätzung in unsere Ausgangs- 
gleichungen einführen zu müssen. 

Wir beschränken uns im folgenden zur Vereinfachung der Formeln auf den Fall der 
senkrechten Anströmung unserer Bodenerhebung: v,=0. Unsere dimensionslosen GIn. (7.4) 
und (7.5) sind dann mit h=h(x) zu ersetzen durch 


d 


Dabei ist h,—=h(x,) und @, ist die mit h, gebildete für den Reibungseinfluß charakte- 
ristische Dimensionslose @, also 
Die Gl. (7.6) für w bleibt natürlich formal ungeändert. 
Die Behandlung der Differentialgleichung (11.1) kann allgemein mit denselben Mitteln 
durchgeführt werden, wie sie bisher angewandt wurden, also durch den Linearisierungs- 
ansatz (8.1): «=1+[;, und zwar wiederum ohne neue Annahme allein als Konsequenz der 


unserer Theorie bereits zugrunde liegenden Voraussetzungen des Kapitels I. An Stelle von 
Gl. (8.2) tritt dann 


(113). 


mit G= z nn G (x). Durch Integration dieser Gleichung von x, bis zu einer variablen 
Stelle « des Strömungsbereichs erhält man an Stelle von (8.5) die Beziehung 


Benutzt man diese Beziehung zur Berechnung von v aus (11.2), so erhält man insgesamt 
wieder formal dieselben Relationen (8.6) für die Geschwindigkeitsverteilung (das würde auch 
noch gelten, wenn v,==0), nämlich 


x 


B 

w(x, 2) - Zu 2) - 


Hieraus folgt, daß auch allgemein der dimensionslosen Störungsgröße £; die ihr früher ge- 
gebene Deutung als Kontur des „Ersatzberges“ zukommt. Insoweit hat sich nichts geändert. 
Damit bleiben auch die Ausdrücke (8.7), (8.8) und (8.9) für Stromlinienverlauf und Strom- 
ablenkungswinkel bestehen, nur daß eben das frühere £; durch einen verbesserten, die Ver- 
änderlichkeit von h berücksichtigenden Ausdruck zu ersetzen ist. Dieser an Stelle von (8.4) 
tretende Ausdruck, das Integral der Differentialgleichung (11.4) lautet: 
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x \ 


Xo 
(Für die praktisch vorkommenden großen Werte von @ kann man nach den schon früher an- 
gegebenen Methoden geeignete Reihenentwicklungen herstellen.) 

Aus diesem Ausdruck (11.7) erkennt man unmittelbar, daß der Ersatz von h (x) (und 
damit von @(x)) durch einen geeigneten Mittelwert als erste Näherung für die Berechnung 
der Stromablenkung berechtigt ist, wenn Ah (x) von geringer Schwankung ist. Dabei Bm u 


keine Rolle, ob in diesen Grenzen plötzliche Schwankungen, d. h. große Werte von u; 


vorkommen, da nur h selbst und Integrationen über diesen Ausdruck maßgebend sind. 
Starke Veränderungen von h können in der Natur auftreten, etwa wenn sich ® im 


Strömungsbereich stark ändert (nach (11.6) also für große d ‚ also z.B. bei Meeresströmungen 


mit ihren kleinen Geschwindigkeiten) oder bei starker Änderung der Bodenrauhigkeit (etwa 
wenn Luft vom Meere herein über stark zerklüftete Gebirgslandschaft hinwegströmt, um ein 
krasses Beispiel zu nennen). Wie man sich hier helfen kann, wenn es notwendig werden sollte, 
sei kurz skizziert. Man zerlege den Bereich B in eine genügende Anzahl von Teilinter- 
in denen man einzeln mit der ge- 
forderten Genauigkeit mit konstantem G@=@G, rechnen darf. Diese Rechnungen werden 
Schritt für Schritt in Strömungsriehtung durchgeführt, wobei für jedes @, zunächst ein im 
voraus zu schätzender Mittelwert in B, eingesetzt wird. (Dazu können die Ausführungen 
des Kapitels III dienen) Wenn sich dieser Wert auf Grund des berechneten ® als zu 
schlecht erweist, kann die Rechnung mit einem nun leicht besser zu schätzenden @, wieder- 
holt werden. Ist die Strömung in B, berechnet, so kennt man die Geschwindigkeiten in 
«=xr, und damit die Anfangsbedingungen für das Intervall B,;,. (Es sei bemerkt daß 
unsere in diesem Kapitel hergeleiteten Formeln nicht unmittelbar hier zu verwenden sind, 
da wir an der Eintrittsstelle des Strömungsbereichs spezielle Aufangsbedingungen angesetzt 
haben. Es wären also zunächst die Ausdrücke bei Anfangsbedingungen allgemeiner Art ein 
für allemal herzustellen, was keiner neuen Schwierigkeit begegnet.) 

Wie schon angekündigt, wird eine Abschätzung lehren, daß sich die Werte von @ bei 
atmosphärischen und ozeanischen Strömungen im allgemeinen in Grenzen bewegen, für die 
nach unserer Theorie der Reibungseinfluß keine bedeutende Wirkung ausübt. Insofern werden 
die eben angestellten Betrachtungen hier wohl kaum eine praktische Rolle spielen. Wesent- 
lich bleibt allein, daß mit unserer Theorie der Reibungseinfluß allgemein erfaßt werden 
konnte. 


III. Anwendung auf Strömungen im Meer und in der Atmosphäre. 


In diesem Kapitel wird zunächst die Frage der Größenordnung der auftretenden Dimen- 
sionslosen, insbesondere der für den Einfluß der Sekundärströmung maßgebenden Größe @ 
bei den’im Meer und in der Atmosphäre vorkommenden Strömungsverhältnissen erörtert. 
Ferner ist zu prüfen, ob die Vernachlässigung der Spiegeleinsenkung in der Gl. (4.6) bei 
diesen Verhältnissen berechtigt war. 


$ 12. Approximation der Geschwindigkeitsverteilung in der bodennahen Reibungsschicht. 
Eine strenge Untersuchung des Einflusses der Sekundärströmung in der bodennahen turbu- 
lenten Reibungsschicht wäre nur möglich, wenn die Gesetzmäßigkeiten beschleunigter turbu- 
lenter Strömungen genau erforscht wären. Bei dem heutigen Stand der Forschung bleibt 
uns nur der Weg einer rohen Abschätzung übrig. 

Die Grundlage dieser Abschätzung geht auf Ausführungen von L. Prandt] zurück. 
Ich entnehme sie der demnächst erscheinenden Neuauflage des Prandtlschen Lehrbuches 
„Abriß der Strömungslehre“. Herr Prof. Prandtl hat mir hierzu das Manuskript des be- 
treffenden Abschnitts: „Reibungswind und ähnliches“ zur Verfügung gestellt. 

Wir betrachten mit Prandt| folgenden verhältnismäßig einfachen Fall. Außerhalb der 
turbulenten Reibungsschicht über einem rauhen, im Mittel ebenen, rotierenden Boden herrsche 
eine Parallelströmung mit der Geschwindigkeit 8*={U*,0}. Die durch ®* erzeugten Coriolis- 
kräfte müssen dann mit einem entsprechenden Druckgefälle im Gleichgewicht stehen. Der 
Druck ist nach (3.6) durch p* =pt —20w’ U* y* gegeben. 

Wir approximieren die mittlere Geschwindigkeitsverteilung in der turbulenten Grenz- 
schicht dieser Strömung durch folgenden Ansatz. Zunächst wird «* durch das gewohnte 
Potenzprofil 

(5) 
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genähert. Die Länge ö (ein Maß für die Grenzschichtdieke) und der Exponent 1/n sind noch 
geeignet zu wählen. Für die Geschwindigkeit v* der durch die Erddrehung hervorgerufenen 
Sekundärströmung senkrecht zu ®* wird angesetzt: 


g* 


Damit ist »*—=0 für *=dÖ und lim tg a. Es ist also « der Winkel, den die Strömung 
unmittelbar am Boden mit der Hauptströmung bildet. 
Da nun die Komponenten r,, r, der Schubspannung am Boden z*=0 sich aus den 


d 
Corioliskräften wegen = —20'0v*, = +2o'ou* zu 
02 02 dy 
2n’tga 
(12.3) 
\u* — u*) d = ow’ U*Ö 
1- 
. berechnen, wird tga wegen tga= = durch n vermöge 
x 


ausgedrückt. 
Es gilt nun, d und » geeignet zu wählen. Dazu bedient sich Prandtl des bekannten 


logarithmischen Gesetzes für die Geschwindigkeitsverteilung einer turbulenten Grenzschicht- 
strömung über ruhenden rauhen Flächen: 


Darin ist k die Höhe der Bodenrauhigkeiten, während die Zahl € von der jeweiligen geome- 
trischen Form und von der Verteilung der Rauhigkeiten abhängt. Wir verweisen auf die 
Messungen von H. Schlichting [7] an verschieden gearteten Rauhigkeiten. Dort wird 
statt k in (12.5) eine reduzierte Rauhigkeitslänge eingesetzt, die äquivalente Sandrauhigkeits- 
höhe k,, mit dem Erfolg natürlich, daß die Zahl, die dann an Stelle von © zu treten hat, 
nämlich C +5,75logk,/k, eine universelle Konstante (= 8,48) ist. Ferner sei auf die über 
verschieden bebautem Gelände ausgeführten Messungen von W. Paeschke [3] verwiesen. 
(Auch hier wird statt k eine anders definierte Länge, eine „wirksame Rauhigkeitshöhe“ in 
(12.5) eingesetzt. Sie stimmte bei den dort betrachteten Rauhigkeiten im allgemeinen mit 
der wirklichen Rauhigkeitshöhe überein und es ergab sich an Stelle von C die Zahl 4,98. 
Diese „wirksame“ Rauhigkeitshöhe kann aber natürlich in ausgearteten Fällen stark ‚von der 
wirklichen Höhe der Rauhigkeiten abweichen. So hat ein Weizenfeld (Höhe 130 cm) nach 
Paeschke dieselbe wirksame Rauhigkeitshöhe wie hohes Grasland (30 em), was man ja 
auch rein anschaulich zu erwarten hat.) 

Bei den hier interessierenden nicht ausgearteten Rauhigkeiten wird C in (12.5) zwischen 5 
und 8,5 variieren. Für unsere Abschätzung müssen wir C festlegen und wählen aus Bequem- 
lichkeitsgründen später C=5,75. (Solange k hinreichend klein gegenüber ö ist, spielt ein 
kleiner Fehler in der Wahl von € keine Rolle. Für eine genauere Betrachtung ist bei dieser 
festen Wahl von € unter k nicht die wirkliche, sondern eben eine auf diesen (-Wert reduzierte 
Höhe zu verstehen, die sich aus log k/k»=1 — (,/5,75 berechnet, wo nun k, die wirkliche 
Rauhigkeitshöhe und C,, die aus Messungen zu ermittelnde zugehörige Zahl des Gesetzes (12.5) 
bezeichnen.) 

Der Bereich der hier zu berücksichtigenden Rauhigkeitshöhen erstreckt sich von einigen 
Zentimetern (Schneefläche etwa) bis zu Größenordnungen von 100 m (zerklüftetes Gebirge). 
Die geeignete Wahl von k für Meeresströmungen dürfte nicht immer leicht fallen. Durch 
Dichteschichtungen etwa können überdies die jeweiligen Rauhigkeiten stark abgemildert in 
Erscheinung treten. 

Das Gesetz (12.5) soll nun in zweifacher Hinsicht verwertet werden. Einmal soll daraus 
das Maß ö der Grenzschichtdieke ermittelt werden. Setzt man in (12,5) 2*=6ö und also 
U* 


ü 
fe 
Ü 
u 
G 
d 
ni 
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und vergleicht man mit (12.3) unter Berücksichtigung von (12.4), so erhält man als Be- 
stimmungsgleichung für 6: 
Es ist nun noch n geeignet zu wählen, d.h. so, daß (12.1) im Mittel gut mit (12.5) zur 
Deckung kommt. Eine geeignete Formulierung dieser Forderung besteht darin, daß die Ver- 


drängungsdicke | — d. i. diejenige Höhe, um welche die äußere Parallel- 


v 
strömung durch Ausbildung der Grenzschicht vom Boden abgehoben wird — in beiden Fällen 
den gleichen Wert haben soll. Dies ergibt 


also n als Funktion von =. Damit ist auch tga als Funktion von 1 nach (12.4) bekannt. 


$ 13. Größenordnung der charakteristischen Dimensionslosen im Bereich der ozeanischen 
und der atmosphärischen’ Strömungen. Die im vorangehenden Paragraphen nach L. Prandtl 
durchgeführte Abschätzung gilt zunächst nur für den vorausgesetzten einfachen Fall einer 
Parallelströmung außerhalb der turbulenten Reibungsschicht über einen im Mittel ebenen 
rauhen Boden. Die Anwendung auf die in unserer Theorie behandelten Strömungen mit 
nicht geradlinigen Isobaren ist aber zulässig, soweit wir uns dessen bewußt bleiben, daß 
die herangezogene Abschätzung nur größenordnungsmäßig in Anspruch zu nehmen ist. 

Die gemäß (4.3) in unsere Theorie eingeführte Länge h berechnet sich wegen 


V2n+1 
(n + Din +2) 


zu 


Zur übersicht- 
lichen Vergegenwär- 
tigung der durch 
unsere Formeln grö- 
ßenordnungsmä- 
ßig gegebenen Zu- 
sammenhänge ist das 
Nomogramm Bild 7 
entworfen worden. 
Es wurde dabei aus 

Bequemlichkeits- 
gründen C=5,75 
gewählt, da es im 
folgenden auf den 
Wert von © nicht 
wesentlich ankommt. 
Ferner wurde im 
Hinblick auf unsere 
Anwendungen = 
7,29.10°° sec”! ge- 
setzt. 

Zunächst ist 
über ö/k (Skala links 
und rechts außen) 
zu verschiedenen 
festen Werten von 4 
U* im Bereich der zum 
uns interessierenden 
Geschwindigkeiten 
der zugehörige Verlauf von ösingp nach (12.7) unter Benutzung von (12.8) als Kurvenschar 
U* — const aufgetragen. Ferner geben die im doppeltlogarithmischen Netz des Bildes schräg 
nach rechts ansteigenden Geraden die Linien konstanter Rauhigkeitshöhen, genauer k sin o 


— Ösiny [m] 


Bild 7. Nomogramm zur groben Abschätzung von U*,k,d,G,h,n,«. 
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—=onst, wieder. Man kann damit in einfacher Weise zu vorgegebenen Werten der Rauhig- 
keitshöhe k, der Geschwindigkeit U* und der geographischen Breite @ den zugehörigen 
Wert von ö entnehmen. 

Auf rechtsstehenden Skalen sind die Werte von h/ö nach (13.1), » nach (12.8) und a 
nach (12.4) zu den in gleicher Höhe liegenden Werten von ö/k angegeben. 

Uns interessiert hier in erster Linie die in Frage kommende Größenordnung der 


IT 


Dimensionslosen G=5 PR Es ist daher noch in einer weiteren rechtsstehenden Skala die 


Größe @G- = 50 nach (12.7) in Verbindung mit (13.1) eine reine Funktion von ö/k: 
B 0) 
aufgetragen. 


Als Nomogramm-Schlüssel ist gestrichelt folgendes Beispiel einer Ablesung eingetragen: 
Zu m/see und ksing=05 m (etwa: k=1m bei 9=30°) liest man ab: ösing 


n=84; h 0,030; G—27-10° (also etwa für 100 7 wird 


G= 27). 

Was die Verhältnisse bei Meeresströmungen anbelangt, so ist U*—=0,5m/sec bereits 
als große Geschwindigkeit zu bezeichnen. In der Atmosphäre entspricht 20 m/see der 
Windstärke 10 nach der Beaufortskala. In Bodennähe würde dem ein frischer Wind ent- 
sprechen. Unser Nomogramm umfaßt bis U*=50 m/sec alle interessierenden Windstärken. 

Als Rauhigkeitshöhe k am Boden der Grenzschicht kommt ihrem Begriffe nach nur eine 
Größe in Frage, die klein ist im Vergleich zur Dicke ö der turbulenten Grenzschicht. Diese 
Rauhigkeiten erzeugen die Turbulenz. (Bodenerhebungen, die groß gegenüber diesen Rauhig- 
keitshöhen sind, geben gemäß den Formeln unserer Theorie die Stromablenkung.). Unser 
Nomogramm ist also nach unten begrenzt, und zwar haben wir bei dem bereits zu kleinen 
Wert ö/k=5 abgebrochen. 

H 


Man kann, wie aus dem Nomogramm ersichtlich ist, annehmen, daß @ >3- 10? B ist, 


im allgemeinen sogar wesentlich größer. Nimmt man beispielsweise za—-r an (bei einer 


Meerestiefe 7=5 km also B=500 km, das sind fast 5 Äquatorgrade; bei der Höhe 7=8 km 
der gleichförmigen Atmospäre B=800 km), so nimmt G@ Werte zwischen 3 und 80 an, im 
allgemeinen also Werte, für die nach $ 9 der Sekundärströmung keine Bedeutung für den 
Verlauf der Hauptströmung mehr beigemessen werden kann. Gewiß ist die Angabe über die 
untere Grenze von @ wie alle hier gemachten Angaben nur größenordnungsmäßig richtig, wir 
müssen also die Möglichkeit des Eintretens von Grenzfällen zulassen, für welche @ hinreichend 
klein ausfällt, um einen merklichen oder gar wesentlichen Einfluß der Reibung in der boden- 
nahen Grenzschicht in Erscheinung treten zu lassen. Nur jedoch, wenn man B wesentlich 
größer als in dem eben angeführten Beispiel annimmt, kann G<1 werden. Bei dem sehr 
großen Wert B=2000 km kann nach unserem Nomogramm @ für Meeresströmungen bei 
entsprechend großen Rauhigkeiten einen Wert unter 1, etwa 0,75 annehmen. Für derart große 
“o' Werte von B aber wollen wir unsere 
| Theorie, welche die Erdoberfläche durch 
% | ihre Tangentialebene ersetzt, nicht in 
| .* Anspruch nehmen. Man bedenke dabei, 

daß die Stromlinienablenkung beim Über- 

schreiten der Bodenerhebung ein be- 

. trächtliches Vielfaches von B erreichen 
10° kann (gerade bei Meeresströmungen, 
Fi siehe das Beispiel am Schluß des $ 9). 

Immerhin darf jedoch nicht vergessen 
m ee werden, daß unsere Abschätzung streng 
ri genommen nur dort angewandt werden 

en darf, wo die B-Strömung weiträumig ge- 
Br nug ist, um an jeder Stelle durch eine 
geradlinig-gleichförmige Strömung für obi- 
Th ge Abschätzung ersetzt werden zu können. 
10° Die Länge h zeigt nach der Ab- 


| 
| 
| 


10 — 


025 


Bild 8. Länge k in Abhängigkeit von der Geschwindigkeit U*. schätzungsformel (13.1) eine recht merk- 


| 
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liche Abhängigkeit von der Geschwindigkeit. Für den ganzen im Nomogramm Bild 7 umfaßten 
Bereich ist diese in dimensionsloser Auftragung in Bild 8 wiedergegeben. Nähert man die 
dargestellte Kurve in der doppeltlogarithmischen Auftragung durch ihre Tangente im Kurven- 
punkt h,/k, U%|o’k, so ergibt sich ein gebietsweise gültiges Näherungsgesetz der Form 


wo im dargestellten Bereich 0,45 <A < 0,85. 


In $ 11 wurde ausgeführt, wie die Veränderlichkeit von h in unserer Theorie berück- 
sichtigt werden kann. In entsprechender Weise kann man eine Korrektur für die Veränderlich- 
keit der Verdrängungsdicke der Reibungsschicht anbringen — die Kontur £ in unserer Theorie 
ist der um die Verdrängungsdicke gehobene Boden, siehe $ 4 —, indem man einen ersten 
gröberen Rechnungsgang mit konstanter Verdrängungsdieke durchführt, aus dem erhaltenen 
Geschwindigkeitsfeld und unserer Abschätzungsformel die Veränderungen der Verdrängungs- 
dicke abschätzt und die Rechnung mit einem so verbesserten { nochmals durchführt. Da die 


Verdrängungsdicke durch u nach (12.1) gegeben ist, wird sie bei Meeresströmungen, wo 


eine solche Korrektur eventuell erwünscht sein könnte, jedoch im allgemeinen klein gegen- 
über Sein. 


Übrigens sei noch erwähnt, daß auch plötzliche Änderungen unseres Parameters @ durch 
sprunghaftes Ändern der Rauhigkeitshöhe k möglich sind, etwa wenn ein Wind vom Meer 
hereinströmend die Küste erreicht, oder vom Boden mit niedrigem Bewuchs her über einen 
Wald, eine Großstadt oder dgl. weht. 


$ 14. Größenordnung der Spiegeleinsenkung. Nach $ 10 ist die durch die Bodenerhebung 
bewirkte Spiegeleinsenkung gegeben durch 


B? 
&Xo 


und es ist eine noch zu bestätigende Voraussetzung unserer Theorie, daß diese bei ozeanischen 
und atmosphärischen Strömungen dem Betrage nach klein bleibt gegenüber |{|. Wie wir 
wissen, ist £;=0({) und wegen $ 13 im allgemeinen sogar [;=/L. 


Zunächst zum Falle der Meeresströmungen. Der Faktor uj?’/g H ist, wenn wir H=3 km 
annehmen und die für ozeanische Strömungen sehr große Geschwindigkeit u; = 60 em/seec zu- 
“2 2 2 
grunde legen, mit g— 981 cm/sec? von der Größenordnung Der Faktor 
ist bei 9=30° und dem für unsere Theorie großen Wert B=500 km etwa 5-10°°? (bei 
B=300 km etwa 1,7-10°?), und das Doppelintegral über £,; sorgt mit dafür, daß auch der 
zweite Summand in (10.6) von der Größenordnung 10°? 0(£) wird. Insgesamt ist also 


Im Falle der atmosphärischen Strömungen ist H=8 km anzusetzen und u, — 20 m/see 
entspricht einem frischen Wind am Boden. Hierfür ist «5?/g H=5-10°®. Ferner ist wieder 
4 ‚2 B? 
mit B=500 km und der Faktor Also auch hier ist 
Jedoch ist überhaupt die Fragestellung hier von geringer Bedeutung, da ja unser Modell der 
Atmosphäre gerade hierin der Wirklichkeit gar nicht entspricht. 


Neben der Spiegeleinsenkung durch den Einfluß der Bodenerhebung ist die Spiegel- 
neigung der ungestörten Oberfläche bei der Parallelströmung 3}; = {uw}, v}}, nach (10.5) durch 


gegeben, abzuschätzen. Denn auch diese muß klein bleiben, da wir H durch 1 ersetzt haben. 
Mit x*, y"= 0 (500 km) (also ein Strömungsbereich von einer waagerechten Erstreckung bis zu 
etwa 1000 km, da &’,y’ vom Zentrum O0 des Bereiches aus gemessen werden) ist im Falle 
der Meeresströmungen — nehmen wir = 0(30 cm/sec), H=5km, 9—=30° — an den 
Begrenzungen des Bereichs |H,—1|=0(10°*), also gewiß im Strömungsbereich als klein 
gegenüber |{]Max (das etwa als 0 (10°) angesetzt werden kann, da wir quadratische Glieder 


| 


302 Görtler, Strömungen über der rotierenden Erde ei er 


in ögegen 1 vernachlässigt haben). Im Falle der atmosphärischen Strömungen wird |H,—1| 
mit 5,05 = 0 (20 m/sec), H=8km, 9=30° etwa 0(10°°). Allein für atmosphärische Strö- 
mungen ist die ganze Fragestellung, wie soeben schon bemerkt, sowieso von geringer Bedeu- 
tung: Wenn wir schon die gleichförmige Atmosphäre als Ersatz der wahren Verhältnisse an- 
nehmen, könnten wir auch von vornherein eine homogene Luftschicht mit horizontalem Spiegel 
(etwa von einem ebenen, natürlich reibungslos zu behandelnden Deckel abgegrenzt) zugrunde 
legen, also H— H vorschreiben und uns damit dieser selbstverschuldeten Schwierigkeit entheben. 


IV. Schlußbemerkungen. 


$ 15. Ausblick. Die in Kapitel II durchgeführte Rechnung für das zweidimensionale 
Problem erfuhr von vorneherein eine wesentliche Vereinfachung dadurch, daß die Grund- 
gleichungen (5.1) und (5.2) sich ohne weiteres je einmal integrieren ließen. Es blieb nur die 
Behandlung einer einzigen nichtlinearen Differentialgleichung erster Ordnung übrig. Im all- 
gemeineren Fall räumlich beliebig ausgedehnter Bodenerhebungen z={(x,y) sind die voll- 
ständigen Gleichungen (5.1) und (5.2) zugrunde zu legen, aber es tritt eine ganz entsprechende 
Erleichterung ein, da man sie in der Gestalt (dimensionslos) 


Div (Rot Bist kalar) . . .... (51) 


schreiben kann. Dabei sind die Geschwindigkeiten ® und B={— »,u} mit irgendeiner für 
den Strömungsvorgang charakteristischen festen Geschwindigkeit V, dimensionslos gemacht und 
R=V,[2®' ist der zu dieser Geschwindigkeit gehörige Trägheitskreisradius. Hiervon ab- 
gesehen sind die Dimensionslosen wie in unserer obigen Theorie gebildet. Die Behandlung 
des Problems allgemeiner Bodengestalt erscheint demnach recht aussichtsreich, und sie wird 
gegenüber der Idealisierung des von uns behandelten „unendlich langen“ Bergkamms wert- 
volle neue Aufschlüsse bringen. Das Problem erfährt zur Zeit seine Bearbeitung. 


In unserer Theorie ist, wie einleitend ausgeführt wurde, eine bewußte Beschränkung auf 
den topographischen Einfluß vorgenommen, um zunächst dessen Wirkung allein auf den 
Strömungsverlauf zu untersuchen. Eine wirkliche Abtrennung dieses Einflusses von sämt- 
lichen anderen, etwa auch von dem planetarischen Einfluß, ist jedoch, solange man die Reibung 
berücksichtigt, nicht möglich (siehe hierüber auch V. W. Ekman [2], S. 435 ff.), d. h. das 
Zusammenwirken beider Einflüsse ist nicht durch Addition der Einzelwirkungen zu erhalten. 
Berücksichtigt man in den Eulerschen Gin. (3.2) die Veränderlichkeit der Breite o, so erhält 
man neben der unveränderten Kontinuitätsgleichung (15.2) eine um ein Zusatzglied erweiterte 
Gleichung an Stelle von (15.1), in welchem der planetarische Einfluß zum Ausdruck kommt. 
An der Gestalt der Gl. (15.1) ändert sich jedoch formal nichts, wenn man nur jetzt berück- 
sichtigt, da& R=V,/2» sin als ortsveränderliche Größe zu betrachten ist. Auch die Unter- 
suchung dieses allgemeineren Problems des Zusammenspiels von topographischen und plane- 
tarischen Einflüssen wird zur Zeit durchgeführt. 


Besonders einfach gestaltet sich selbstverständlich die Theorie des reibungslosen 
Falles (k =0), der nach dem Ergebnis unserer obigen Abschätzung des Reibungseinflusses be- 
sonderes Interesse verdient. Hier empfiehlt es sich, neben der Kontinuitätsgleichung (15.2) 
an Stelle der nichtlinearen Gl. (15.1) den folgenden Wirbelsatz zu verwenden, den wir zum 
Schlusse kurz ableiten wollen. 


Da das zum Felde der Hauptströmungsgeschwindigkeiten parallele Vektorfeld B(1— JA 
divergenzfrei ist, kann man durch Einsetzen von B=A(1— £)"' in (15.1) sofort schließen, daß 


ist, also gilt längs den einzelnen Stromlinien der Hauptströmung ® der Wirbelsätz: 


oder dimensionsbehaftet: 


Rot®=— 
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(C, C* auf Stromlinien der V-Strömung konstant). Aus dieser Gleichung erkennt man un- 
mittelbar das Zusammenspiel der topographischen und der planetarischen Wirkungen. 
Sofern man bewußt von dem planetarischen Einfluß absieht und also @ = const setzt, er- 
gibt sich bei Annahme einer drehungsfreien Anströmung vor dem Berg einfach 
B 


Rott 


V. Zusammenfassung. 


Es wird eine allgemeine Theorie der stationären Strömung einer Flüssigkeitsschicht 
konstanter Dichte über weitausgedehnte schwache Bodenerhebungen auf der rotierenden 
Erde entwickelt mit dem Ziel, zur Klärung der Frage nach dem Einfluß der Bodentopo- 
graphie auf den Verlauf der Strömungen im Meer und in der Atmosphäre beizutragen. Dabei 
wird der Einfluß der Sekundärströmung, welche durch das Absinken der Corioliskräfte in 
der bodennahen Reibungsschicht hervorgerufen wird, mitberücksichtigt. Von dem Einfluß 
der Erdkrümmung und von weiteren für den Verlauf ozeanischer und atmosphärischer Strömun- 
gen wesentlichen Einflüssen wird abgesehen. 

Die mathematische Formulierung der Theorie führt zu verhältnismäßig einfachen Grund- 
gleichungen, die zunächst im zweidimensionalen Fall (Strömung über langgestreckte Gebirgs- 
kämme etwa) allgemein in vollständiger Strenge gelöst werden. Die hergeleiteten einfachen 
Ausdrücke für Geschwindigkeitsverteilung, Isobaren- und Stromlinienverlauf gehen aus jenen 
hervor, die man ohne Berücksichtigung der Reibung erhält, indem man in den letzteren 
überall die jeweilige Bodenerhebung Z durch eine in einfacher Weise aus £ hervorgehende, 
gedachte Bodenerhebung [; ersetzt. Die Spiegeleinsenkung ist in allen praktisch interes- 
sierenden Fällen für den Strömungsverlauf bedeutungslos. 

Die Größenordnung der für die Strömungsvorgänge charakteristischen Parameter 

R H\ 

B läßt 
sich bei dem heutigen Stand unseres Wissens über die Gesetze turbulenter Strömung ab- 
schätzen. Auf Grund dieser Abschätzung kann man schließen, daß der Einfluß der Sekundär- 
strömung auf den Verlauf der Hauptströmung bei den praktisch vorkommenden Verhältnissen, 
sowohl im Meer wie in der Atmosphäre, im allgemeinen eine völlig untergeordnete Rolle 
spielt; er kann jedoch eventuell in Grenzfällen stark in Erscheinung treten und den Verlauf 
der Strömung entscheidend beeinflussen. 


(insbesondere der für den Reibungseinfluß entscheidenden Dimensionslosen G = 
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Nachtrag: 


Nach Abschluß dieser Arbeit bin ich der Aufforderung nachgekommen, in einem zu- 
sammenfassenden Bericht kurz über die neueren Arbeiten aus dem Kaiser Wilhelm-Institut 
für Strömungsforschung, Göttingen, zum Problemkreis der atmosphärischen und ozeanischen 
Strömungen zu berichten. In diesem Bericht sind auch die Grundgedanken und einige Er- 
gebnisse der obigen ausführlichen Mitteilung veröffentlicht worden. Siehe: 


[8] H.Görtler: Neuere Beiträge zur Dynamik atmosphärischer und ozeanischer Strömungen. Naturwiss. 
29. Jahrg. (1941), S. 473—479. 
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Zur rechnerischen Auflösung einer Gleichung vierten Grades. 
Von Helmut Heinrich in Breslau. 


Für eine reelle Gleichung vierten Grades wird aus der Tatsache, daf ihre 
kubische Resolvente stets eine positive reelle Wurzel besitzt, ein einfaches, schnell 
arbeitendes Verfahren zur gleichzeitigen Bestimmung aller ihrer Wurzeln 
(auch der komplexen) entwickelt und ein Rechenschema dafür bereitgestellt. 


Hat man von einer Gleichung höheren Grades auch die komplexen Wurzeln zu be- 
stimmen, so wendet man in der Regel das Graefte-Verfahren an. Daß man bei diesem zu- 
nächst nur die absoluten Beträge der Wurzeln erhält, muß in Anbetracht der Überlegungen, 
die dazu notwendig sind, um daraus die Vorzeichen und Real- und Imaginärteil zu gewinnen, 
ohne Zweifel als ein Mangel bezeichnet werden, der sich vor allem dann bemerkbar macht, 
wenn die Rechnung, wie es in vielen Betrieben der Fall ist, von mathematisch ungeschulten 
Kräften durchgeführt werden soll. Darüber hinaus schleichen sich erfahrungsgemäß leicht 
Vorzeichen- und Stellenwertfehler ein, die eine dauernde umständliche Kontrolle der Rechnung 
erforderlich machen '). 

Es soll deshalb in den folgenden Ausführungen wenigstens zur numerischen Auflösung 
von Gleichungen vierten Grades mit reellen Koeffizienten (Gleichungen fünften Grades lassen 
sich olıne Schwierigkeit durch Abspaltung einer stets vorhandenen reellen Lösung darauf 
zurückführen) ein Verfahren angegeben werden, das als wesentlich einfacher und leistungs- 
fähiger als das Graeffe-Verfahren bezeichnet werden darf und auch von Hilfskräften ohne 
weiteres vollständig durchgerechnet werden kann. Es liefert bei einem nur geringen Auf- 
wand an Arbeitszeit sämtliche vier Wurzeln, auch die komplexen, auf einmal, 
und zwar unmittelbar einschließlich Vorzeichen und Real- und Imaginärteil. Die zu leistende 
Rechenarbeit besteht im wesentlichen in der Bestimmung einer — von vornherein als vor- 
handen feststehenden — positiven reellen Wurzel einer Gleichung dritten Grades (der kubi- 
schen Resolvente der verkürzten Gleichung vierten Grades) und in der Ausführung einiger 
weniger elementarer Rechnungen nach einem festen Schema. 


Das Verfahren. 
Lautet die gegebene Gleichung 


so bringen wir sie zunächst dadurch, daß wir die linke Seite mittels Hornerschen 
Schemas nach Potenzen von 


entwickeln, auf die verkürzte Form 
Aus ihren Koeffizienten a, b, e bilden wir die Zahlen 
und stellen mit ihnen die kubische Gleichung 


auf. Da ihr absolutes Glied y negativ ist, hat sie stets mindestens eine positive reelle Wurzel. 
Unsere erste und wichtigste Aufgabe besteht darin, eine solche positive Lösung von (5) zu 
bestimmen, am besten, wenn mehrere vorhanden sind, die größte, weil man bei ihr am ehesten 
eine gute Konvergenz des Newtonschen Näherungsverfahrens erwarten kann, vor allem, 
wenn man sich ihr von rechts her nähert. 

Mit der Bestimmung dieser einen Wurzel z der kubischen Gleichung — die Rechnung 
erfolgt wieder mit Hilfe des Hornerschen Schemas — ist die Hauptarbeit geleistet. Denn 
kennt man z, so hat man nur noch 


1 
und 
.....0 
zu berechnen und erhält mit den vier möglichen Vorzeichenkombinationen in 
..... 


alle vier Wurzeln der verkürzten Gl. (3) und aus ihnen mit Hilfe von 2) Din von (1). 


1) Buange-König: Numerisches Rechnen, Berlin 1924, S. 164 ff. 
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Zweifellos ist es ein Vorzug des Verfahrens, daß sich die Genauigkeit, mit der die eine 
Wurzel 2 berechnet ist, unmittelbar auf alle vier Wurzeln der Gleichung vierten Grades 
überträgt. 

Die Grundlage des Verfahrens. 


Das Verfahren beruht auf einem Satz der allgemeinen algebraischen Theorie der Gleichung 
vierten Grades, der besagt, daß man die Wurzeln der verkürzten Gleichung vierten Grades (3) 
durch Wurzelausdrücke, die die Lösungen oder wenigstens eine Lösung der kubischen Re- 
solvente enthalten, ausdrücken kann’). Daß diese Tatsache bisher vom praktischen Rechner 
nicht ausgenutzt worden ist, mag daran liegen, daß ihre in den Lehrbüchern der Algebra 
meist in sehr abstrakten Zusammenhängen gegebene Darstellung natürlicherweise auf die 
praktisch-mathematischen Bedürfnisse keine Rücksicht nimmt und daher das Ergebnis in einer 
für die numerische Behandlung zunächst nicht geeignet erscheinenden Form liefert. 

Nach einer schon von Euler angegebenen Methode kann man s durch u+o-+-w er- 


setzen und w so bestimmen, daß +mw=—2a und urw=—Jb ist. Dann kann 
man zeigen, daß «*, vo’, w? die Wurzeln der kubischen Resolvente (5) sind. Sind sie bekannt, 
so hat man daraus «, v, w mit den geeigneten Vorzeichen — es muß urw=—b sein! — 


zu bestimmen, und die vier Wurzeln von (3) sind durch 


(9) 


gegeben. Daraus kann man dann weiter herleiten, daß, wenn u?” eine Wurzel der Resolvente 
ist, sich alle vier s; durch sie in der Form 


ausdrücken lassen. Daraus läßt sich übrigens umgekehrt ausrechnen, daß 


1 
ist, womit in Evidenz bestätigt ist, daß man mit (10) alle Wurzeln von (3) erfaßt hat. 
Die Resolvente (5) hat nun wegen des negativen absoluten Gliedes y = — b? stets eine 


positive reelle Wurzel «?, und daher kann in (10), wo man ja nur eine Wurzel der Resolvente 
braucht, « reell und positiv angenommen werden. Diese für die Zwecke des nume- 
rischen Rechnens entscheidende Tatsache ist bisher offenbar wenig beachtet worden. 
Der Grund dafür liegt wohl nicht zuletzt darin, daß man bei der oben skizzierten bisher üb- 
lichen Ableitung von den drei an Stelle von s neu eingeführten Unbekannten «, ® und w 
annehmen muß, daß sie auch komplex sein können. Verzichten wir aber auf die dabei er- 
reichte, an sich sehr schöne Symmetrie und setzen mit nur zwei Unbekannten & und n, die 
wir jedoch als reell voraussetzen, 


....... (MB). 


Den einen Freiheitsgrad, den wir (gegenüber zweien bei der anderen Methode) durch den 
Ansatz (12) gewonnen haben, benutzen wir dazu, um (13) in die beiden Gleichungen 


(14,) 


aufzuspalten. Dabei haben wir beim zweiten Teil auf den Faktor Y» verzichten können. 
Denn ist 7=0, so reduziert sich (13) wieder auf f(&$)=0, und wir sind wieder am Ausgangs- 
punkt angelangt. Nach der Aufspaltung in (14) bleibt trotzdem der Fall 7=0 von Bedeutung. 
Er besagt dann, daß f(£) und f’(£) gleichzeitig verschwinden und somit & eine mehrfache 
Wurzel ist. 

& können wir von Null verschieden annehmen, da &=0 b=f’(0)=0 nach sich zieht 
und für s? eine mit elementaren Mitteln lösbare rein quadratische Gleichung liefert. 

Setzen wir nun aus (14,) 


000 


so geht (3) über in 


2) z. B. Bieberbaceh-Bauer: Algebra, Leipzig 1928, insbes. S. 212f. 
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Rechenplan für die Gleichung +Ax&°+Bz?+Cz+D-0. 


I. Verkürzung der Gleichung und Berechnung der Koeffizienten der Resolvente. 


. . b 
a 2a a? 
a=2a y=—b 


II. Bestimmung der Wurzel der kubischen Resolvente +a2’+ßz+y=0. 


zo 1 a y 
ö, 
Ö, 
Ö, v 
Ö, Ö, . 
Ö, 
v 
us 
Ö, Ö, 


III. Berechnung der vier Wurzeln der gegebenen Gleichung aus z. 


2 

A, z+2ua 

4, 
A, A, 

A, A+b 

; 


_ 
| 
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in (14,) ein, so erhalten wir nach einfacher Rechnung 
1 
—4— 2a a? e)+ (16). 


Die bereits mehrfach hervorgehobene Tatsache, daß R(z) = 0 mindestens eine positive Wurzel 
besitzt, zeigt uns, daß R(42°?)=0 mit reellem & erfüllbar ist. Man braucht dazu nur für eine 
solche positive Wurzel z der Resolvente in Übereinstimmung mit (6) 4° —=z zu setzen. Die 
in (10) angeschriebenen vier Wurzeln s von (3) erscheinen dann in der besonders einfachen 
Form (8), auf der wir unser Verfahren aufgebaut haben. 


Praktische Durchführung und Zahlenbeispiel. 

Die Rechnung erfolgt zweckmäßig nach dem beigefügten Plan, der bei Kenntnis des 
Hornerschen Schemas ohne weiteres verständlich ist. Das Newtonsche Verfahren zur 
Verbesserung eines Näherungswertes z, von 2 wird in der Form verwendet, daß die linke 
Seite der Resolventengleichung fortlaufend nach Potenzen des Fehlers entwickelt wird. Das 
hat den Vorteil, daß man auch bei der Bestimmung von z auf mehrere Dezimalen weitgehend 
mit dem Rechenschieber auskommt. 

Im folgenden ist die vollständige Rechnung für die Auflösung der an anderer Stelle®) 
mit dem Graeffe-Verfahren behandelten Gleichung 


+52 +52 -10=0 
wiedergegeben. 
1 —3 2 7 
3 er 12 
49 
I. 13 49 
3 48 
1 16 
3 12 
4 
7 28 u 1/28 — 0,036 
0,086 0,086 0,2583 1017 
7,036 28,253 
0,086 0,255 
7,072 
7,108 28,508 0,017 0,017: 28,508 — — 0,0006 
- 0,0006 — 0,01 — 0,426 — 0,017 
7,107. 28.082 0,000 
z—= 3.0354 
A, 1,3054 
A, +iye 0,8711 - 0,8711 
A, A,A, 0,9019 — 0,9019 
A, 4+7 7,9019 6,0981 
A + VA 1505895 —1,50595 1,3229 — 13229 
8 | | 08711315095 | 04518 | 2,140 
| s+1 1,8711 +1,50595 | 1,4518 — 1,1940 299 


3) Willers: Methoden der praktischen Analysis, Berlin 1928, S. 209 ff. 
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Ergänzung zur Theorie der Grenzlinien 
adiabatischer Potentialströmungen. In mei- 
nem Aufsatz!) über dieses Thema (vgl. insbesondere 
Abschnitt 4) konnte der Beweis, daß die durch die 
Gl. (24) und (14) definierte Grenzlinie eine Rückkehr- 
kante der Stromlinien ist, nicht in voller Allgemein- 
heit zeführt werden: denn das benötigte Nıchtver- 
schwinden des Ausdruckes (30) mußte im Grunde 
zusätzlich vorausgesetzt werden. Ein Verschwinden 
von (30) blieb, wenn auch ausnahmsweise, immerhin 
möglich. 

Im weiteren Verlauf meiner Untersuchungen (vgl. 
die Abschnitte 5 und 6) hatte sich die grundlegende 
Bedeutung der Tatsache herausgestellt, daß der 
Grenzhodograph kein charakteristischer Hodograph 
ist. Diese Eigenschaft kann man leicht an den 
bekannten Beispielen von Grenzlinien bestätigen. 
Da ich diese Tatsache aber nicht als zusätzliche 
Voraussetzung einführen wollte, hatte ich einen 
allgemeinen Beweis dafür angesetzt. Dabei war 
die Übertragbarkeit der Richtungsübereinstimmung 
zweier Kurven aus dem Geschwindigkeits- in den 
Ortsplan benutzt. Diese Übertragbarkeit gilt auch 
noch auf dem Grenzhodographen, wie man nach 
(14) ausrechnen kann, es sei Br daß neben der 

ya 


Grenzgleichung gleichzeitig ist. In 
pi 


diesem Falle nimmt nämlich av den unbestimmten 


Wert . an. Da dieser Ausnahmefall bei dem er- 


wähnten Beweis hineinspielt, so ist die eingeschla- 
zene Beweisführung zu bemängeln. 

Es ist bemerkenswert, daß die beiden aufge- 
zeigten Lücken in der allgemeinen Theorie der 
Grenzlinien sich im Zuge einer und derselben Über- 
legung schließen lassen. 

Wir beginnen mit der Neufassung des Beweises 
dafür, daß der Grenzhodograph kein charakteristi- 
scher Hodograph ist. In den Gleichungen (35a) 
und (35b) waren als notwendige Bedingungen für 
das Zusammenfallen des Grenzhodographen mit 
einem charakteristischen Hodographen entweder 


ia ” . (35 a) 
oder 

— ya 

d = . (35 b) 


erkannt worden. Von da ab schließen wir jetzt so 


Der Richtungstangens des Grenzhodo- 


graphen wird unmittelbar aus der Grenzgleichung 


weiter. 


1 
F 
d 
dgq 
ausgerechnet. Es ist 
(z N) q 2 
—)» — - ’ 
ya a? ur 
1 1 


1 


Die erste Möglichkeit (35a) läßt sich nun sofort 
ausschließen: denn (35a) bedeutet konstantes y (9,9). 


1) Z. angew. Math. Mech., Bd. 21 (1941), S. 140 bis 152. 


definiert also eine Stromlinie. Die Grenzlinie 
kann aber als Machsche Welle nicht mit einer 
Stromlinie zusammenfallen. Setzt man nun (35e) 
unter Benutzung von (35d) in (35a) ein, so kann 
man die eben erkannte Unmöglichkeit auch so 
formulieren, daß 


1 1 
yaayat a? q 


|, 2 1 | 
Rz) 
(« —1 2 


wird. Die linke Seite dieser Ungleichung ist aber 
mit (30) identisch. Da nunmehr das Nichtver- 
schwinden von (30) bewiesen ist, so ist gemäß der 
Bemerkung am Eingang dieser Mitteilung jetzt 
ohne jede Einschränkung gezeigt, daß die Grenz- 
linie eine Rückkehrkante der Stromlinien ist. 

Zum Ausschluß der zweiten Möglichkeit (35b: 
setzen wir (35c) unter Berücksichtigung von (35d) 
in (35b) ein. Danach würde 


1 
2yaa ya 


2 ya: 1 


erforderlich sein. Das dritte Glied der linken Seite 
verschwindet wegen der Grenzgleichung. Die beiden 
ersten Glieder werden nach der Tschapligin- 
schen Differentialgleichung (15) umgeformt, so daß 
sich schließlich als Bedingung 


1 
q? 


q +) +2 
ergibt. Da die linke negativ ist, 


sc ist ein Widerspruch aufgewiesen und damit auch 
die zweite Möglichkeit ausgeschlossen. Es gilt 
daher allgemein, daß der Grenzhodograph kein 
charakteristischer Hodograph ist. 
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Die Grundgleichungen der Segeltheorie 
unter Berücksichtigung der Segelneigung 
in Längs- und Querrichtung. 

Die Prandtlsche Tragflügeltheorie hat auch auf 
dem Gebiet des Segelns zu einer neuen Betrach- 
tungsweise der Vorgänge am segelnden Boot ge- 
führt. Im Zusammenhang damit wurde die Ein- 
führung dimensionsloser Beiwerte für Auftrieb und 
Widerstand vorgenommen, die in der Segeltheorie 
‚unächst die Aufstellung elementarer Gleichungen 

für Vor- und Quertrieb ermöglichten. 


NE F- 
(F= Segelfläche: q= Staudiuck). 


F- 


Die Beiwerte c, und cs ergeben sich aus den Be- 
ziehungen 
Cy sin a® — a* 
und 
= Ca 608 a® + sin a*, 


o* ist der Winkel zwischen relativem Wind unil 
Fahrtrichtung. Diese Gleichungen hätten allgemeine 
Gültigkeit, wenn Vor- und “Quertriebskräfte und 
etwa außer diesen entstehende Kräfte unabhängig 
von der wechselnden seitlichen Bootsneigung (Krän- 


XUM 
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gung) wären. Dies ist aber nicht der Fall. In 
Wirklichkeit treten, je nach Anlage des Segels und 
Krängung des Bootes veränderliche aufwärts, oder 
abwärts gerichtete Kräfte hinzu. Ihre Mitberück- 
sichtigung ist daher erforderlich. Sie geschieht 
nun derart, daß man eine Querebene durch den 
relativen Wind annimmt und auf diese Ebene das 
Segel innerhalb einer senkrecht zur Kraftresultie- 
renden des Segels liegenden Ebene gekippt denkt. 
In der Querebene des relat. Windes liegt nunmehr 
die Projektion der Resultierenden (der Auftrieb) in 
einem Winkel +» zur Horizontalen. je nach 
ursprünglichen Anlage und späterer Krängung des 
Segels, während der Widerstand senkrecht zur 
Querebene steht. Bei ebenen unendlich glatten 
Flächen vollzieht sich dies Kippen in der Ebene 
der Fläche selbst. Die weiter angegebenen Be- 
ziehungen gelten für diesen Fall streng. Handelt 
es sich aber um gewölbte Flächen, oder Profilsegel 
mit oder ohne Verwindung, so wird hierbei die Ver- 
nachlässigung einer durch die Krängung entstehen- 
den, oder sich verändernden „Pfeilstellung“ in den 
Kauf genommen. Eine rechnerische Berücksichti- 
sung der letzteren ist nicht möglich: es läßt sich 
aber auf Grund vorhandenen Versuchsmaterials und 
Überlegungen theoretischer Art sagen, daß ihr 
Einfluß in den meisten Fällen vernachlässigbar 
klein ist. Man darf daher annehmen, daß die Rieh- 
tung der Resultierenden im Raume durch das Kip- 
pen des Segels in einer zu ihr senkrechten Ebene 
unverändert bleibt. 

Die Grundgleichungen nehmen auf Grund trigono- 
metrischer Rechnungen nunmehr folgende Form an: 


Cy = Cu sin a* — COS a* 
und 
= Ca 608 a* + sin a*, 


Für gilt 
tg a’ 

= 
Hierin bedeuten « Winkel zwischen Bootsmittel- 
linie und relativem Wind, «’ Krängungswinkel und 
v Mastneigung nach achtern. Pf’ ist der Winkel 
zwischen Bootsmittellinie und der Schnittlinie der 
Mastbaum-Ebene des Segels mit der horizontalen 
Ebene. 

1 


cos a’ 
tg v 


ist der Winkel der Resultierenden des Segels mit 
der Lotrechten auf derjenigen Sehne des Segels, 
auf welche die Segelpolare bezogen ist (in der 
Regel ist dies der Baum). 

Zur Berechnung des Quermomentes bedient man 
sich eines Beiwertes 


608 (a’ — y) cosy" 


cs’ 


y berechnet sieh aus tgy=tg3-sinr. 

Der Auftriebsbeiwert des Segels muß jeweils auf 
Grund der gegebenen Segelstellung nach der Polare 
des Segels erst festgestellt werden, da sich bei 
Krängung des Bootes der Einfallswinkel ö des 
Windes zur Mast-Baum-Ebene ändert. Hierbei ist 


sinö=yY 1— cos? - sin? (a — y) - sin (a — #). 


Nach der Polare wird auch der Winkel $ für Winkel 
ö berechnet, und gilt positiv bei Neigung nach 
vorne gegenüber der Lotrechten auf der Bezugs- 
sehne. 


Ist $# konst. der günstigste Einstellwinkel des 
Segels in aufrechter Lage, so ergeben sich beim 
Krängen des Bootes Vortriebsverluste. Ihre Ur- 
sache liegt einmal im Abfall des Auftriebes in- 
folge Verringerung des Einfallswinkels des Windes 
und andrerseits in einer Verringerung der horizon- 
talen Komponente desselben wegen Neigung des 
Auftriebes um <&n nach abwärts. Je nach der 
Polare des Segels und der gewählten Einstellung 
desselben, sowie der vorhandenen Mastneigung 
nach achtern, wirken sich beide Effekte verschieden 
aus. Es besteht daher generell kein Zusammen- 
hang mit dem cos. des Krängungswinkels, wie 
bisher angenommen wurde. Von den es und r’- 
Werten gilt analoges. Bild 1 und 2 zeigen dieses 
auf Grund der Polare eines Hochsegels, die im 
Windkanal gewonnen wurde. Mastneigung nach 
achtern war nicht vorhanden. Die Kurven sind 
errechnet. Aerodynamisch bessere Aggregate, Profil- 
segel usw. verhalten sich ähnlich: die Vortriebs- 
verluste sind — immer bezogen auf die günstigste 


Stellung in aufrechter Lage — kleiner. 
S / 
N / 
/ 
05 


Bild 1. Vortriebsverluste zweier Hochsegel in Abhängig- 
keit von der Krängung. 


Segell_]r=0; 
Einfallswinkel des Windes zum Boot: «= 30°. 


PER 
N 

0% 

cos 
02 Cy 
@-0 70° 20° 30° 40° 50° 


Bild 2. Vor- und Quertriebskräfte eines Hochsegels in 
Abhängigkeit vom Krängungswinkel «*. 


Einfallswinkel des Windes zum Boot: @a=30%; \r=0. 


Die Beibehaltung der für aufrechte Bootslage 
besten Einstellung des Segels auch in geneigter 
Lage erweist sich allgemein als am vorteilhaftesten. 
Ein schärferes Anholen des Segels würde zwar X ö 
wieder vergrößern, jedoch würde dem nunmehr 
vergrößerten Krängungsmoment eine größere Nei- 
gung entsprechen, und der Verlust wäre im allge- 
meinen noch größer. Entscheidend ist hierbei der 
Verlauf der Stabilitätskurve des Bootes. Von 
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„weifelhaftem Vorteil für den Vortrieb und bedenk- 
lich im Hinblick auf die Stabilität erweist sich die 
Verwendung selbsttätig sich auf konstanten An- 
stellwinkel einstellender Segel; letzteres besonders 
deswegen, weil die natürliche Abnahme des Quer- 
momentes bei der Krängung sich als größer er- 
weist, als allgemein angenommen wird. 

Eine andere Arbeitsweise zeigen Segel mit nach 
achtern geneigter Drehachse (Mast). Der Anstell- 
winkel ö wächst während der Krängung bis zu 
sin a’ 
tgv 
entspricht, um von dort langsamer abzunehmen, 
als bei Segeln mit lotrechter Achse. Dieser Ver- 
lauf der Winkeländerung gibt eine Erklärung da- 
für, daß die Vorsegel länger vollstehen, als das 
Großsegel. Die Vortriebsverluste sind wesentlich 
geringer und erreichen den Wert Null ebenfalls bei 


tg 


liegt. 

Es besteht demnach ein Zusammenhang zwischen 
Krängungswinkel und günstigster Neigung des 
Mastes, oder der Vorsegelachsen nach achtern, und 
es erweist sich, daß rankere Fahrzeuge auch ge- 
neigtere Segelachsen erfordern: besonders die unter 
starker Belastung stehenden Vorsegel. 

Die Momente solcher Segel sind bei gleichen Vor- 
trieben geringer als bei lotrechten Masten. — Das 
Verhältnis von Quer- zu Vortrieb wird durch Krän- 
gung bis etwa 30° kaum beeinflußt, darüber hin- 
aus verschlechtert. Anch diese nachteilige Wir- 
kung vermindert sich sehr wesentlich durch Nei- 
gung des Segels nach achtern. 

Vorliegende Betrachtungen „eben den Grund- 
gleichungen der Segeltheorie eine allgemeine Form. 
Sie gestatten nunmehr die Arbeit des Segels in 
allen Lagen zu übersehen und sind auch unver- 
ändert auf die Bewegung der Bootsflosse (Kiel) 
im Wasser anwendbar. 

Damit ist ein umstrittenes Problem der quantita- 
tiven Behandlung zugänglich gemacht und aufge- 
hellt worden. Mit Hilfe von Segel- und Boots- 
polare ist es jetzt möglich, das gesamte Kräfte- 
spiel am segelnden Fahrzeug zu verfolgen. 
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einer Schräglage, die der Beziehung tg3’—= 


—, da nunmehr die Kraftresultierende 


Über genäherte Lösungen von Differen- 
tialgleichungen mit Hilfe von Minimal- 
ansätzen '). 

Für ein einfaches Beispiel aus der Balkenbiegung 
wird die Näherungslösung erst mit Hilfe des Mini- 
malansatzes der potentiellen Energie, dann mit 
Hilfe eines anderen gelegentlich angewandten 
Minimalansatzes gefunden und die Ergebnisse mit 
der richtigen Lösung verglichen. Die z. T. sonder- 
baren Ergebnisse werden zum Schlusse untersucht 
und erklärt. 

Als Beispiel ist ein Balken genommen, der in 
den Punkten —1 und = +1 gestützt ist. 
Die Belastung q (x) sei symmetrisch, also q (x) — 
4-8); für” gegebene Belastung wird die Durch- 


biegung gesucht: ist das negative Biege- 


dr 
moment ay ist die Belastung: 
i astung: 
4 
=0 


Die gewählte Näherungslösung sei /(r). 
der potentiellen Energie gibt: 


einem Dresdener Vortrag. 


Der zweite Minimalansatz ist: 


1 
f _ 
Die Randbedingungen sind: 


Min 


(Durchbiegung gleich null) 
_ _, 
dz?’ dr? 


(Biegemomente gleich null). 


Nur die richtige Lösung macht das Integral der 
(3) zum Minimum, das den Wert null hat. 

Die Untersuchung ist im weiteren aus Symme- 
triegründen nur für die rechte Balkenhälfte 0<r 
< 1 durchgeführt. Als Balkenbelastung ist gewählt 
nach Bild 1: 


mit 
für 
für 
und 
für 0OEr<09—: 
für 
(X) = 0 für 09<z<si1. 


03 - 
7 


Gewählte Belastung. 


Bild 1. 


Wir lassen e—0 gehen und erhalten an der 
Stelle x — 0,9 die Einzelkraft von der Größe — 1.2. 
Es werde erst die genaue Lösung ermittelt, um 
später hiermit die Näherungslösungen vergleichen 
zu können. Für die Teilbelastung g,(z) erhält 


man die Durehbiegung yı: für diese gilt: 


für 0<r<09 


2 
0522 0,495 


dx? 
Y,= 0,0416... — 0,2475.2? + 0.2058375 


und für 


Y = 0,15 23 — 0,45 2° + 0,1215. 2 + 0,1785. 
Für die Teilbelastung %, (2) erhält man mit 


die Durchbiegung y,, wobei negatives 
Moment zu bestimmen ist: 
für 0o<r<09 

d’y 


= 0.06 2? — 0,0598 
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und für 09<r<i1 


d? 
—12r+12 
— 022 40,61? — 0,486.2+ 0,086 . 
Die gesamte Durchbiegung = ist 


überall positiv. In der Mitte für ©=0 erhält man 
Ym = 0.2058375 — 0.0598 — 0,1460375.. 


Als Näherung für beide Minimalansätze neh- 
men wir die Durchbiegung der Belastung qı (X), 
die wir mit einem noch frei wählbaren Beiwert 
multiplizieren: 


= a (0,0416... 2? 
für 0O<r<0,9 und 

f(z) = a(0,15.2° — 
für 09<Sr<1. 


— 0,2475 2? + 0,2058375) 


0.45 2? -+ 0,1215. + 0,1785) 


Der Minimalansatz der potentiellen Energie gibt. 


wenn für q(@)= qM(X) + die betreffenden 
Werte eingesetzt werden und nur über die 


Strecken integriert wird. wo diese nicht gleich 
null sind: 
0,9 1 
\a? (0, 502 —0 495)? (0. 
v 0,9 
0,9 
2 \ a (0,0416... 2 — 0,2475 2? + 0,2058375) dr 


—2 a (0,0416... 

+ 
— 0,2475 2? + 0,2058375)d Min. 

Das letzte Integral erstreckt sich über die 


Länge e, der Klammerausdruck zibt an der Stelle 
x — 0,9 den Wert 0,0327. Das Integral gibt: 


— 1,2 


a-0,0327 — a-0,03924 . 
Aus dem Minimalansatze folgt hiermit: 
a?-0,130032 — 2 a-0.130032 + 2 a-0.03824 = Min. 
Das günstigste @op erhält man durch Differen- 
tiation nach a: 
0,03924 


Aopt = 1— 0.130032 1 


Mit diesem Werte wird die Durchbiegung in der 
Mitte: 


/m (0) = 0,2058385 — 0,062116 — 0,143721. 


Dieser Wert ist etwa 2% kleiner als der genaue 
Wert. 

Der zweite Minimalansatz gibt mit derselben 
Näherungsfunktion für die Durchbiegung 
0,9 0,9 1 
0 v 

Das dritte Integral hat auf die Bestimmung von 


— 0,27935. 


PL 


@opt keinen Einfluß; das Integral \ ag(z)dx gibt 


den Wert «(09 — 123 = — 0,3 a. 
Hiermit erhält man: 
0.9 a? + 0,6a = Min 
und durch Differentiation 
Ap=—083... 


Die Durchbiegung nach dieser Nährungsmethode 
nimmt überall negative Werte statt posi- 
tiveran. 

Um diese sonderbaren Ergebnisse zu erklären. 
formen wir die Gl. (2) und 3) um. 

Gl. (2) gibt, wenn wir darin auf Grund zwei- 
maliger partieller Integration 

d?y d?/(z) 
setzen und das für das Minimum belanglose Glied 


2 
dx hinzufügen: 


dr 
( 

f| ir dze=Min . (2a). 
0 
Gl. (3) hingegen gibt mity (z) = Tr, 

1 

des dre=Min. (3a). 

Im ersten Falle wird durch «@,: das Fehler- 


integral der zweiten Ableitung, im zweiten Falle 
das Fehlerintegral der vierten Ableitung zum Mi- 
nimum gemacht. Am günstigsten wäre ein Mini- 
malansatz, bei dem das Fehlerintegral der Funktion 
selbst zum Minimum gemacht wird, wenn man 
sich für die Funktionswerte interessiert. Je höhere 
Ableitungen man annähert, desto größere Abwei- 
chungen werden im allgemeinen für die Funktions- 
werte eintreten. 

Es lassen sich auch Beispiele konstruieren, bei 
denen der Minimalansatz der potentiellen Energie 
ganz falsche Näherungswerte für y gibt. 

Bei Benutzung der Minimalansätze ist darum 
stets Vorsicht geboten. Nur der Ausdruck, der 
selbst zum Minimum gemacht wird, wird bei gün- 
stigen Näherungsansätzen gut angenähert. 
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Neuer Beweis des Entwicklungssatzes der 
Vektoralgebra. 

1. A. Föppl| (in seiner „Einführung in die Max- 
wellsche Theorie der Elektrizität“, Leipzig 1894). 
und nach ihm viele andere, haben den Entwick- 
lungssatz der Vektorrechnung bewiesen, indem sie 
auf die Zerlegung nach Komponenten eines recht- 
winkligen Dreibeins zurückgriffen. Diesem mehr 
rechnerischen Verfahren hat LagallyM. in seiner 
ausgezeichneten Darstellung der Vektorrechnung 
(Vorlesungen über Vektorrechnung, Leipzig 1934) 
einen anderen anschaulicheren Beweis entgegen- 
gestellt, der indessen doch, wie sein Verfasser 
selbst bemerkt, auch nicht kürzer ist. Im folgenden 
gebe ich die Herleitung des Satzes in einer neuen 
Form, die wohl allen billigen Ansprüchen an Ein- 
fachheit und Anschaulichkeit entsprechen dürfte. 

2. In dem Entwicklungssatz: EX (AX B)—=(B-E) 
A—(A-E)B ist nämlich ein ebenso einfacher wie 
schöner Umklappungssatz der ebenen Geometrie 
enthalten, der allerdings gar nicht oder nur wenig 
bekannt sein dürfte und dessen Beweis ich daher 
voranstelle, wobei ich mich auf Bild 1 beziehe. 

Von € sind auf die Geraden MM’ und NN’ die 
Lote CA und CB gefällt. Über OC als Durchmesser 
ist der Kreis durch A und B gelegt. Das Dreieck 
AOB ist um die Winkelhalbierende des Winkels 
AOB umgeklappt. In dem umgeklappten Dreieck 
4'OB’ steht nun — und das ist der Inhalt des Um- 
klappungssatzes — A’B’ zu OC senkrecht. 


Zum Beweis beachte man, daß (wegen der Um- 
klappung und wegen des Umfangswinkelsatzes): 
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Bild 2 


Bild 3. 


ist, so daß das bei B rechtwinklige Dreieck OBC 
mit Dreieck OE4’ in den Winkeln übereinstimmt, 
weshalb <7OE4’=90° ist und 4’B’ zu OC senk- 
recht steht. 

3. Da nun die Formeln: EX(AXB)=(B-E) 
A—(A-C)B, in allgemeinen Vektoren, und c x 
(aX b)=(b-c)a— (a-c) b,in Einheitsvektoren, durch 
Division oder durch Multiplikation mit ABC (dem 
Produkte der Beträge der drei Vektoren) in ein- 
ander überführbar sind, genügt es, die Formel für 
Einheitsvektoren herzuleiten. 

Wir setzen a, b und c (Bild 2) bei O zu einem 
Dreibein zusammen und errichten dann, senkrecht 
zur Ebene von a und b, den Vektor d=aXb, 
sowie, zur Ebene von c und # senkrecht, den Vektor 
E=cX(aX b), der in die Ebene von a und b 
fällt. Es gilt dann, unter den Winkel von a 
nach b, unter ö aber den von c nach # verstanden: 


und 
= sind. 


4. Wie können wir E konstruieren? 
Wir fällen von C, dem Endpunkte von ct, folgende 
Lote: CC, auf die Ebene von a und b; 
CC’ auf 9; CA auf a und schließlich CB auf b. 
Dadurch erhalten wir C,4 senkrecht zu a und 
C, B senkrecht zu b. Nun klappen wir noch Drei- 
eck OAB um die Halbierende des Winkels x um 
Im umgeklappten Dreieck O4’B’ ist A’B’, als 
Vektor, gleich E=cxX(aXb): denn erstens ist 
4’B’ (im richtigen Sinne) zur Ebene von c und 9 
senkrecht, wie man mittels des Umklappungssatzes 
erkennt; zweitens ist auch dem Betrage nach 
4AB=|E=sing sind; denn es ist: 

0C,=CC' =0C sin ö= sin Ö. 

Da aber OC, Durchmesser des um O4AC,B be- 

schreibbaren Kreises ist, besitzt die 4’ B’ gleiche, 


y gegenüberliegende Sehne AB die Länge sin ö 
sing=|E. 

Andererseits ist, weil die Projektion eines Ein- 
heitsvektors auf einen anderen durch ihr Punkt- 
produkt dargestellt wird: 


0B=0B=(b.e), 
04=04'=(a:c). 


Daher ist: Vektor 04’=04Ab=(a-c)b und Vektor 
OB'=0OBa=(b-c) a, so daß wir finden: 


AB=6=cX(aXb) 
OB — 04 = (b-)a— b. 


5. Sollte C, etwa nicht in den Winkel oder 
in dessen Scheitelwinkel fallen, so zieht man zum 
Beweise noch die zu 4’B’in bezug auf O spiegelige 
4" B” zu Hilfe, die man an Stelle von 4’B’ erhält, 
wenn man Dreieck AOB um die Halbierende der 
Nebenwinkel von p umklappt. Diese Strecke 4” B” 
ist A’ B’ entgegengesetzt gleich und stellt — € vor. 
An Bild 3 läßt sich das nachprüfen. 


Aschaffenburg. 0. Hofmann. 259. 


Linearverbindung kompressibler ebener 
Strömungsfelder. Die Potentialgleichung der 
stationären ebenen kompressiblen Strömung ist 
nichtlinear. Man kann sie jedoch durch Übergang 
zum Geschwindigkeitsbild linearisieren 1) und hat 
auf dieser Grundlage verschiedene Integrations- 
methoden zur Ermittlung kompressibler Strömungs- 
felder entwickelt ?2). Im folgenden wird auf ein sich 


1) Vgl. A. Busemann: Z. angew. Math. Mech. Bd. 17 
(1937), S- 73 bis 79 und die dort angegebenen früheren 
Arbeiten von P.Molenbroek, A.Tschapligin, D.Ria- 
bouehinski und B. Demtschenko. 

2) A. Busemann:a.a. C. Ferrari: Aerotecnica 
Band 18 (1938), S. 400 bis ur: Hsue-Shen-Tsien: J. Aero- 
naut. Sciences Bd. 6 (1939) S. 399 bis 407; F. Ringleb: 
Z. angew. Math. Mech. Bd. 20 (1940) S. 185 bis 198. 
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unmittelbar aus der Linearisierung ergebendes In linearverbundenen Strömungsfeldern ent- 


geometrisches Verfahren hingewiesen, welches aus 
zwei bekannten Strömungsfeldern eine 1-parametrige 
Menge „linearverbundener“ neuer Strömungsfelder 
zu konstruieren gestattet. Das Verfahren wird zur 
Herleitung eines neuen Beispiels eines Strömungs- 
feldes mit Durchgang durch die Schallgeschwindig- 
keit 3) 4) 5) benützt. 


1. Linearisierung der Potential- 
gleichung. 


Wir schreiben die Potentialgleichung der statio- 
nären ebenen kompressiblen Strömung in der Form 


2 


Hierbei bedeuten die unteren Zeiger partielle Dif- 
ferentiationen und die Schallgeschwindigkeit c ist 
eine Funktion von gx?—+gy?: diese Funktion ist 
durch die Adiabatengleichung ?=P (0) festgelegt 
und lautet beispielsweise für vollkommene Gase 


—1 


(k = Verhältnis der spezifischen Wärmen, ® = 
Ruhedruck, 0, = Ruhedichte). 

Der Übergang von der Strömungsebene (x, y- 
Ebene) in die Ebene des Geschwindigkeitsbildes 
(u, v-Ebene) erfolgt am einfachsten mittels der 
Legendre-Transformation 


v=yy, (3) 


Sie verwandelt die nichtlineare Potentialgleichung 
(1} in die lineare Differentialgleichung t) 


2 2 2 
Du v (1 Pu (4), 


wobei jetzt c eine Funktion von u? + v2 ist. 


um 


2. Konstruktion linearverbundener 
Strömungsfelder. 


Wegen der Linearität der Gleichung (4) lassen 
sich ihre Lösungen linear überlagern, d. h. mit ®, 
und 8, ist auch 


ein Integral von (4). Aus (3) und (5) kommt 
) 
. . (6) 
(u, v) | 


und man hat somit folgende Konstruktionsvor- 
schrift (Bild 1, 3): 

2 kompressible Strömungsbereiche mit derselben 
Adiabatengleichung seien punktweise 
einander zugeordnet, derart, daß in entsprechenden 
Punkten P,, P, die Geschwindigkeitsvektoren 
nach Größe und Richtung übereinstimmen. Es er- 
gibt sich ein neuer kompressibler Strömungsbereich 
mit derselben Adiabatengleichung, wenn man alle 
Verbindungsstrecken P,P, in einem festen Ver- 
hältnis A (A — beliebiger Parameter) innen oder 
außen teilt und jedem Teilpunkt P wieder den Ge- 
schwindigkeitsvektor iw zuweist. Die 3 Strömungs- 
felder der Punkte P,, P,. P bezeichnen wir als 
„linearverbunden“. 


3) W. Tollmien: Z. angew. Math. Mech. Bd. 17 (1937) 
S. 117 bis 136. 

4) F. Ringleb: a.a.O. und Deutsche Mathematik Bd. 5 
(1941) S. 377 bis 384. 

5) H.Görtler:Z.angew. Math. Mech. Bd. 19 (1939) S. 325 
bis 337 und Bd. 20 (1940) S. 254 bis 262. 


sprechen sich: 

a) die Isotachen und Isoklinen, d. h. die Kurven 
konstanten Betrags bzw. konstanter Richtung des 
Geschwindigkeitsvektors; 

b) die Unterschall- und Überschallbereiche und 
die etwaigen Durchgangslinien durch die Schall- 
geschwindigkeit; 

e) die Machschen Netze der Überschallbereiche. 
Die Machschen Kurven sind hierbei parallelbezogen, 
d. h. sie haben in entsprechenden Punkten parallele 
Tangenten. Die Machschen Netze sind reziprok zu 
einem gemeinsamen charakteristischen Netz im 
Geschwindigkeitsbild in ähnlicher Weise, wie ver- 
schiedene Lagepläne demselben reziproken Kräfte- 
plan zugeordnet sein können. 

Die Stromlinien und die Potentiallinien ent- 
sprechen sich bei der punktweisen Zuordnung 
linearverbundener Strömungsfelder im allgemeinen 
nicht. Die Potentialwerte sind gemäß (6) linear 
kombiniert. 

Kompressible Strömungsfelder können im Über- 
schallbereich .„Grenzlinien“ haben, auf die erstmals 
W.Tollmien?) hingewiesen hat und die neuerdings 
unter dem Namen .„Stoßlinien“ von F. Ringleb*®) 
untersucht wurden. An diesen Grenzlinien flanscht 
das Strömungsfeld von einem Blatt auf ein anderes 
Blatt der doppelt überdeckt zu denkenden Ebene 
über. Bei der Linearverbindung von Strömungs- 
feldern durch Zuordnung von Punkten mit gleichem 
Geschwindigkeitsvektor entsprechen die Grenz- 
linien einander i. a. nicht, sondern einer Grenzlinie, 
an der beispielsweise 2 Blätter des P,-Feldes zu- 
sammenhängen, ist im P,- und im P-Feld eine 
Kurve zugeordnet, deren beiderseitige Umgebungen 
zum gleichen Blatt gehören. 

Wenn sich die kompressible Strömung mit 0 — 
eonst (Machsche Zahl M — 0) zur inkompressiblen 
Strömung spezialisiert, werden die Isotachen-Iso- 
klinen-Netze isotherm und die Zuordnung durch 
l:inearverbindung geht in eine konforme Abbildung 
über. 

3. Beispiele. 

Wir benützen die Linearverbindung zur Kon- 
struktion von Strömungsfeldern mit Durchgang 
durch die Schallgeschwindigkeit. Dabei legen wir 
die Adiabatengleichung der idealen Gase zugrunde 
und bezeichnen den Geschwindigkeitsbetrag mit w. 
Dann ist nach (2) die Schallgeschwindigkeit des 


ruhenden Gasese,— Vk»,/0, und der bei Expan- 
sion ins Vakuum eintretende Höchstwert der 
Strömungsgeschwindigkeit 


Co 1 bzw. W max = (7), 
wenn wir zur Abkürzung W — w/c, setzen. Dem 
Durehganz durch die Schallgeschwindigkeit ent- 
spricht 

4 
k+1 
a) Linearverbindung von Quellströmung 

und Wirbelströmung. 


Die kompressible Quellströmung ergibt 
durch Integration von 


bzw. W* = (8). 
sich 


Radiusvektor r und Geschwindigkeitsbetrag w bzw. 
W sind verknüpft durch die Beziehung 


P7°: 
"= rin \/ 9). 
k+1 


- 
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wobei der kleinste Radiusvektor 7,n beliebig ge- 

wählt werden kann. Der Kreis r=[r„in, dem 

i entspricht, ist sowohl Grenz- 


linie als auch Durchgangslinie vom Unterschall 
zum Überschall; der Außenbereich r > rin ist 
doppelt überdeckt, vn zwar hat man im Unterschall 


I<W<W*= % Überschall 


1 
<< 
Kreisinnere "< rin läßt sich die -kompressible 
Quellströmung nicht fortsetzen. 


Für die kompressible Wirbelströmung folgt aus 


r dr 


und im 


In das 


=0 


die Beziehung 


(10) 


mit dem-wieder beliebigen BE Radiusvektor 
Für W=W* = 


er als Durchgangslinie vom 


erhält man den 


Kreis r*=1r,in \ 
Unterschallbereich r >r* zum Überschallbereich 
>r >ruin: Ebenso wie die Quellströmung läßt 
sich auch die kompressible Wirbelströmung nicht 
in das Kreisinnere r< rin fortsetzen. 

Die Isotachen sowohl der Quellströmung als auch 
der Wirbelströmung sind die konzentrischen Kreise 
r — econst. Bei der Linearverbindung einer Quell- 


und Wirbelströmung werden je 2 Isotachen-Kreise 
(9) und (10) mit gleichem W einander zugeordnet, 


Bild 3. Linearverbindung von Quelle und Tollmien- 


scher Spiralströmung. 


w 
a Bild I. Linearverbin- 
dung von Quelle und 
Wirbel 
/ De 
/ 


N 


Bild 2. Stromlinien der Linearver- 
bindung von Quelle und Wirbel. 


/ 


wobei die entsprechenden Radienvektoren zuein- 
arder senkrecht stehen (Bild 1). Die von den Teil- 
punkten P erzeugten Isotachen sind wieder Kreise 
und man erhält als Ergebnis der Linearverbindung 
m „Wirbelquellströmung“, wie sie bereits von 

J. Taylor®) angegeben wurde. Auch bei der 
Wirbelquelle kann die Strömung ins Innere eines 
Kreises r — rin nicht fortgesetzt werden. Dieser 
Kreis ist hier jedoch lediglich Grenzlinie und liegt 
innerhalb des Kreises r — r*, an dem die Strömung 
vom Unterschall zum Überschall übergeht. 


In Bild 2 ist ein Ausschnitt einer durch Linear- 
verbindung von Wirbel und Quelle entstandenen 


6) G.J. Taylor: Z.angew. Math. Mech. Bd. 10 (1930) S. 341 


EIS 


Bild 4. Stromlinien der Linearverbindung 
von Quelle und Tollmienscher Spiral- 
strömung. 
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Wirbelquellströmung dargestellt. Die mit Pfeilen 


gekennzeichneten, durch Drehung um ausein- 
ander hervorgehenden spiralenförmigen Kurven 


sind Stromlinien; die Pfeilrichtung weist vom 
Unterschall zum Überschall. Die Grenzlinie Kreis 
FT=1rnin) ist stark ausgezogen, die Durchgangs- 
linie vom Unter- zum Überschall (Kreis <= 
ist strichpunktiert. Die im ringförmigen Überschall- 
bereich verlaufenden gestrichelten Kurven sind 
Machsche Linien. Sie beginnen an der Durchgangs- 
linie r — r* senkrecht zu den Stromlinien mit Rück- 
kehrpunkten und berühren die Grenzlinie r — Tin - 


b) Linearverbindung von Quellströmung 
und Tollmienscher Spiralströmung. 


Die Tollmiensche Spiralströmung 3) hat die kon- 
zentrischen, durch Drehung auseinander hervor- 
zehenden logarithmischen Spiralen 

. (11) 
— Polarkoordinaten, a — tg = konstante Nei- 
unge der Spiralen gegen die Radienvektoren) als 
lsotachen. Längs jeder Spirale 7 — const ist nicht 
nur der Betrag W der Geschwindigkeit, sondern 
auch der Neigungswinkel « des Geschwindigkeits- 
vektors gegen die Spirale konstant. W und «a sind 
(laher Funktionen der Veränderlichen 7: vgl. hierzu 
W. Tollmien 3). 


Bei der Linearverbindung einer Quellströmung 
mit einer Spiralströmung wird je ein Isotachenkreis 
(9; und eine Isotachenspirale (11) einander zuge- 
ordnet, wobei entsprechende Radienvektoren den 
konstanten Winkel #-a miteinander einschließen 
(Bild 3). Die von den Teilpunkten P erzeugten 
Isotachen sind spiralenförmige Kurven, welche in 
unendlich vielen Windungen einen Kreis (— in 
Bild 3 punktiert —) asymptotisch umschlingen. Das 
durch Linearverbindung sich ergebende Strömungs- 
feld stellt ein neues Beispiel einer kompressiblen 
Strömung mit Durchgang durch die Schhllgeschwin- 
digkeit dar; während in den bekannten Beispielen 
der Quellströmung. Wirbelströmung, Wirbelquell- 
strömung und Tollmienschen Spiralströmung die 
Isotachen stets von den Stromlinien unter kon- 
stanten Winkeln durchsetzt werden, sind bei dem 
neuen Beispiel die Isotachen nicht Isogonaltrajek- 
torien der Stromlinien. 


In Bild 4 ist ein Ausschnitt der dureh Linear- 


verbindung von Quell- und Spiralströmung ze- 
wonnenen neuen Strömung dargestellt, wobei die 


Bezeichnungen dieselben sind wie im analogen 
Bild 2. 
Aachen. R. Sauer. 254 


Funktionsgewichte in der 
Ausgleichsrechnung. 

Kurze Inhaltsangabe: 

Bei der Errechnung irgendwelcher Zustands- 
größen eines Vorganges aus den zugehörigen Meß- 
werten läßt sich eine Genauigkeitssteigerung da- 
durch erzielen, daß man mehr Meßgeräte ansetzt, 
als es zur einmaligen Bestimmung der betr. Zu- 
standsgröße erforderlich ist. Es werden Methoden 
entwickelt zur Errechnung eines Mittelwertes 
unter Berücksichtigung aller Einzelmessungen und 
zur Abschätzung seiner Genauigkeit. 


Aufgabenstellung. 

Zur Bestimmung eines (z. B. zeitlich veränder- 
lichen) Funktionswertes x seien im gleichen Zeit- 
punkt rn Größen aı, @e, ... a, gemessen. Es sei 
weiterhin angenommen, daß der Funktionswert 


schon durch n—1 dieser Meßgrößen vollkommen 


bestimmt ist. Durch Auswahl von je n—1 dieser 
n. Argumente a, erhält man jeweils einen Funk- 
tionsausdruck 


) 


Es gibt n Funktionen z,. Wären nun die Argu- 
mente fehlerlos, so würden alle Werte z, überein- 
stimmen; da aber die a, mit Meßfehlern behaftet 
sind, werden sich im allgemeinen auch rn verschie- 
dene Werte x, ergeben. Um aus diesen den wahr- 
scheinlichsten Wert z zu ermitteln, sind die ein- 
zelnen ., derart mit Gewichten zu belegen, daß 
die zuverlässigsten Funktionswerte den größten 
Einfluß auf die Mittelwertbildung = haben. 

Definition der Funktionsgewichte. 

Die mittleren quadratischen Fehler der einzelnen 
Meßgrößen a, seien mit zu, tr, . . 
der mittlere quadratische Fehler der Funk 
tion 7, sei 

Das Einzelgewicht eines Meßwertes a, aus 
einer Reihe gleichartiger (aber im Allgemeinfall 
nicht gleichwertiger) Messungen kann, wie in der 
Ausgleichsrechnung gezeigt wird'), dem irgendwie 
festzustellenden Quadrat des Einzelfehlers um- 
gekehrt proportional gesetzt werden. Man kann 
also schreiben 


k? 
k? 
Gy . . . . . . . . . . (2), 


wobei g, das Gewicht der Meßgröße a,., G, das 


Gewicht des Funktionswertes 7, darstellt. 

Der wahrscheinlichste Funktionswert x ist dann 

Gu 


Die Beziehung zwischen G, und 9, wird durch 


das Gaußsche Fehlerfortpflanzungsgesetz her- 
gestellt. Es gilt 
und aus (1) (2) (4) folgt 
da,/ 


v 


Die in diesen Formeln auftretenden Argu- 
mentgewichte g, (gegenüber den Funktions- 
gewichten G,) sind von den MeBßfehlern der ein- 
zelnen Argumente a, abhängig, gestatten also, 
2. B. die Meßgenauigkeit eines Instrumentes zu 
berücksichtigen. Die Argumentgewichte g, werden 
hier als von vornherein gegeben angenommen. 


Bereehnung der Ableitungen. 


Für die Bestimmung der Funktionsgewichte 
nach (5) werden sämtliche Ableitungen aller Funk- 
tionswerte 7, =f*’(...) benötigt. Die Berech- 
nung dieser Ableitungen kann nun z. B. durch 
Differentiation der einzeln aufgestellten Funktions- 
ausdrücke geschehen; im folgenden soll ein Ver- 
fahren angegeben werden, das in der Praxis die 
Bildung der Ableitungen oft wesentlich erleichtert. 


1) Hütte Bd. 1. 26. Aufl., S. 1056. Kohlrausch: 
Praktische Physik. Lpz., 1935, 17. Aufl.. S.21.— Helmert: 
Die | der Methode der kleinsten 
Quadrate, 1907, ?. Aufl. 
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Dabei wird nicht mehr von sämtlichen Einzel- 
funktionen © ausgegangen, sondern es wird nur 
eine von ihnen, etwa die am einfachsten gebaute, 
zur Bildung der Ableitungen herangezogen. Dies 
sei die Funktion 


(Bei Benutzung nur einer Grundfunktion ist der 


Index am Funktionszeichen f fortgelassen.) Zwi- 
schen den Argumentwerten 4,, @s, ... 4, besteht 
nun eine Bedingungsgleichung 


Eine solche folgt z. B. aus der Gleichsetzung ver- 
schiedener Funktionswerte x, oder auch aus geo- 
metrischen Beziehungen. Aus der Bedingungs- 
eleichung (7) und der Ausgangsfunktion (6) er- 
geben sich umgekehrt alle Funktionswerte 7,. 
Daher müssen sich auch die Ableitungen sämt- 
licher Funktionswerte aus dieser einen Ausgangs- 
funktion ünd der Bedingungsgleichung folgern 
lassen. Aus Gleichung (7), die ja alle Größen a, 
enthält, läßt sich jedes «a, als Funktion der übri- 
gen Argumentwerte errechnen. Es sei nach (7) 
insbesondere 


Setzt man diesen Ausdruck in (6) ein, so erhält 
man eine Funktion, die nicht mehr von a, ab- 
hängt, also den Funktionswert z,. Es ist somit 
p(dı, 

Ay), dx +1, .. 


2, =lla,4,,... 


Ay), 


und für die Ableitungen von x, nach «a, erhält 


man 
of 
1 a, day 
da, 
day ar 
d Ay day 


Mittlerer Fehler des Mittelwertes. 


Entscheidend ist nun zu wissen, wie genau der 
nach (3) (5) (8) gefundene Mittelwert bestimmt ist, 
Es ist also der mittlere Fehler des Mittelwertes x 
festzustellen. Das Gewicht des Mittelwertes 7 ist 
Gı+@:+...+6,: für den mittleren Feh- 
ler M des Mittelwertes gilt also ent- 
sprechend (2) 


G+6-+ +6 
in Verbindung mit (2) folglich 


Für die Praxis ist es wichtig, den mittleren 
Fehler des Mittelwertes möglichst herabzudrücken. 
Dies erreicht man nach (9) durch Vergrößerung 
der Anzahl der zur Berechnung herangezogenen 
Werte. Bei Langbasis-Luftaufnahmen mit Kino- 
theodoliten, wofür dieses Verfahren entwickelt 
wurde, ist diese Genauigkeitssteigerung besonders 
bemerkenswert. da hierbei ja eine Erhöhung der 
Genauigkeit durch mehrfach wiederholte 
Messungen nicht möglich ist. Bei Hinzunahme 
eines dritten Basispunktes für die Vermessung 
eines Luftzieles lassen sich — wie man leicht 


nachreehnen kann — zur Berechnung des Mittel- 
wertes der Höhe des Luftzieles sogar 19 Funk- 
tionswerte heranziehen. 


Erweiterung. 


Sind die Funktionswerte 2 durch weniger als 
n —1 Größen, etwa durch n— m Größen a, be- 
stimmt, so lassen sich aus den n Meßwerten «a, 


je nach der Art der Aufgabe bis zu (») verschie- 


dene Funktionswerte errechnen. Die Funktions- 
werte seien auch in diesem Falle stets mit den 
Indizes der nicht in der betreffenden Funktion 
vorkommenden Meßgrößen «@, gekennzeichnet: also 
z. B. für n=7 und m=3 


Die Beziehung (5) geht dann über in 
1 
\ 
— day 
Die Ableitungen errechnen sich bei einer Aus- 
gangsfunktion 
und m Bedingungsgleichungen 
Fir 
1, 43, n (7%), 


m) 


wie durch H. Molitz bewiesen wurde, nach 


u 


da, da, da; 
da, da, da; R 
89, 
arwarwarw arm 
dA, da, da, da, | 


worin 4 die Hauptunterdeterminante der großen 


Determinante ist, also der Faktor von ar. Der 
v 
wahrscheinlichste Funktionswert ist jetzt 


Ergebnis. 

Als Ergebnis der Untersuchungen wird fest- 
gestellt: Falls Überbestimmung eines Funktions- 
wertes durch eine Bedingungsgleichung vorliegt, 
ist sein wahrscheinlichster Wert gegeben durch 


1 
mit @, = Gewicht des Funktionswertes 7, aus (5) 


und (8), wobei 9, das durch die Beziehung "= 


definierte Gewicht des v-ten Meßwertes darstellt und 
4, als mittlere Abweichung der Meßgröße a, ge- 
geben ist. Im Allgemeinfall gelten entsprechend 
die Beziehungen (3*) mit k=1, (5*) und (8*). 
Berlin. H. Knobloch. 239 
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BUCHBESPRECHUNGEN 


Festschrift der Mathematischen Gesellschaft in 
Hamburg zu ihrem 250jährigen Bestehen 1940. 
Erster Teil: Geschichte der Gesellschaft von 1890 
bis 1940. (Mitteilungen der Mathematischen Gesell- 
schaft in Hamburg, Bd. VIII, Teil 1.) 73 S. Ham- 
burg 1940. 

Del. Zweiter Teil: Abhandlungen (Mitteilungen der 
Mathematischen Gesellschaft in Hamburg, Bd. VII, 
Teil 2.) 239 S. Hamburg 1940. 

Der erste Teil der Festschrift enthält eine Ge- 
schichte dieser allein schon durch ihr Alter be- 
merkenswerten Mathematischen Gesellschaft über 
den Zeitraum von 1890 bis 1940. Der zweite Teil 
enthält fachwissenschaftliche Aufsätze, hauptsäch- 
lich aus dem Gebiet der Reinen Mathematik. 

Dresden. Rellich. 29 


WILHELM BLASCHKE, Mathematik und 
Leben. (Hamburger Mathematische Einzelschrif- 
ten, 27. Heft 1940.) 13 S. Leipzig und Berlin 1940, 
Verlag B. G. Teubner. Preis geh. 1 M. 


Vergnügte Plauderei über die Mathematik und 
ihren Leumund. 


Dresden. Rellich. 269 


Dr.-Ing. EUGEN GRÜNAGEL, Flüssigkeits- 
bewerung in umlaufenden Radial- 
rädern. (V DI-Forschungsheft 405, Beilage zu 
„Forschung auf dem Gebiete des Ingenieurwesens“ 
Ausg. B. Bd. 11, November/Dezember 1940.) 21 8. 
m.52 Bildern und 1 Zahlentafel. Berlin 1940, VDI- 
Verlag G.m.b.H. Preis 5 M. 

Der erste Teil des Heftes berichtet über Ver- 
suche, die Strömung in umlaufenden Radialrädern 
sichtbar zu machen. Das geschieht in der be- 
kannten Weise durch eine Versuchseinrichtung, in 
der Wasser mit freier Oberfläche strömt, auf die 
Aluminiumpulver gestreut wird. Die Aufnahmen 


der Strömungen in einer größeren Zahl von 
Pumpen- und Turbinenrädern sind sehr gut ge- 
lungen und geben ein deutliches Bild der Verhält- 


nisse. Leider sind überhaupt keine Zahlenangaben 
über die Drehzahlen und Strömungsgeschwindig- 
keiten gemacht; es ist anzunehmen, daß diese 
Werte außerordentlich klein waren. Der Verfasser 
nimmt zwar an, daß der Einfluß der Reynoldsschen 
Zahl nicht sehr groß sein wird. Jedoch sind die 
Unterschiede gegenüber einer normalen Strömungs- 
maschine so groß, daß selbst die qualitative. Über- 
einstimmung mit den hier auftretenden Strömungen 
nicht als sicher angenommen werden kann. 

Im zweiten Teil der Arbeit wird der Versuch ge- 
macht, die Strömungsvorgänge rechnerisch zu ver- 
folgen, wenn in der Beschauflung Ablösungen «in- 
treten und sich Totwasserräume bilden. Es wird 
für die Rechnung angenommen, daß im Schaufel- 
kanal ein Teil von einem geordnet strömenden 
Flüssigkeitsstrahl und der übrige Teil von Tot- 
wasser ausgefülit ist und daß zwischen diesen Be- 
reichen eine eindeutige Grenze besteht. In bezug 
auf die Reibung in dieser Grenzfläche wird die be- 
kannte Annahme über die turbulente Scheinreibung 
gemacht, die formal mit der Reibung an festen 
Wänden übereinstimmt. Es wird gezeigt, daß mit 
(diesen Annahmen die Leistungsübertragung und 
die Verluste in der Strömungsmaschine rechnerisch 
verfolgt werden können, vor allem auch in den 


Spalten zwischen Laufrad und Leitapparat. An 


dieser Stelle bestanden bisher noch ganz erheb- 
liche Unklarheiten. Über die Größe der turbulenten 
Scheinreibung liegen noch keine brauchbaren Meß- 
ergebnisse vor, so daß eine zahlenmäßige An- 
wendung der Rechnung nicht möglich ist. — Die 


Einheitliehkeit der Arbeit leidet etwas unter den 
zu häufigen Hinweisen auf andere Arbeiten des 
Verfassers. 


Augsburg. E. Sörensen. 273 

Dr. MOHAMED HILAL, Beiträge zur 
Theorie und Berechnung von Balken- 
brücken aus Eisenbeton. (Mitteilungen 
aus dem Institut für Baustatik an der Eidg. Tech- 
nischen Hochschule in Zürich, herausgegeb. v. d. 
Professoren Dr. M. Ritter und Dr. F. Stüssi, Mit- 
teilung Nr. 11.) 208 S. m. 128 Abb. Zürich und 
Leipzig 1940, Verlag A. G. Gebr. Leemann & Co. 
Preis brosch. 9 M. 

Die Arbeit behandelt einige Fragen, deren 
Klärung für die Beurteilung des Spannungs- und 
Verschiebungszustandes der Balkenträger für 
Eisenbetonbrücken nützlich ist. Die erste Unter- 
suchung beschäftigt sich mit dem Trägerrost und 


der lastverteilenden Wirkung der Querträger 
zwischen mehreren Hauptträgern. Der Verfasser 
ergänzt eine kurze Schilderung der bekannten 


Lösungen dieses wichtigen bautechnischen #Pro- 
blems durch eine Ableitung der statischen Zu- 
sammenhänge nach den Ansätzen, die zur Be- 
rechnung des durchlaufenden Trägers auf elastisch 
senkbaren Stützen verwendet werden. Der zweite 
\bscehnitt behandelt die Eigenspannungen des 
durchlaufenden Trägers mit elastisch drehbaren 
Stützen bei Temperaturänderung und beim Schwin- 
den des Baustoffs. Hierbei wird auch der Unter- 
schied des Schwindmaßes im bewehrten und unbe- 
wehrten Bereich des Trägers beachtet. Darauf 
folgen Untersuchungen über den Einfluß des ver- 
änderlichen Trägheitsmomentes mit verschiedener 
analytischer Formulierung auf die Biegemomente 
des beiderseits elastisch eingespannten Trägers. 
Zum Schluß beschäftigt sich der Verfasser mit 
dem Spannungszustand “des durchlaufenden Trägers, 
ohne die für die klassische Elastizitätstheorie 
übliche lineare Abhängigkeit zwischen Spannung 
und Dehnung des Betons anzunehmen. Die Auf- 
gaben besitzen wissenschaftliche Bedeutung. Sie 
sind zwar zum großen Teil bereits gelöst, von dem 
Verfasser jedoch eingehend behandelt und durch 
Tabellen, Zahlenrechnung und Beispiele erläutert 
worden, so daß auch das Studium dieser Ergeb- 
nisse nützlich ist. 


Dresden. K. Beyer. 293 
Vorspannung im Eisenbetonbau. 
Grundlagen, Ziel, Zweck und Anwen- 
dung. Beiträge von Prof. Dr.-Ing. L. PISTOR, 
Dipl. -Ing. R. OPPERMANN, Dr. Ing. W. PASSER, 
Öberbaurat Dr.-Ing. e. h. Dr. F. v. EMPERGER. 
mit einem Geleitwort von Prof. Dr.-Ing. A. KLEIN- 
LOGEL. 80 S. m. 68 Textabb. Berlin 1940, Verlag 
Wilhelm Ernst & Sohn. Preis brosch. 3,80 M. 


Die Festigkeit der Trägerroste und Stabwerke 
aus Eisenbeton wird im konstruktiven Ingenieur- 
bau nach den bekannten Ansätzen zur Navierschen 
Balkenbiegung mit versagender Zugzone des Be- 
tonquerschnitts nachgewiesen, da die ersten Zug- 
risse je nach der Art der Bewehrung oft bereits 
schon bei Bruchteilen der endgültigen Belastung 
eintreten. Es ist daher verständlich, daß schon 
seit längerer Zeit versucht worden ist, den ein- 
zelnen Bauteilen Druckkräfte durch eine Vorspan- 
nung der Rundeisen zu überlagern, um der Riß- 
gefahr und den damit verbundenen Nachteilen vor- 
zubeugen. Diese Bemühungen sind jedoch erst 
durch die Verwendung von Rundstahl mit hoher 
Streckgrenze erfolgreich gewesen, um bei der 


Vorspannung Dehnungen zu erzielen, die nicht 
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(dureh das Schwinden des Betons und das Kriechen 
des Baustoffs unter Last ausgeglichen werden. 
Diese mit der Vorspannung von Eisenbeton und 
mit der Verwendung von hochwertigem Baustahl 
verbundene theoretischen und konstruktiven Be- 
trachtungen werden in einzelnen Aufsätzen be- 
handelt, dig bereits in Fachzeitschriften veröffent- 
licht oder vorgetragen worden sind. L. Pistor, 
Müncher, gibt eine allgemeine Schilderung der ge- 
senwärtig verwendeten bautechnischen Methoden. 
. Oppermann, Frankfurt a. M., schildert die 
Versuche mit Spannbetonträgern der Bauart 
Freyssinet und die daraus gewonnene Er- 
kenntnis über die Verbundwirkung von Beton und 
Stahl. Darauf folgt die Beschreibung der Bauaus- 
führung einer Straßenbrücke mit unterspannten 
Balken durch W. Passer, Berlin, in der zur 
statischen Berechnung und der Ausrüstung der- 
artiger Eisenbetonkonstruktionen Stellung genom- 
men wird. Zum Schluß erörtert F.v.Emperger, 
Wien, Ziel und Zweck der Vorspannung im Eisen- 
betonbau und die gleichzeitige Verwendung von 
vorgespannten und nicht vorgespannten Stabein- 
lagen. Der Leser gewinnt mit diesen Aufsätzen 
eine gute Übersicht über den gegenwärtigen Stand 
der Technik und die Ziele, die der Eisenbetonbau 
gegenwärtig zur Behebung der Rißgefahr und zur 
Einsparung von Baustahl verfolgt. Entscheidend 
ist hier jedoch nicht die Theorie, sondern die 
Bauausführung. 

Dresden. K. Beyer. 2% 


Dr. phil. BRUNO TIEDEMANN, Über Boden- 
untersuchungen bei Entwurf und 
Ausführung von Ingenieurbauten. 
40 S. m. 22 Textabb. und 5 Tafeln. Berlin 1941. 
Verlag Wilhelm Ernst & Sohn. Preis brosch. 2 M. 

Nach der im Reichsgebiet baupolizeilich und bau- 
verwaltungsseitie kürzlich verbindlich eingeführten 
„Richtlinien für die zulässige Belastung von 
Flach- und Pfahlgründungen“ DIN 1054 sind bei 
allen größeren Ingenieurbauten eingehende Boden- 
untersuchungen vorgeschrieben und dabei die Zu- 
ziehung einer anerkannten Versuchsanstalt oder 
eines anerkannten Fachmannes empfohlen worden. 
Der Ingenieurgeologe der Preußischen Versuchs- 
anstalt für Wasser-, Erd- und Schiffbau, Dr. B. 
Tiedemann, behandelt in dem vorliegenden 
Heft — einem Sonderdruck aus der Zeitschrift 
„Die Bautechnik“ — die Art dieser Gemeinschafts- 
arbeit zwischen Entwurfsbearbeiter und Bauleitung 
einerseits und der Fachstelle andererseits. Ins- 
besondere wird der Teil der Untersuchungen be- 
sprochen, die bei Großbauten der Baustelle zu- 
fallen; das sind Bohrungen und Schürfungen, die 
nicht dem „Bohrmeister“ allein überlassen werden 
dürfen, die richtige Benennung der Bodenarten 
nach bodenphysikalischen Gesichtspunkten, wofür 
Anleitungen gegeben werden, und die sachgemäße 
Entnahme von ungestörten Bodenproben, sowie 
ihre Beschriftung, Verpackung, Aufbewahrung oder 
Versendung, um die sich ebenfalls ein geschulter 
Ingenieur kümmern muß. Für Großbaustellen wird 
die Einriehtung einer örtlichen Bodenprüfstelle für 
(die einfacheren Untersuchungen empfohlen und An- 
leitung auch dafür gegeben, um die Versuchs- 
anstalten nicht unnütz zu belasten. 

Das ausgezeichnete Heft ist aus der jahrzehnte- 
langen Erfahrung einer führenden Erdbauversuchs- 
anstalt entstanden und im richtigen Augenblick 
erschienen. Es bedarf deshalb keiner weiteren 
Empfehlung. 


Naumburg/Saale. Dr. Seheidig. 276 


GÜNTHER BRIEGLEB, Atome und Ionen. 
(Hand- und Jahrbuch der Chemischen Physik, her- 
ausgegeben von A. Eucken und K. L. Wolf, Bd. 2, 


Abschnitt IA.) 404 S. m. 150 Fig. im Text. Leipzig 
1940, Akademische Verlagsgesellschaft Becker & 
Erler, Kom.-Ges. Preis brosch. 38.50 M. 

Der Beitrag behandelt die allgemeinen Eigen- 
schaften der Atome (wie Isotopie, Loschmidtsche 
Zahl, Atom- und Ionenradien), die Elementarteil- 
chen (Wellen- und Teilehennatur der Materie und 
des Lichtes, Elektron, Positron, Proton, Neutron. 
Meson), die Atomkerne (Bau, Kräfte und Energien, 
Umwandlungen) und die Atomhülle (Einelektronen- 
systeme, Vektormodell, Atomterme). Wenn auch 
für die meisten Teilgebiete schon gute zusammen- 
fassende Darstellungen bestehen, so ist doch der 
vorliegende Beitrag zu begrüßen als eine Übersicht 
über die experimentellen Methoden und Ergebnisse 
des Gesamtgebietes der einfacheren Materieteilchen. 
Darüber hinaus bringt der Abschnitt über die 
Elementarteilchen viele wertvolle Zusamnienstel- 
lungen über Nachweis- und Meßmethoden, die ohne 
ihn nur mit großer Mühe zu erhalten wären. 

Leipzig. F. Hund. 271 


Dr. HANS GEORG MÖLLER, o. Prof. a.d. Univ. 
Hamburg, Grundlagen und mathema- 
tische Hilfsmittelder Hochfrequenz- 
technik (Lehrbuch der drahtlosen Nachrichten- 
technik Bd. I, herausgegeb. von Nicolai v. Kor- 
shenewsky und Wilhelm T. Runge). XVI+ 2938. 
m. 353 Abb. Berlin 1940, Verlag Julius Springer. 
Preis geb. 28,80 M. 

Der vorliegende Band bringt auf dem engen 
Kaum von noch nicht 300 Seiten eine ziemlich 
vollständige Theorie der Hochfrequenztechnik. 
Nur die ausgesprochen mathematische Behand- 
lungsweise erlaubt eine derart gedrängte und da- 
bei exakte Darstellung eines so umfassenden 
Wissensgebietes. Darüber hinaus sind in einem 
Anhang nicht nur die komplexe Rechenweise und 
die Vektorenrechnung als besonderes mathemati- 
sches Werkzeug des Hochfrequenztechnikers, son- 
dern sogar die wichtigsten Grundgesetze der Elek- 
trizitätslehre kurz umrissen. Der Stoff des Haupt- 
teiles beginnt mit dem einfachen Schwingungs- 
kreis und ‚führt über die gekoppelten Kreise zu 
den Filtern. Es folgt, als eigenstes Arbeitsgebiet 
des Verfassers etwas ausführlicher behandelt, die 
Physik der Elektronenröhre und ihre Anwendung 
als Verstärker, Gleichrichter und Generator, von 
der ursprünglichen Rückkopplungsschaltung bis 
zu den Barkhausen- und Habann-Generatoren. Ein 
dritter Abschnitt behandelt die Antenne und die 
Ausbreitungsgesetze der elektromagnetischen 
Wellen. Der gesamte Lehrstoff ist in über 200 
kleine Abschnitte und Unterabschnitte aufgeteilt. 
Die wichtigsten Ableitungen und Formeln sind je- 
weils durch kurze Zahlenbeispiele ergänzt. 

Das Gesamtlehrbuch wird neben dem vorliegen- 
den Einführungsband 5 weitere Bände umfassen, 
von denen der nächste, über Antennen und Wel- 
lenausbreitung, bereits erschienen ist, während 
die weiteren sich verteilen auf Elektronenröhren 
— Verstärker und Empfänger — Sender — Fern- 
sehtechnik. 

Nach wissenschaftlichen Lehrbüchern der draht- 
losen Nachrichtentechnik besteht seit langem ein 
dringendes Bedürfnis, denn nachdem vor Jahr- 
zehnten die ersten Lehrbücher der klassischen 
drahtlosen Telegraphie erschienen waren, fanden 
sich während der folgenden Zeit stürmischer Ent- 
wicklung kaum Autoren für eine zusammen- 
fassende Behandlung des stark angewachsenen 
Gebietes. Erst in letzter Zeit ist versucht worden, 
diese für die Ausbildung des Ingenieurnachwuch- 
ses sehr nachteilige Lücke auszufüllen. Der Sprin- 
gersche Verlag bringt sein Lehrbuch als Gemein- 
schaftswerk mehrerer Verfasser heraus. Daß diese 
Schriftenreihe durch den einführenden Band eines 
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namhaften Physikers und Lehrers eröffnet wird, 
der selbst die gesamte Entwicklung des Gebietes 
mit erlebt und mit beeinflußt hat, bedeutet einen 
vielversprechenden Anfang. Ebensosehr ist es zu 
begrüßen, daß dem ins Breite und Einzelne gehen- 
den Charakter des Sammelwerkes mathematische 
Kürze und Beschränkung auf das methodisch 
Wichtigste im Einführungsband die Waage hält. 
Diese schwierige Aufgabe ist dem Verfasser in er- 
staunlichem Maße gelungen. Für den Unterricht 
wird das Buch wertvolle Dienste leisten und dem 
Studierenden wie dem tätigen Ingenieur, beson- 
ders dem mathematisch begabten, von großem 
Nutzen sein. 


Dresden. Wolman. 321 


Dr. G. JOOS, Prof. an d. Univ. Göttingen, und 
Dr. TH. KALUZA, Prof. an d. Univ. Göttingen, 
Höhere Mathematik für den Prakti- 
ker, an Stelle einer 6. Aufl. des Lehrbuchs der 


Differential- und Integralrechnung von H. A. 
Lorentz. 2. verb. Auf. XI+368 S. m. 
85 Abb. Leipzig 1940, Verlag Johann Ambrosits 


Barth. Preis geb. 24,50 M. 


In der zweiten Auflage werden die Hyperbelfunk- 
tionen stärker berücksichtigt als bisher. neu hinzu- 
eekommen sind eine Skizze der Zylinderfunktionen 
(Besselsche Differentialgleichung, Reihenent- 
wicklung, Integraldarstellung) und eine Anleitung 
zur formelmäßigen Darstellung praktisch wichtiger 
Kurven. Im übrigen sind Licht und Schatten noch 
so verteilt wie in der ersten Auflage. An dem 
Licht haben sich die erfreut, die nur an den An- 
wendungen der Mathematik interessiert sind, auf 
den Schatten haben die Mathematiker hingewiesen, 
z.B. R.Rothe in der Besprechung der ersten Auf- 
laget. Im Vorwort der vorliegenden zweiten 
findet sich nun eine Bemerkung, die an die 
Rothesche Kritik anknüpft, aber an ihrem Kern 
vorbeigeht: Rothe hat keineswegs bemängelt. 
daß die Verfasser dem Praktiker Existenzbeweise 
vorenthalten wollen, sondern vielmehr, daß dem 
Leser ein in wesentlichen Punkten unvollständiger 

haliig vorgesetzt wird, und das ist wirklich eine — 
noch dazu ganz überflüssige — Unkorrektheit, die 
mit der Sünde, eine Näherungsformel außerhalb 
ihres Gültigkeitsbereiches anzuwenden, durchaus 
vergleichbar ist! 


Freiberg/Sa. 


„Beweis“ der Gleichung als stich- 


Grüß. 297 


1) Z. angew. Math. Mech. Bd. 18 (1938), S. 315 bis 316. 


Dr.-Ing. e. h. Dr. techn. e.h. FRITZ EMDE, Prof. 
a. d. Techn. Hochschule Stuttgart, Tafeln ele- 
mentarer Funktionen. XI+181 S. m. 
83 Textfig. Leipzig und Berlin 1940, Verlag B. G. 
Teubner. Preis geb. 12 M. 


In der dritten Auflage der allbekannten Funk- 
tionentafeln von Jahnke und Emde sind die 
elementaren Funktionen ausgeschieden worden. 
Die dadurch entstandene Lücke hat der Heraus- 
geber durch ein neues Werk geschlossen, das den 
elementaren Funktionen eine umfangreiche und 
vielfach neuartige Behandlung und Darstellung zu- 
teil werden läßt. Von den Inhaltserweiterungen 
gegenüber der zweiten Auflage der Funktionentafeln 
mögen hier die folgenden aufgeführt werden: Die 
Faktorentafel, die eine Zerlegung der natürlichen 
Zahlen bis zu 10% leistet; Lösungsverfahren für die 
Gleichungen zweiten bis vierten Grades mit aus- 
führlichen Anleitungen und graphischen Dar- 
stellungen; neue umfangreichere Tafeln für die 
Kreis- und Hyperbelfunktionen; Tafeln und Figuren 
zur Ermittlung der Stromverdrängung; Kurven- 


bilder der Tschebyschewschen Polynome: ein 
logarithmisches Diagramm der mechanischen Ein- 


heiten; und schließlich die sehr beachtenswerten 
Ausführungen über Zahlenrechnungen. Darüber 


hinaus aber wird dem Benutzer auch in den über- 
nommenen. und ausgestalteten Teilen noch vieles 
andere Neue und Wichtige dargeboten. 

Die Bezeichnungsweise ist sorgfältig gewählt: 
nur müßte es auf S. 160 unten folgerecht „Kehr- 
wert“ statt „Umkehrung“ heißen. 

Das Werk besitzt die bekannten Vorzüge der 
Funktionentafeln in erhöhtem Maße: dazu trägt 
auch das größere Format bei, das der Deutlichkeit 
der Formeln und der Figuren zugute kommt. Das 
Buch ist ein gediegenes Anschauungs- und Rechen- 


hilfsmittel auf allen Anwendungsgebieten .Jer 
Mathematik. 
Dresden. P. E. Böhmer. 288 
Grund- und Zeitfragen der Ver- 


sicherungswirtschaft und Versiche- 
rungswissenschaft. (Schriftenreihe des In- 
stituts für Versicherungswissenschaft der Uni- 
versität Köln, Bd. 1.) Herausgegeben von Prof. Dr. 
WaltherRohrbeck. 1248. Berlin 1941, Ver- 
lag Julius Springer. Preis brosch. 6,40 M. 

Das 124 Seiten starke Heft enthält die Vorträge, 
die bei der Eröffnung des Instituts für Versiche- 
rungswissenschaft an der Universität Köln gehalten 
wurden. Diese Vorträge lassen die Wandlungen 
erkennen, welche die Auffassungen über Gegen- 
stand und Ziel der behandelten Wissenschafts- 
cebiete im Laufe der letzten Jahre erfahren haben. 

Dresden. P. E. Böhmer. 28 


U. PESONEN, Messungen der Grund- 
linien Örebro und Värnamo in Schwe- 
den im Jahre 1938. (Baltische Geodätische 
Kommission, Sonderveröffentlichung Nr. 8.) 26 S. 
Helsinki 1941. 


Einer Anregung der Baltischen Geodätischen 
Kommission entsprechend wurden im Sommer 1938 
unter der Leitung von J. Bonsdorff zwei 
Grundlinien gemessen: die eine in der Nähe der 
Stadt Örebro, die andere südlich von der Stadt 
Värnamo. Die Beobachter waren U. Pesonen 
und P. Kalaja vom Finnischen Geodätischen 
Institut. 

Die Eichung der benützten 24 m langen Iavar- 
drähte erfolgte vor und nach der Messung der 
Grundlinien auf der bei Nummela in Finnland ge- 
legenen, 36 Drahtlängen großen Eichstrecke. Die 
Grundlinie bei Örebro ist etwas mehr als 7 km 
lang; diejenige bei Värnamo ist nicht ganz 7 km 
lang. Für die Messungen wurden acht Invardränte 
des Finnischen Geodätischen Instituts verwendet. 
Die mittleren Fehler der beiden rund 7 km langen 
Grundlinien liegen zwischen #1.0 und +1,5 mm. 


Stuttgart. P. Werkmeister. 2% 


Tafeln zur Übertragung geogra- 
phischer Koordinaten auf. dem inter- 
nationalen Erdellipsoide im Bereich 
350 bis 71° Breite. (Baltische Geodätische 
Kommission, Sonderveröffentlichung Nr. 9.) 53 8. 
m. Auszug. Helsinki 1941. 

Die von der Baltischen Geodätischen 
Kommission als Sonderveröffentlichung Nr. 9 
herausgegebenen Tafeln beziehen sich auf die geo- 
dätische Hauptaufgabe, bei der es sich um die 
Übertragung von geographischen Koordinaten auf 
dem Erdellipsoid mit Hilfe einer geodätischen Linie 
von bekannter Richtung und Länge handelt. Die 
Tafeln gelten für das Internationale Erdellipsoid 
mit der großen Halbachse gleich 6378388 m und 
der Abplattung 1:297. Unter den verschiedenen 
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Nachrichten 


Z. angew. Math. Mech. 
Bd. 21 Nr. 5 Okt. 1941 


Formeln für die Koordinatenübertragung wurden 
die „alten, gut bewährten und in den meisten 
Staaten wohl angewandten“ Formeln von 0. 
Sehreiber mit den von L. Krüger angege- 
benen Abänderungen gewählt. 

Stuttgart. P. Werkmeister. 29 


Schriftleitung folgende 


Ferner sind bei der 
Besprechung 


Bücher eingegangen (ausführliche 
bleibt vorbehalten): 


Internationalen 
Brückenbau und 
35 S. Zürich 


Mitteilungen der 
Vereinigung für 
Hochbau, Nr. 1, 20. Januar 1941. 
1941. 


B. N. DELONE und D. K. FADDEEW, Theo- 
rie der Irrationalitäten dritter Stufe. 
(Arbeiten des mathematischen Institutes B. W. A. 
Steklow XL) 340 S. Moskau und Leningrad 1940, 
Verlag Akademie der Wissenschaften, USSR. Preis 
30 Rubel. 


Dr €. KLOUCEK, Das Prinzip der fort- 
zeleiteten Verformung als Weg zur 
Ausschaltung der Unbekannten aus 
dem Formänderungsverfahren. 391 S. 
Berlin 1941, Verlag Wilhelm Ernst & Sohn. Preis 
brosch. 18 M. 


J. P. AREND, Die wirkliche Mechanik 
des Geschehens im Lebenden und 
Toten. 456 S. Zürich und Leipzig 1941, Rascher 
Verlag. Preis geb. 10 M. 


Dr. CLEMENS SCHAEFER, o.ö. Prof. der Physik 
a. d. Univ. Breslau Einführung in die 
Maxwellsche Theorie der Elektrizi- 
tät und des Magnetismus. 4. verb. Aufl. 
142 S. m. 33 Textfig. Berlin 1941, Verlag Walter 
de Gruyter & Co. Preis geb. 6 M. 


Prof. G. MAHLER, Physikalische For- 
melsammlung. Neubearbeitett von Prof. 
K. Mahler, Studienrat a. d. Schubartoberschule 
Aalen in Württemberg. 7. verb. Aufl. (Sammlung 
Göschen Bd. 136). 152 S. m. 69 Fig. Berlin 1941, 
Verlag Walter de Gruyter & Co. Preis geb. 1,62 M. 


Dr.-Ing. habil. FRANZ KARAS, Elastische 
Formänderung ‘und Lastverteilung 
beim Doppeleingriff gerader Stirn- 
radzähne. (VDI-Forschungsheft 406, Beilage zu 
„Forschung auf dem Gebiete des Ingenieurwesens“, 
Ausgabe B, Bd. 12, Januar/Februar 1941.) 233 S. 
m. 13 Bildern und 4 Zahlentafeln. Berlin 1941, 
VDI-Verlag G.m.b.H. Preis brosch. 5 M. 


Dr.-Ing. BRUNO ECK, Technische Strö- 
mungslehre. VII-+ 252 S. m. 278 Abb. Ber- 
lin 1941, Verlag Julius Springer. Preis geb. 13,50 M. 


NACHRICHTEN 


Arbeitstagung 
der Deutschen Mathematiker-Vereinigung in Jena. 

Vom 19. bis zum 24. Oktober 1941 fand in Jena 
eine überaus gut besuchte Arbeitstagung der 
Deutschen Mathematiker-Vereinigung statt. In der 
Gruppe „Angewandte Mathematik“ wurden die 
folgenden Vorträge gehalten: 

E. Hölder (Braunschweig): Grundsätzliches 
über die Behandlung der Anfangs- und Rand- 
bedingungen bei einer nicht linearen hyperboli- 
schen Differentialgleichung mit zwei unabhängi- 
gen Veränderlichen (an Beispielen der Gasdyna- 
mik), 

W. Damköhler (Jena) u. H. Grell (Augs- 
burg): Zur Integration der Differentialgleichung 
der stationären ebenen adiabatischen Gasströmung, 

E. Hopf (Leipzig): Mathematische Grundpro- 
bleme der Hydrodynamik, 


G. Schulz (Berlin): Wahrscheinlichkeiten ab- 


hängiger Ereignisse, 

H. Gebelein (Eßlingen); Anwendung der 
Kombinatorik auf den Aufbau der Statistik, 

H. Hornich (Wien): Zur Theorie des Risikos, 

K. Ludwig (Hannover): Zur Biegung des 
Kreisringausschnittes, 

A. Walther (Darmstadt): Neuere Verfahren 
der praktischen Mathematik, 

L. Viötoris (Innsbruck): Zur Theorie des In- 
tegraphen, 

Th. Zech (Darmstadt); Zum Abklingen von 
nicht linearen Schwingungen, 

L. Collatz (Karlsruhe): Natürliche Schritt- 
weite bei numerischer Integration von Differen- 
tialgleichungssystemen, 

L. Collatz (Karlsruhe): Eigenwertberechnung 
bei Gleichungen der Form L(f) (f)=0, 

K.Bögel (Blankenfelde): Gedämpfte harmonische 
unter Einfluß der trockenen Rei- 
ung, . 

Bückner (Berlin: Näherungslösungen der 
linearen Differentialgleichung erster Ordnung. 


Auch in den anderen Gruppen wurden eine 
Reihe von Fragen behandelt, die für den für die 
Anwendungen interessierten Mathematiker von 
Wichtigkeit sind. Erwähnt seien vor allem die 
Vorträge von 

E. Kamke (Tübingen): Zur Festlegung und 
Berechnung der Eigenwerte bei gewöhnlichen 
Differentialgleichungen durch die sog. direkten Me- 
thoden der Variationsrechnung, 

F. Rellich (Dresden): Vollständigkeit gestör- 
ter Eigenfunktionen, 

K. Schröder (Berlin): Untersuchungen über 
biharmonische Funktionen: Zur allgemeinen 
Theorie der ebenen und räumlichen biharmoni- 
schen Randwertaufgabe; die Randwertaufgabe für 
die rechteckige Kontur, 

J.L. Krames (Wien): Über die geometrischen 
Grundlagen der Luftphotogrammetrie. 

Am Mittwoch Nachmittag fand die Jahresver- 
sammlung des Mathematischen Reichsverbandes 
statt, in der insbesondere die Richtlinien für die 
Ausbildung der Diplom-Mathematiker besprochen 
wurden, über die Herr Oberregierungsrat Dr. Füh- 
rer berichtete. 320 


Persönliches. 

Der außerplanmäßige Prof. a. d. Univ. Göttingen 
Dr. Maximilian Schuler wurde zum a. o. Prof. 
für angewandte Mechanik ernannt. 

Der Dozent a. d. Techn. Hochschule Brünn 
Dr. W. v. Koppenfels wurde zum a. o. Prof. 
der Mathematik an dieser Hochschule ernannt. 

Am 26. April 1941 verstarb der o. Prof. für Strö- 
mungslehre und angewandie Mathematik, Direktor 
des Aerodynamischen Institutes der Techn. Hoch- 
schule Aachen C. Wieselsberger. 

Ferner verstarb in Berlin-Lichterfelde der frühere 
Öbservator und Prof. am Astronomischen Rechen- 
institut der Universität Berlin Dr. Jean Peters, 
der Herausgeber bekannter Tafelwerke. 
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